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ABSTRACT 



Studies of the structure of excited baryons are key to the N* program at Jefferson Lab. Within the first year 
of data taking with the Hall B CLAS12 detector following the 12 GeV upgrade, a dedicated experiment will 
aim to extract the N* electrocouplings at high photon virtualities Q 2 . This experiment will allow exploration 
of the structure of N* resonances at the highest photon virtualities ever yet achieved, with a kinematic reach up 
to Q 2 = 12 GeV 2 . This high-Q 2 reach will make it possible to probe the excited nucleon structures at distance 
scales ranging from where effective degrees of freedom, such as constituent quarks, are dominant through the 
transition to where nearly massless bare-quark degrees of freedom are relevant. In this document, we present a 
detailed description of the physics that can be addressed through iV* structure studies in exclusive meson elec- 
troproduction. The discussion includes recent advances in reaction theory for extracting N* electrocouplings 
from meson electroproduction off protons, along with QCD-based approaches to the theoretical interpretation 
of these fundamental quantities. This program will afford access to the dynamics of the non-perturbative strong 
interaction responsible for resonance formation, and will be crucial in understanding the nature of confinement 
and dynamical chiral symmetry breaking in baryons, and how excited nucleons emerge from QCD. 
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I. THE CASE FOR NUCLEON RESONANCE STRUCTURE STUDIES AT HIGH PHOTON VIRTUALITIES 

How virtual photons couple to ground state nucleons in forming excited nucleon states N*, in a distance- 
dependent way, opens the door to the most challenging sector of the Standard Model, that is, the strong interaction 
in the regime of large quark-gluon coupling or in the non-perturbative regime. Key to understanding the nature 
of the strong force in the non-perturbative regime is extracting the electromagnetic amplitudes for the transitions 
between ground and excited nucleon states, i.e. the j v NN* electrocouplings, as they evolve with photon virtuality 
Q 2 HH21- These electrocoupling studies are the necessary first steps in understanding how QCD generates most of 
the matter or mass in the real world, namely, mesons, baryons, and atomic nuclei. 

The non-perturbative strong interaction is enormously challenging. The degrees of freedom are not asymptotically- 
free current quarks and gauge gluons. The non-perturbative interaction of quarks and gluons is entirely different 
from that which exists within the perturbative Quantum ChromoDynamics (pQCD) realm and it gets quite com- 
plicated with all current quarks and gauge gluons becoming "dressed" by a cloud of virtual gluons and qq pairs. 
In the regime of large quark-gluon coupling, the dressing of current bare quarks by dressed gluons gives rise 
to a momentum-dependent dynamical mass and structure of dressed quarks; and these are the effective degrees 
of freedom employed in constituent quark models. In the evolution of the strong interaction from the pQCD 
regime of almost point-like and weakly coupled quarks and gluons (distances < 10~ 17 m) to the non-perturbative 
regime, where dressed quarks and gluons acquire dynamical mass and structure (distances w 10 -15 m), two major 
non-perturbative phenomena emerge: a) quark-gluon confinement and b) Dynamical Chiral Symmetry Breaking 
(DCSB). They both are completely outside of the pQCD scope. 

More than 98% of the hadron mass is generated non-perturbatively through DCSB processes, while the Higgs 
mechanism accounts for less than 2% of the light-quark baryon masses. The studies of the Dyson-Schwinger 
Equation of QCD (DSEQCD) have shown that dressing processes in the large quark-gluon running coupling 
regime are responsible for quark-gluon confinement. How confinement and DCSB emerge from QCD remains 
a challenging problem in present-day hadron physics, which is being addressed through studies of the "f v NN* 
electrocouplings. Extracting the j v NN* electrocouplings gives information on the dressed-quark mass, structure, 
and non-perturbative interaction, which is critical in exploring the nature of quark/gluon confinement and DCSB 
in baryons. The non-perturbative strong interaction is responsible for the formation of all individual N* states as 
bound systems of quarks and gluons. Experimental studies of the structure of all prominent N* -states, in terms of 
the ^ V NN* electrocoupling evolution with Q 2 carried out in close connection with QCD-based theories, offer a 
promising means of delineating the nature of the strong interaction in the non-perturbative regime. 

In particular, the extraction of ^ V NN* electrocouplings from the data on meson electroproduction off nucleons 
ll8T nTI serve to promote understanding the strong force. In May of 2012 the 6-GeV program with the CEBAF 
Large Acceptance Spectrometer (CLAS) in Hall B at Jefferson Lab was successfully completed. Among the many 
data runs with photons and electrons were several dedicated experiments focusing on hadron spectroscopy and 
hadron structure. CLAS was a unique instrument formed of a set of detectors and designed for the comprehensive 
exploration of exclusive meson electroproduction off nucleons. CLAS afforded excellent opportunities to study the 
electroexcitation of nucleon resonances in detail and with great precision. The CLAS detector has contributed the 
lion's share of the world's data on meson photo- and electroproduction in the resonance excitation region l9l [T0l[T2l - 
16 1 . For the first time detailed information from sets of unpolarized cross sections and different single- and double- 
polarization asymmetries have become available for many different meson photo- and electroproduction channels 
off protons and neutrons. 

The electroexcitation amplitudes for the low-lying resonances P33(1232), Pn(1440), P/i3(1520), and Sii(1535) 
were determined over a wide range of Q 2 in independent analyses of the n + n, 7r°p, i]p, and ir + ir~p electroproduc- 
tion channels Ifl4l [TTl [181 . Two of them, the P n (1440) and Z? 13 (1520) electrocouplings, have become available 
through independent analyses of single- and charged-double-pion electroproduction channels. The successful 
description of the measured observables in these exclusive channels resulting in the same j v NN* electrocou- 
pling values confirms their reliable extraction from the experimental data. These results have recently been 
complemented by still preliminary electrocouplings of high-lying resonances with masses above 1 .6 GeV |[T9l 
in the Tr + ir~p electroproduction channel, which is particularly sensitive to high-mass resonances. An alternative 
resonance electrocoupling extraction in a combined multi-channel analysis of the Nn, Nrj, and KY channels 
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within the framework of the advanced Excited Baryon Analysis Center Dynamic Coupled Channel (EBAC-DCC) 
approach is in progress [20, 21 1 and is further described in Chapter |ILD" 
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FIG. 1. (Left) The A1/2 electrocoupling of the Pn (1440) excited state from the analyses of the Nir electroproduction data 1171 
(circles), tt + tv~p electroproduction data 1 13] (triangles), and preliminary results from the iv + ir~p electroproduction data at Q 2 
from 0.5 to 1.5 GeV 2 1121 (squares). The photocouplings are taken from RPP |22| (open square) and the CLAS data analysis 
1231 (open triangle). Predictions from relativistic light front quark models [24 25 1 are shown by black solid and dashed lines, 
respectively. The absolute value of the meson-baryon cloud contribution as determined by the EBAC-DCC coupled-channel 
analysis |26| is shown by magenta thick solid line. (Right) The A1/2 electrocoupling of the Di3(1520) state. The data symbols 
are the same as in the left panel. The results of the hypercentral constituent quark model |27 1 and the absolute value of meson- 
baryon dressing amplitude 1 26] are presented by the black thin and magenta thick solid lines, respectively. 



The CLAS data on the P 33 (1232), Pn(1440), P>i 3 (1520), and Sii(1535) electrocouplings, when compared to 
the predictions of relativistic light front quark models l24l |25| and the results on the N* meson-baryon dressing 
amplitudes from the advanced EBAC-DCC coupled channel analysis H26I . shed additional light on the relevant 
components in the structure of N* -states at different distance scales. It was found that the structure of nucleon 
resonances in the mass range below 1.6 GeV is determined by contributions from both: a) an internal core of 
three dressed quarks and b) an external meson-baryon cloud. As an example, these contributions to the structure 
of the Pn(1440) and P/i 3 (1520) states are shown in Fig. |T| The absolute values of the meson-baryon-dressing 
amplitudes are maximal for Q 2 < 1.0 GeV 2 . They decrease with Q 2 and in the region of Q 2 > 1.0 GeV 2 , there is 
a gradual transition to where the quark degrees of freedom dominate, as is demonstrated by a better description of 
the P 11 (1440) and P 13 (1520)electrocouplings that are obtained within the framework of quark models. 

At photon virtualities of Q 2 > 5.0 GeV 2 , the quark degrees of freedom are expected to dominate the N* 
structure 0. This expectation is supported by the present analysis of the high-Q 2 behavior of the j v pN* elec- 
trocouplings 1T71 shown in Fig. |5j where the electrocoupling values scaled with Q 3 , as expected by constituent 
counting rules, are plotted. The indicated onset of scaling seen for Q 2 > 3.0 GeV 2 is likely related to the pref- 
erential interaction of the photon with dressed quarks, while interactions with the meson-baryon cloud results in 
strong deviations from this scaling behavior at smaller photon virtualities (greater distances). 

Therefore, in the "f v pN* electrocoupling studies for Q 2 > 5.0 GeV 2 , the quark degrees of freedom, as expressed 
by the N* structure, will be directly accessible from experiment with only small or even negligible contributions 
from the external meson-baryon cloud. This will mark the first time this new and fully unexplored region in the 
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FIG. 2. The A 1/2 electrocouplings of Pn (1440) (triangles), Di 3 (1520) (squares), 5x1(1535) (circles), and Fi 5 (1685) (stars) 
scaled with Q 3 from the CLAS data analysis 1171 . 




FIG. 3. (Left) The mass function for the s-quark evaluated within the framework of LQCD |28| (points with error bars) 
and DSEQCD [29, 30 1 (solid lines) for two values of bare masses, 70 MeV and 30 MeV, are shown in green and magenta, 
respectively. The chiral limit of zero bare quark mass, which is close to the bare masses of u and d quarks, is shown in red. 
Momenta p < 0.4 GeV correspond to the confinement, while those at p > 2. GeV correspond to the regime which is close to 
pQCD. The areas that are accessible for mapping of the dressed quark mass function by the "/ V NN* electrocoupling studies 
with 6 GeV and 1 1 GeV electron beams are shown to the left of blue solid and red dashed lines, respectively. (Right) Available 
(filled symbols) and projected CLAS12 |7| (open symbols) A 1 /2 electrocouplings of the Pn(1440) excited state. 
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electroexcitation of nucleon resonances can and will be investigated. A dedicated experiment on the N* studies in 
exclusive meson electroproduction off protons with the CLAS12 detector ( E12-09-003) is scheduled to take place 
within the first year after the completion of the JLab 12-GeV Upgrade Project [7| to Hall B. By measuring the 
differential cross sections off protons in the exclusive single-meson and double-pion electroproduction channels, 
complemented by single- and double-polarization asymmetries in single-meson electroproduction, this experiment 
seeks to obtain the world's only foreseen data on the electrocouplings of all prominent N* states in the still 
unexplored domain of photon virtualities up to 12 GeV 2 . As an example, the projected Ain electrocoupling values 
of Pn(1440) at photon virtualities from 5 to 12 GeV 2 are shown in the right panel of Fig. [5] A similar quality of 
results is expected for the electrocouplings of all other prominent N* states. The available reaction models for the 
extraction of the resonance electrocouplings have to be extended toward these high photon virtualities with the goal 
to reliably extract the ^f v pN* electrocouplings from the anticipated data on meson electroproduction off protons. 
In particular, the new reaction models have to account for a gradual transition from meson-baryon to quark degrees 
of freedom in the non-resonant reaction mechanisms. The current status and prospects of the reaction model 
developments are discussed in the Chapter|H| 

At present, there are two conceptually different approaches, which are used for interpreting the experimental 
results on resonance electrocouplings. Both these of approaches start from the first principles of QCD; they are: 
a) Lattice QCD (LQCD) and b) Dyson-Schwinger Equation of QCD (DSEQCD). Recent progress and the future 
prospects of these two approaches in N* structure studies are outlined in Chapters [Hi] to |VJ respectively. 

On the left hand side of Fig.[3]are plotted the dressed-quark masses as a function of momentum running over the 
quark propagator as calculated by DSEQCD |29 30 1 and LQCD |28 1 for different values of the bare-quark mass. 
The sharp increase from the mass of almost undressed current quarks (p > 2 GeV) to dressed constituent quarks 
(p < 0.4 GeV) clearly demonstrates that the dominant part of the dressed quark and consequently the hadron 
mass in toto is generated non-perturbatively by the strong interaction. The bulk of the dressed quark mass arises 
from a cloud of low-momentum gluons attaching themselves to the current-quark in the regime where the running 
quark-gluon coupling is large, and which is completely inaccessible by pQCD. The region where the dressed 
quark mass increases the most further represents the transition domain from pQCD (p > 2 GeV) to confinement 
(p < 0.4 GeV). The solution of the DSE gap-equation OT1 shows that the propagator pole in the confinement 
regime leaves the real-momentum axis and the momentum squared p 2 of the dressed quark becomes substantially 
different than that for the dynamical mass squared [M(p)] 2 . This means, that the dressed quark in the confinement 
regime will never be on-shell as it is required for a free particle when it propagates through space-time. Hence 
dressed quarks have to be strongly bound, locked inside the nucleon, and confinement becomes a property of the 
dressed quark and gluon propagators. These dressing mechanisms are responsible for DCSB. 

The DSEQCD studies 12 [31] have established that the quark-core contribution to the electromagnetic transi- 
tion amplitudes from the ground to excited nucleon states are determined primarily by the processes as depicted 
in Fig. [4] The momentum-dependent dressed-quark mass affects all quark propagators and the virtual photon in- 
teraction with the dressed-quark electromagnetic currents affords access to the dynamical quark structure. The 
Schwinger interaction of the virtual photon with the transition currents between diquark and two-quark states elu- 
cidates the very details of the strong interaction between and among dressed quarks. The value of momentum 
that is carried away by a single quark can roughly be estimated by assuming equal sharing of the virtual photon 
momentum among all three dressed quarks. Under this assumption it is reasonably straightforward to see that the 
measurements of ^ V NN* electrocouplings at 5 GeV 2 < Q 2 < 12 GeV 2 will be able to span nearly the entire 
range of quark momenta where the transition from confinement to pQCD occurs, as is seen Fig. [3] Therefore, the 
DSEQCD analyses on the -f v NN* electrocouplings of all prominent N* states expected from CLAS12 will offer 
a unique way to explore how quark-gluon confinement emerges from QCD and how more than 98% of each in- 
dividual nucleon resonance mass is generated non-perturbatively from the nearly massless current quarks through 
DCSB. 

Lattice QCD opens up an altogether different way for conceptualizing and interpreting the *y v NN* electro- 
couplings. Recent advances have shown the promising potential of LQCD in describing the resonance <y v NN* 
electrocouplings from first principles of QCD, in that LQCD starts from the QCD Lagrangian. Proof-of-principle 
results on the Q 2 -evolution of the F^ 11 (Q 2 ) and F^ 11 (Q 2 ) form factors for the transition from the ground state 
proton to the excited Pi i (1440) state have recently become available employing unquenched LQCD calculations 
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FIG. 4. The dressed quark interactions for the quark-core contribution to the electromagnetic transition amplitudes (<y v NN* 
electrocouplings) from the ground nucleon state of four-momentum Pi to excited nucleon states of four-momentum Pj in the 
DSEQCD approach 1 3 1 ] . Solid lines and double-solid lines stand for dressed quarks and the superposition of scalar and axial- 
vector diquark propagators, respectively. The T vertices describe the transition amplitudes between two-quark and diquark 
states, while the X -vertices represent the Schwinger interaction of the virtual photon with the transition current between the 
diquark and two-quark states. The ipi and ipf amplitudes describe the transitions between the intermediate diquark quark states 
and the initial nucleon and final N* -states, respectively. 



see Fig. 26 The corresponding experimental values of the Ff 11 (Q 2 ) and F^ 11 (Q 2 ) form factors are com- 
puted from the CLAS results ifTTl on the 7„pPn (1440) electrocouplings. 

Despite the simplified basis of the projection operators used in these computations, along with relatively large 
pion masses of s» 400 MeV, a reasonable description of the experimental data from CLAS IfTTl [181 was still 
achieved. In the future, when the LQCD evaluation of -f v NN* electrocouplings will become available, a realistic 
basis of the projection operators in a volume or box of relevant size with the physical pion mass, we will be able 
to compare the experimental electrocoupling results for all prominent N* states with LQCD. And the goodness 
of this comparison will allow us to answer the most challenging question in the Standard Model of whether QCD 
is in fact the fundamental theory of strong interactions and whether QCD is indeed sufficient in accounting for 
the full complexity of non-perturbative mechanisms, which are responsible for generating the ground and excited 
hadron states from the quarks and gluons degrees of freedom. Consistent results of those calculated by the two 
conceptually different frameworks of DSEQCD and LQCD in extracting the j v NN* electrocouplings, will offer 
compelling evidence for predicting the distance-dependent behavior of resonance electrocouplings. And consistent 
they should be as both are predicated on the first principles of QCD. 

In general, studies of the N* structure at high Q 2 address most of the fundamental issues of present-day hadron 
physics: 



1. How does nature achieve confinement? 
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2. How is confinement tied into dynamical chiral symmetry breaking, which describes the origin of more than 
98% of all visible mass in the universe? 

3. Can the full complexity of the N* structure be described based on the fundamental QCD Lagrangian? 

The current state of knowledge on the structure of nucleon ground and excited states makes for a strong moti- 
vation to systematically study the electrocouplings for all prominent excited baryons. These studies are necessary 
for fully accessing the complexity of the quark/gluon interactions. These very interactions are responsible for 
the formation of a given resonance, wherein a unique set of quantum numbers will characterize each individual 
resonance. 

DSEQCD studies have revealed the sensitivity of j v NN* electrocoupling to diquark correlations in baryons. 
This qq-paii correlation is generated by a non-perturbative strong interaction, which is responsible for meson 
formation, and can be described by the finite sizes of quasi-particles formed of paired quarks. In the DSEQCD 
approach, the two-quark assembly is described by either a superposition of pseudo-scalar and axial-vector diquarks, 
or a superposition of scalar and vector diquarks. It turns out that the relative contributions of the possible diquark 
components strongly depend on the quantum numbers of the N* state. Furthermore, the amplitudes shown in 
Fig. |4] which describe the transitions between the intermediate diquark-quark state and the initial ground ip^ or 
the final N* state ipf, are strongly dependent on the quantum numbers of both the initial nucleon and the final 
N* states. Once again, information on the electrocouplings of as many N* states as possible is needed for fully 
separating and identifying the mechanisms in the electroexcitation of nucleon resonances. And this information is 
further required for gaining access to the dressed-quark mass function and thereby the dynamical quark structure 
from data delineating the Q 2 evolution of resonance electrocouplings. 

Recent LQCD studies of the N* spectrum and structure 1 3 1 have further demonstrated the need for photo- and 
electrocoupling data over the the full range in the excited nucleon spectrum. Specifically, LQCD results of the 
full spectrum of excited proton state have recently become available. The excited proton state structure was de- 
termined in terms of contributions of three-quark configurations, which are described by vectors in the irreducible 
SU s /(6)xO(3) group representations . It has been found that the structure of N* states having masses less than 
1.75 GeV is dominated by no more than two SU s /(6)xO(3) configurations. However, the structure of several 
higher- lying nucleon excitations (M > 1.75 GeV) represent a superposition of many different configurations, 
which is another clear piece of theoretical evidence that both the low- and high lying-resonance electrocouplings 
must be measured. 

LQCD results 1 3 ] further predict the contributions for particular configurations in the resonance structure that 
should couple strongly to the glue. Such N* states, moreover, with these dominant contributions from such 
configurations will represent baryon hybrids. The observation of such hybrids would eliminate the distinction 
between gluon fields as the sole carriers of the strong interaction and quarks as the sole sources of massive matter. 
The proposed search for hybrid N* s [ 8 ] opens up yet another door in the N* program with the CL AS and CL AS 1 2 
detectors. Based on the LQCD results [3|, the hybrid N* masses are expected to be heavier than 1.9 GeV. The 
hybrid states may be seen as an overpopulation of the SU s /(6)xO(3) multiplets. However, the hybrid N*s should 
have the same quantum numbers as regular N*s, hence only through measuring the Q 2 evolution of the7„pA* 
electrocouplings, which is a reflection of a specific configuration of the baryon's structure, can the hybrid nature 
of the excited state be established. The high-Q 2 regime is of particular interest, since it is in this region where the 
contribution from quark and gluon degrees of freedom to the N* structure is expected to dominate. 

For all these reasons, data on electrocouplings for all prominent N* states are important for understanding how 
each individual N* structure emerges from QCD. As quantum numbers characterize the baryon resonances, we will 
separate each N* state, not only by mass, but by tagging its spin, angular momentum, radial, and parity. Moreover, 
since the nucleon structure for both the ground and excited states is generated by the non-perturbative quark- 
gluon interaction, future experiments on elastic form factors as well as Generalized Parton Distribution (GPD) 
and Transverse Momentum Distribution (TMD) structure functions will necessarily be incomplete if the studies 
are limited to the nucleon ground state structure. It requires studying the full complexity of the non-perturbative 
interactions as they generate nucleons from quarks and gluons. Therefore, all these experiments, combined, with 
coordinated studies on the ground and N* states will yield an understanding of confinement and DCSB in baryons. 

Results on j v pN* electrocouplings at high Q 2 will further allow access the parton distributions expressed in 
excited nucleon states within the framework discussed in the ChapterJV] Present-day knowledge on parton distri- 
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butions in baryons is limited to the ground state nucleons only. New information on the parton degrees of freedom 
in excited nucleon states will allow us to further develop the GPD concept, extending it to the transition between 
the ground and excited nucleon states and offering opportunities to map out the N — > N* transition densities in 
three-dimensional space. 

The need to analyze information on j v pN* electrocouplings for all prominent N*, requires extending the scope 
of theoretical analyses by making use of constituent quark models. Currently these models are the only available 
tool for physics analyses to link all the different resonance electrocouplings. Despite the shortcoming of lacking a 
connection to the fundamental QCD Lagrangian, quark models offer valuable information on resonance structure 
in the physics analysis of the j v pN* electrocoupling data. And such studies can guide the development of QCD- 
based approaches. The prospects for physics analyses of "f v pN* electrocouplings at high photon virtualities within 
the framework of quark models is discussed in the Chapter fVTl| 

A new theoretical tool for hadron physics, "Light-Front Holography", derived from mapping the dynamics of 
AdS-space to physical space-time at fixed light-front time r = t + z/c, has led to new insights into the color- 
confining, non-perturbative dynamics and the internal structure of relativistic light-hadron bound states. The 
AdS/QCD formalism is relativistic and frame-independent. Hadrons are described as eigenstates of a light-front 
Hamiltonian with a specific color-confining potential. A single parameter k sets the mass scale of the hadrons. 
The hadronic spectroscopy of the light-front holographic model gives a good description of the masses of the ob- 
served light-quark mesons and baryons. Elastic and transition form factors are computed from the overlap of the 
light-front wavefunctions. Many predictions of light-front holography for baryon resonances can be tested at the 



12 GeV JLab facility. An outline of this new method is given in section VII. D 

A strong collaboration between experimentalists and theorists is therefore required and indeed, has been estab- 
lished for achieving the challenging objectives in pursuing N* studies at high photon virtualities. Three topical 
Workshops 0314351 have been organized by Hall B, the Theory Center at Jefferson Lab, and the University of 
South Carolina to foster these efforts and create opportunities to facilitate and stimulate further growth in this field. 
This overview is prepared based on the presentations and discussions at these dedicated workshops with the goal 
to develop: 

1. reaction models for the extraction of the j v pN* electrocouplings from the data on single-meson and double- 
pion electroproduction off protons at photon virtualities from 5.0 to 12.0 GeV 2 by incorporating the transi- 
tion from meson-baryon to quark degrees of freedom into the reaction mechanisms; 

2. approaches for the theoretical interpretation of j v pN* electrocouplings, which are capable to explore how 
N* states are generated non-perturbatively by strong interaction and how these processes emerge from the 
QCD. 
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N*, A* 
states 


Branching fraction 

Ntt [%] 


Branching fraction 

iV7r7T [%] 


Prominent in 
Ntt exclusive 
channels 


Prominent in the 
■K + ir~p exclusive 
channel 


P 33 (1232) 


100 





* 




Pn(1440) 


60 


40 






Di 3 (1520) 


60 


40 






5n(1535) 


45 


< 10 


* 




531(1620) 


< 25 


75 




* 


5n(1650) 


75 


< 15 






Fi 5 (1680) 


65 


35 




* 


D 33 (1700) 


< 15 


85 






Pis(1720) 


< 15 


>70 




* 


F 35 (1905) 


< 10 


90 




* 


^37(1950) 


40 


>25 







TABLE I. Ntt and Nnir branching fractions for decays of excited proton states that have prominent contributions to the 
exclusive single- and/or charged-double-pion electroproduction channels. The values are taken from [22| or from the CLAS 
data analyses 1 171 1 1 91 . Symbols * mark most suitable exclusive channel(s) for the studies of particular N* state. 



II. ANALYSIS APPROACHES FOR EVALUATION OF NUCLEON RESONANCE ELECTROCOUPLINGS 
FROM THE CLAS DATA: STATUS AND PROSPECTS 

II.A. Introduction 

The CLAS detector at Jefferson Lab is a unique instrument, which has provided the lion's share of the world's 
data on meson photo- and electroproduction in the resonance excitation region. Cross sections and polarization 
asymmetries collected with the CLAS detector have made it possible for the first time to determine electrocouplings 
of all prominent N* states over a wide range of photon virtualities (Q 2 < 5.0 GeV 2 ) allowing for a comprehensive 
analysis of exclusive single-meson ( ir + n, ir°p, rjp, and KY) reactions in electroproduction off protons. Fur- 
thermore, CLAS was able to precisely measure 7r + 7r~p electroproduction differential cross sections owing to the 
nearly full kinematic coverage of the detector for charged particles. 

With the advent of the future CLAS 12 detector, the Q 2 reach will be considerably extended for exploring the 
nature of confinement and Dynamical Chiral Symmetry Breaking in baryons for our N* structure studies. Indeed, 
the CLAS 12 detector will be the sole instrument worldwide that will allow for performing experiments to determine 
the j v NN* electrocouplings of prominent excited proton states as listed in the Table [I] These will be the highest 
photon virtualities yet achieved for N* studies with photon virtualities in the range between 5.0 and 10 to 12.0 
GeV 2 , where the upper Q 2 boundary depends on the mass of excited proton state. The primary objective of 
the dedicated experiment E12-09-003, Nucleon Resonance Studies with CLASH Q, is to determine the ^ V NN* 
electrocouplings from the exclusive-meson electroproduction reactions, TT + n, tt°p, and tt + tt~p, off protons. The 
CLAS 12 experiment E12-09-003 Q was approved for 40 days of running time and it is scheduled to start in the 
first year of running with the CLAS 12 detector, that is right after the Hall B 12 GeV upgrade. This experiment 
represents the next step towards extending our current N* Program with the CLAS detector |9][l0l; it will make 
use of the 1 1-GeV continuous electron beam that will be delivered to Hall B of Jefferson Lab. We remark that the 
maximum energy of the electron beam to Hall B will be 1 1 GeV for the 12 GeV upgrade to JLab. 

The data from the ir + n, ir°p, and tt + tt~p electroproduction channels will play a key role in the evaluation 
of "f v NN* electrocouplings. The two primary Ntt and tt + tt~p exclusive channels, combined, account for ap- 
proximately 90% of the total cross section for meson electroproduction in the resonance excitation region of 
W < 2.0 GeV. Both single- and charged-double-pion electroproduction channels are sensitive to N* contribu- 
tions, as can be seen in Table [I] Further, these two channels offer complementary information for cross checking 
the N* parameters derived in the fits to the observables in the exclusive channel. 

A necessary first step towards extracting the *y v NN* electrocouplings, in a robust way, requires that we em- 
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ploy independent analyses of the single- and charged-double-pion electroproduction data within the framework of 
several different phenomenological reaction models. A reliable separation of resonant and non-resonant contri- 
butions, moreover, is crucial for evaluating the *y v NN* electrocouplings within each of the frameworks provided 
by these approaches. Independent analyses of major meson electroproduction channels offer a sensitive test as 
well as a quality check in separating the contributions of resonant and non-resonant mechanisms. Consistent ex- 
tractions of the N* parameters among different channels are imperative. The ir + n, ir°p, and ir + ir~p exclusive 
electroproduction channels, for example, have entirely different non-resonant mechanisms. Clearly, the value of 
the ^ V NN* electrocouplings must be analysis independent and must remain the same in all exclusive channels, 
since resonance electroexcitation amplitudes do not depend on the final states that will populate the different ex- 
clusive reaction channels. Therefore, consistency in ascertaining the values of the <y v NN* electrocouplings from 
a large body of observables, as measured in 7r + n, ir°p and ir + ir~p electroproduction reactions, will give good 
measure of the reliability of the extraction of these fundamental quantities. In the next section we shall review 
the current status and prospects for developing several phenomenological reaction models with the primary ob- 
jective of determining ^y v NN* electrocouplings from independent analyses of the single- and charge-double -pion 
electroproduction data (7r + n, ir°p, and ir + Tr~p). 



II.B. Approaches for independent analyses of the CLAS data on single- and charged-double-pion electroproduction 

off protons 

Several phenomenological reaction models lfT71 - [T9l I36H381 were developed by the CLAS Collaboration for 
evaluating the j v NN* electrocouplings in independent analyses of the data on n + n, ir°p, and ir + ir~p electropro- 
duction off protons. These models were successfully applied to single-pion electroproduction for Q 2 < 5.0 GeV 2 
and W < 1.7 GeV. The ^ V NN* electrocouplings were extracted from the CLAS Tr + ir~p electroproduction data 
for the kinematical ranges of Q 2 < 1.5 GeV 2 and W < 1.8 GeV EIUED. ^ CLAS data on sin gle-pion 
exclusive electroproduction were also analyzed within the framework of the MAID [ 1 1| and the SAID [39 40) ap- 
proaches. These reaction models have allowed us to access resonant amplitudes by fitting all available observables 
in each channel independently and within the framework of different reaction models. Consequently, the "f v NN* 
electrocouplings, along with the respective Nir and iV7r7r hadronic decay widths, were determined by employing 
a Breit-Wigner parameterization of the resonant amplitudes. 



II.B.l. CLAS Collaboration approaches for resonance electrocoupling extraction from the data on single-pion 

electroproduction off protons 

Analyses of the rich CLAS data samples have extended our knowledge considerably of single-meson electro- 
production reactions off protons, i.e. the ir + n and ir°p channels. A total of nearly 120,000 data points have been 
collected on reactions arising from unpolarized differential cross sections, longitudinally-polarized beam asymme- 
tries as well as from longitudinal-target and beam-target asymmetries with a nearly complete coverage in the phase 
space for exclusive reactions ifTTl . The data were analyzed within the framework of two conceptually different ap- 
proaches, namely: a) the unitary isobar model (UIM) and b) a model, employing dispersion relations [36, 37). All 
well-established N* states in the mass range Mjy* < 1.8 GeV were incorporated into the Nir channel analyses. 

The UIM follows the approach detailed in Ref. [11 1. The Nir electroproduction amplitudes are described as a 
superposition of N* electroexcitations in the s-channel and non-resonant Born terms. A Breit-Wigner ansatz, with 
energy-dependent hadronic decay widths BT1I . is employed for the resonant amplitudes. Non-resonant amplitudes 
are described by a gauge-invariant superposition of nucleon s- and u-channel exchanges and in the t -channel by 
7r, p, and co exchanges. The latter are reggeized; this allows for a better description of the data in the second- 
and the third-resonance regions, whereas for W < 1.4 GeV, the role of Regge-trajectory exchanges becomes 
insignificant. The Regge-pole amplitudes were constructed using the prescription delineated in Refs. 11421 l43l 
allowing us to preserve gauge invariance of the non-resonant amplitudes. The final-state interactions are treated as 
irN rescattering in the K-matrix approximation [36). 
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In another approach, dispersion relations relate the real and imaginary parts of the invariant amplitudes, which 
describe Ntt electroproduction in a model-independent way ll36l . For the 18 independent invariant amplitudes 
describing the j v N — > Ntt transition electromagnetic current, dispersion relations (17 unsubtracted and 1 sub- 
tracted) at fixed t are employed. This analysis has shown that the imaginary parts of amplitudes are dominated by 
resonant contributions for W > 1.3 GeV. That is to say, in this kinematical region, they are described solely by 
resonant contributions. For smaller W values, both the resonant and non-resonant contributions to the imaginary 
part of the amplitudes are taken into account based on an analysis of irN elastic scattering and by making use of 
the Watson theorem and the requisite dispersion relations. 

For either of the these approaches, the Q 2 evolution of the non-resonant amplitudes is determined by the be- 
havior of the hadron electromagnetic form factors, which are probed at different photon virtualities. The s- and 
it-channel nucleon exchange amplitudes depend on the proton and neutron electromagnetic form factors, respec- 
tively. The i-channel tt, p, oj exchanges depend on pion electromagnetic form factors and p(uj) — > irj transition 
form factors. The exact parameterization of these electromagnetic form factors as a function of Q 2 that are em- 
ployed in the analyses of the CLAS single-pion electroproduction data can be found in Ref. [ 17 1. These analyses 
have demonstrated that for photon virtualities Q 2 > 0.9 GeV 2 , the reggeization of the Born amplitudes becomes 
insignificant in the resonance region for W < 1.9 GeV. Consequently, at these photon virtualities, the background 
of UTM was constructed from the nucleon exchanges in the s- and u-channels and in the i-channel through tt, p 
and uj exchanges. In the approach based on dispersion relations, we additionally take into account the Q 2 evolution 
of the subtraction function f su h(Q 2 , t). The subtraction function was determined using a linear parameterization 
over the Mandelstam variable t and through fitting two parameters to the data in each bin of Q 2 IfTTI . Employ- 
ing information on the Q 2 evolution of hadron electromagnetic form factors from other experiments or from our 
fits to the CLAS data, we are able to predict the Q 2 evolution of the non-resonant contributions to single-pion 
electroproduction in the region of Q 2 , where the meson-baryon degrees of freedom remain relevant. 
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These two approaches provide a good description of the Nit exclusive channel observables in the entire range 
covered by the CLAS measurements: W < 1.7 GeV and Q 2 < 5.0 GeV 2 , resulting in x 2 /d.p. < 2.0 for 
Q 2 < 1.0 GeV 2 and x 2 /d.p. < 3.0 at Q 2 from 1.5 to 4.5 GeV 2 flT). Exclusive structure functions a T +ea L , a T T, 
<Jlt, and <Jlt' were derived from the measured CLAS cross sections and polarization asymmetries. An example 
of the description of the structure function moments is shown in Fig. [5] The results from these two approaches 
further provide information for setting the systematical uncertainties associated with the models. And a consistent 
description of a large body of observables in the Nit exclusive channels, obtained within the respective frameworks 
of two conceptually different approaches, lends credibility to a correct evaluation of the resonance contributions. 

II.B.2. Evaluation ofy v NN* resonance electrocouplings from the data on charged-double-pion electroproduction off 

protons 

The it + tt~p electroproduction data measured with the CLAS detector \12, 13 1 provide information on nine 
independent one-fold differential and fully-integrated cross sections in a mass range of W < 2.0 GeV and for 
photon virtualities in the range of 0.25 < Q 2 < 1.5 GeV 2 . Examples of the available ir + ir~p one-fold differential 
cross-section data for specific bins in W and Q 2 are shown in Figs. [6] and [7] Analysis of these data have allowed 
us to establish which mechanisms contribute to the measured cross sections. The peaks in the invariant mass dis- 
tributions provide evidence for the presence of the channels arising from j v p —> Meson + Baryon —> n + n~p 
having an unstable baryon or meson in the intermediate state. Pronounced dependences in angular distributions 
further allow us to establish the relevant i-, u-, and s-channel exchanges. The mechanisms without pronounced 
kinematical dependences are identified through examination in various differential cross sections, with their pres- 
ence emerging from correlation patterns. The phenomenological reaction model JM lTT~8l [191 138 45 1 was developed 
in collaboration between Hall B at Jefferson Lab and the Skobeltsyn Nuclear Physics Institute in Moscow State 
University. The primary objective of this work is to determine the j v NN* electrocouplings, and corresponding 
7rA and pp partial hadronic decay widths from fitting all measured observables in the tt + tt~p electroproduction 
channel. 

The amplitudes of the j v p —> Tr + ir~p reaction are described in the JM model as a superposition of the 7r~ A ++ , 
7T + A°, pp, 7r + Z3j 3 (1520), 7r + F{ ) 5 (1685), and 7r~P^ + (1600) sub-channels with subsequent decays of unstable 
hadrons to the final ir + ir~p state, and additional direct 2ir -production mechanisms, where the final ir + ir~p state 
comes about without going through the intermediate process of forming unstable hadron states. 

The JM model incorporates contributions from all well-established N* states to the irA and pp sub-channels 
only. We also have included the 3/2 + (1720) candidate state, whose existence is suggested in the analysis lfl2ll 
of the CLAS ir + ir~p electroproduction data. In the current 2012 JM model version [ 18], the resonant amplitudes 
are described by a unitarized Breit-Wigner ansatz as proposed in Ref. |46l : it was modified to make it fully con- 
sistent with the parameterization of individual N* state contributions by a relativistic Breit-Wigner ansatz with 
energy-dependent hadronic decay widths [47 1 employed in the JM model. After unitarization, the Breit-Wigner 
ansatz accounts for transitions between the same and different N* states reflected in the dressed-resonant propa- 
gators, making resonant amplitudes consistent with the restrictions imposed by unitarity fT8l . Quantum number 
conservation in strong interactions allows the transitions between Z?i3(1520)/Z?i3(1700), S l n(1535)/S , n(1650), 
and 3/2+ ( 1 720) /Pi 3 ( 1720) pairs of N* states incorporated into the JM model and are listed in the Table [j] We 
found that use of the unitarized Breit-Wigner ansatz has a minor influence on the j v NN* electrocouplings, but it 
may substantially affect the N* hadronic decay widths determined from fits to the CLAS data. 

Non-resonant contributions to the 7rA sub-channels incorporate a minimal set of current-conserving Born terms 
ll38l l47l . They consist of t-channel pion exchange, s-channel nucleon exchange, u-channel A exchange, and 
contact terms. Non-resonant Born terms were reggeized to preserve current conservation, as proposed in Refs. 11421 
|43| . The initial- and final-state interactions in ttA electroproduction are treated in an absorptive approximation, 
with the absorptive coefficients estimated from the data from irN scattering ll47l . Non-resonant contributions to 
the 7r A sub-channels further include additional contact terms that have different Lorentz-invariant structures with 
respect to the contact terms in the sets of Born terms. These extra contact terms effectively account for non- 
resonant processes in the irA sub-channels beyond the Born terms, as well as for the final-state interaction effects 
that are beyond those taken into account by absorptive approximation. Parameterizations of the extra contact terms 
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FIG. 6. Fits to the CLAS ep -> e'vr+Tr-p data [T2) within the framework of JM model [19]|38) at W = 1.71 GeV and Q 2 =0.65 
GeV 2 . Full model results are shown by thick solid lines together with the contributions from 7r~A ++ (dashed thick lines), 
pp (dotted thick lines), 7r+A° (dash-dotted thick lines), 7r + D 13 (1520) (thin solid lines), and 7r + f 1 ° 5 (1685) (dash-dotted thin 
lines) isobar channels. The contributions from other mechanisms described in the Section [Tl.B.2| are comparable with the data 
error bars and they are not shown in the plot. 



in the irA sub-channels are given in Ref. [38). 

Non-resonant amplitudes in the pp sub-channel are described within the framework of a diffractive approxima- 
tion, which also takes into account the effects caused by p-line shrinkage [48 1. The latter effects play a significant 
role in the N* excitation region, and in particular, in near-threshold and sub-threshold p-meson production for 
W < 1.8 GeV. Even in this kinematic regime, when the non-resonant parts of the pp sub-channel become small, 
the pp sub-channel may affect the one-fold differential cross sections due to the contributions from nucleon reso- 
nances that decay into the pp final states. Therefore, a reliable and credible treatment of non-resonant contributions 
in the pp sub-channel becomes important for ascertaining the electrocouplings and corresponding hadronic param- 
eters of these resonances. The analysis of the CLAS data [ 12lH3ll has revealed the presence of the pp sub-channel 
contributions for W > 1.5 GeV. 



16 



W=1.71 GeV, Q^O.95 GeV 2 




1 1.5 

7T+ p mass, GeV 



0.1 0.3931 0.6862 0.9793 

tt+ 7T— mass, GeV 



1.5 

7T- p mass, GeV 




re — ,deg b 



20 - 



10 




1? 7T+ ,deg b 



100 200 

•# p.deg 




Pill! 
i i i 

_l I I I I I L 



200 

a(n — p)(n+ pf),deg 



200 

a (n+ p)(tt— pf),deg 



200 

a (pf p)(7T— TT + ),deg 



FIG. 7. Resonant (blue bars) and non-resonant (green bars) contributions to differential cross sections obtained from the CLAS 
data 1 12| fit within the framework of the JM model at W= 1.71 GeV, Q 2 =0.95 GeV 2 . Red lines show the fit results. 



The ir+Dl 3 (152Q), 7T + Ff 5 (1685), and tt~ P^ 3 + (1600) sub-channels are described in the JM model by non- 
resonant contributions only. The amplitudes of the 7r + _D5 , 3 (1520) sub-channel were derived from the non-resonant 
Born terms in the irA sub-channels by implementing an additional 75 -matrix that accounts for the opposite parities 
of the A and Z?i3(1520) [49 1 . The magnitudes of the 7r + Z)!j ) 3 (1520) production amplitudes were independently fit 
to the data for each bin in W and Q 2 . The contributions from the ir + Z?j , 3 (1520) sub-channel should be taken into 
account for W > 1.5 GeV. 

The 7r + -F 1 ° 5 (1685) and tt-P 33 + (1600) sub-channel contributions are seen in the data [12] at W > 1.6 GeV. 
These contributions are almost negligible at smaller W. The effective contact terms were employed in the JM 
model for parameterization of these sub-channel amplitudes l45l |49l . Magnitudes of the 7r + F° 5 (1685) and 
7r _ P 3 t 3 + (1600) sub-channel amplitudes were fit to the data for each bin in W and Q 2 . 

A general unitarity condition for ir + ir~p electroproduction amplitudes requires the presence of so-called direct- 
2n -production mechanisms, where the final ir + ir~p state is created without going through the intermediate step 
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of forming unstable hadron states 1501 . These intermediate-stage processes are beyond those aforementioned con- 
tributions from two-body sub-channels. Direct 2n production amplitudes were established for the first time in the 
analysis of the CLAS 7r + 7r~p electroproduction data [ 3 8 J . They are described in the JM model by a sequence of 
two exchanges in t- and/or u- channels by unspecified particles. The amplitudes of the 27r-production mechanisms 
are parameterized by an Lorentz-invariant contraction between spin-tensors of the initial and final-state particles, 
while two exponential propagators describe the above-mentioned exchanges by unspecified particles. The magni- 
tudes of these amplitudes are fit to the data for each bin in W and Q 2 . Recent studies of the correlations between 
the final-hadron angular distributions have allowed us to establish the phases of the 2tt direct-production ampli- 
tudes f5T). The contributions from the 2ir direct-production mechanisms are maximal and substantial (rj 30% ) for 
W < 1.5 GeV and they decrease with increasing W, contributing less than 10% for W > 1.7 GeV. However, even 
in this kinematical regime, 2ir direct-production mechanisms can be seen in the ir + ir~p electroproduction cross 
sections due to an interference of the amplitudes from two-body sub-channels. 

The JM model provides a reasonable description of the tt + tt~p differential cross sections for W < 1.8 GeV 
and Q 2 < 1.5 GeV 2 with a % 2 /d.p. < 3.0, accounting only for statistical uncertainties in the experimental data. 
As a typical example, the nine one-fold differential cross sections for W — 1.71 GeV and Q 2 = 0.65 GeV 2 , with 
fits, are shown in Fig. |6j together with the contributions from each of the individual mechanisms incorporated into 
the JM description. Each contributing mechanism has a distinctive shape for the cross section as is depicted by 
the observables in Fig. [6] Furthermore, any contributing mechanism will be manifested by substantially different 
shapes in the cross sections for the observables, all of which are highly correlated through the underlying-reaction 
dynamics. The fit takes into account all of the nine one-fold differential cross sections simultaneously and allows 
for identifying the essential mechanisms contributing to ir + Tr~p electroproduction off protons. Such a global fit 
serves towards understanding the underlying mechanisms and thereby affording access to the dynamics. 

This successful fit to the CLAS ir + Tr~p electroproduction data has further allowed us to determine the reso- 
nant parts of cross sections. An example is shown in Fig. [7] The uncertainties associated with the resonant part 
are comparable with those of the experimental data. It therefore provides strong evidence for an unambiguous 
separation of resonant/non-resonant contributions. A credible means for separating resonances from background 
was achieved by fitting CLAS data within the framework of the JM model and it is of particular importance in the 
extraction of the j v NN* electrocouplings, as well as for evaluating each of the excited states decay widths into the 
7rA and pp channels. A special fitting procedure for the extraction of resonance electrocouplings with the full and 
partial irA and pp hadronic decay widths was developed, thereby allowing us to obtain uncertainties of resonance 
parameters and which account for both experimental data uncertainties and for the systematical uncertainties from 
the JM reaction model [18]. 



II.C. Resonance electrocouplings from the CLAS pion electroproduction data 



Several analyses of CLAS data were carried out on single- and charged-double-pion electroproduction off pro- 
tons within the framework of fixed-i dispersion relations, the UIM model, and the JM model as described in Sec- 
which have provided, for the first time, information on electrocouplings of the Pi i (1440), 



II.B.l 



and II.B.2 



tions 

Z?i3(1520), and Fi5(1685) resonances from independent analyses of n + n, ir°p, and ir + ir~p electroproduction 
channels lfT7TIT9l . The electrocouplings of the i-ii(1440) and Z?i3(1520) resonances determined from these chan- 
nels are shown in Figs. [8] and [9] They are consistent within uncertainties. The longitudinal Si/ 2 electrocouplings 
of the ,013(1520), S , ii(I535)/75 3 i(1620), Su(1650), Fi 5 (1685), £> 33 (1700), and Pi 3 (1720) excited proton states 
have become available from the CLAS data for the first time as well ifTTl [191 . 

Consistent results on <y v NN* electrocouplings from the P 11 (1440), Z?i3(1520), and F 15 (1685) resonances 
that were determined from independent analyses of the major meson electroproduction channels, 7r + n, ir°p, and 
7r + 7r~p, demonstrate that the extraction of these fundamental quantities are reliable as these different electropro- 
duction channels have quite different backgrounds. Furthermore, this consistency also strongly suggests that the 
reaction models described in sections II. B. 1 and II.B.2 will provide a reliable evaluation of the j v NN* electrocou- 
plings for analyzing either single- or charged-double-pion electroproduction data. It therefore makes it possible to 
determine electrocouplings for all resonances that decay preferentially to the Nn and/or Ntttt final states. 
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FIG. 8. Ax/2 (left) and S1/2 (right) electrocouplings of the Pn(1440) resonance determined in independent analyses of the 
CLAS data on Ntt (circles) [171 , and tt + tt~p (triangles) |19| electroproduction off protons. Squares and triangles at Q 2 =0 
GeV 2 correspond to [22 1 and the CLAS Ntt [ 23] photoproduction results, respectively. 
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FIG. 9. Ax/2 (left), Si/2 (middle), and A- A / 2 (right) electrocouplings of the Di3(1520) resonance determined in independent 
analyses of the CLAS data on Ntt (circles) 1171 . and tt + tt~p (triangles) |19| electroproduction off protons. Squares and 
triangles at Q 2 =0 GeV 2 correspond to [22 ] and the CLAS Ntt [23 1 photoproduction results, respectively. 



The studies of Ntt exclusive channels are the primary source of information on electrocouplings of the N* 
states with masses below 1.6 GeV fl7l . The reaction kinematics restrict the P33(1232) state to only the Ntt 
exclusive channels. The Pn(1440) and Di3(1520) resonances have contributions to both single- and double-pion 
electroproduction channels, which are sufficient for the extraction of their respective electrocouplings. Analysis of 
the 7r + 7r~p electroproduction off protons allows us to check the results of Ntt exclusive channels for the resonances 
that have substantial decays to both the Ntt and A^Tr channels. 

For the Sii(1535) resonance, the hadronic decays to the final state is unlikely (see Table [jj. Therefore, the 
studies of this very pronounced A^7r-electroproduction resonance become problematic in the charged-double-pion 
electroproduction off protons. On the other hand, the Sii(1535) resonance has a large branching ratio to the ttN 
and r]N channels. And sincel999, this resonance has been extensively studied at JLab over a wide range of Q 2 
up to 4.5 and 7 GeV 2 for the channels A^tt and N-q, respectively in the electroproduction off protons (see Fig. 10 1. 
For Ntj electroproduction, the Sii(1535) strongly dominates the cross section for W < 1.6 GeV and is extracted 
from the data in a nearly model-independent way using a Breit-Wigner form for the resonance contribution |[T4l 
These analyses assume that the longitudinal contribution is small enough to have a negligible effect 
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FIG. 10. Transverse electrocoupling A1/2 of the 7*p — > Sn(1535) transition. The full circles are the electrocouplings 
extracted from Nir electroproduction data 1 17|. The electrocouplings extracted from Nrj electroproduction data are: the stars 
1521 . the open boxes 1141 . the open circles 1 53 1 , the crosses |54|, and the rhombuses 1361 1371 . The full box and triangle at 
Q 2 = correspond to |22 | and the CLAS Nn 1231 photoproduction results, respectively. 



on the extraction of the transverse amplitude. This assumption is confirmed by the analyses of the CLAS Nir 
electroproduction data IfTTl . Accurate results were obtained in both reactions for the transverse electrocoupling 
Ax/2', they show a consistent Q 2 slope and allowed for the determination of the branching ratios to the Nir and 
Nrj channels IfTTl . Transverse A1/2 electrocouplings of the 5n(1535) extracted in independent analyses of Nrr 
and Nrj electroproduction channels are in a reasonable agreement, after taking into the systematical uncertainties 
of the analysis IfTTl into consideration. Expanding the proposal, Nucleon Resonance Studies with CLAS 12 JTJ, by 
further incorporating Nr] electroproduction at high Q 2 would considerably enhance our capabilities for extracting 
self-consistent and reliable results for the 5*11(1535) electrocouplings in independent analyses of the Nir and Nr] 
electroproduction channels. 

The charged-double-pion electroproduction channel is of particular importance for evaluation of high-lying res- 
onance electrocouplings, since most N* states with masses above 1.6 GeV decay preferentially by two pion emis- 
sion (TablefTji. Preliminary results on the electrocouplings of the 531(1620), 5n(1650), Fi 5 (1685), D 33 (1700), 
and Pi 3 (1720) resonances were obtained from an analysis of the CLAS ir + Tr~p electroproduction data 1 12 1 within 
the framework of the JM model |l9l . As an example, electrocouplings of the -P33(1700) resonance were deter- 



mined from analysis of the CLAS ir + ir~p electroproduction data and are shown in Fig. 11 in comparison with the 
previous world's data taken from Ref. Il55l . The £>33(1700) resonance decays preferentially to A^7r7r final states 
with the branching fraction exceeding 80%. Consequently, electrocouplings of this resonance determined from the 
A^7r electroproduction channels have large uncertainties due to insufficient sensitivity of these exclusive channels 
contributing to the Z?33(1700) resonance. The CLAS results have considerably improved our knowledge on elec- 
trocouplings of the 53i(1620), 5n(1650), Fi B (1685), Z) 33 (1700), and P 13 (1720) resonances. They have provided 
accurate information on the Q 2 evolution of the transverse electrocouplings, while longitudinal electrocouplings 
of these states were determined, again, for the first time. 

Most of the Af* states with masses above 1.6 GeV decay preferentially through channels with two pions in the 
final state, thus making it difficult to explore these states in single-pion electroproduction channels. The CLAS 
KY electroproduction data lfT31 fT6ll may potentially provide independent information on the electrocouplings of 
these states. At the time of this writing, however, reliable information on KY hadronic decays from A^*s are not 
yet available. The N* hyperon decays can be obtained from fits to the CLAS KY electroproduction data lfT5lfT6l . 
which should be earned out independently in different bins of Q 2 by utilizing the Q 2 -independent behavior of 
resonance hadronic decays. The development of reaction models for the extraction of -f v NN* electrocouplings 
from the KY electroproduction channels is urgently needed. Furthermore, complementary studies of the KY 
decay mode can be carried out with future data from the Japan Proton Accelerator Research Complex (J-PARC) 
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FIG. 11. Electrocouplings of Z?33(1700) resonance Awj (left), Sx/j (middle) and ^3/2 (right) determined in analyses the 
CLAS 7r + 7r~p electroproduction data | 38| and world data on Nir electroproduction off protons 1551 . 



and through J/ty decays to various NN* channels at the Beijing Electron Positron Collider (BEPC). 

Most of the well-established resonances have substantial decays to either the Nir or Nirir final states. There- 
fore, studies of A^7r and ir + ir~p electroproduction off protons will allow us to determine the electrocouplings of 
all prominent excited proton states and such studies will mark the first step in the evaluation of resonance electro- 
couplings in the unexplored regime of photon virtualities ranging from 5 to 12 GeV 2 . 



II.D. Status and Prospect of Excited Baryon Analysis Center (EBAC) 

II.D.l. The case for a multi-channel global analyses 

Interactions among different hadronic final states are termed final-state interactions (FSI). In exclusive-meson 
electroproduction, for example, FSI represent a key issue both in terms of the extraction as well as in the physical 
interpretation of the nucleon resonance parameters. In the reaction models for analyses of different exclusive- 
meson electroproduction channels, as detailed above, FSI are treated phenomenologically for each specific reac- 
tion. Analyses of exclusive hadroproduction have allowed us to establish explicitly the relevant mechanisms for 
hadron interactions among the various final states for different exclusive photo- and electroproduction channels in 
terms of the meson-baryon degrees of freedom. The information on meson electroproduction amplitudes comes 
mostly from CLAS experiments. These results on meson-baryon hadron interaction amplitudes open up addi- 
tional opportunities for the extraction of resonances, their photo- and electrocouplings, as well as their associated 
hadronic decay parameters. These parameters can be constrained through a global analysis of all exclusive-meson 
electroproduction data from different photo- and electroproduction channels as analyzed within the framework of 
coupled-channel approaches 156-381. 

These approaches have allowed us to explicitly take into account the hadronic final-state interactions among 
the exclusive meson electroproduction channels and to build up reaction amplitudes consistent with the restric- 
tions imposed by a general unitarity condition. Another profound consequence of unitarity is reflected by the 
relations among the non-resonant meson production mechanisms and the contributions from meson-baryon dress- 
ing amplitudes (i.e. the meson-baryon cloud) to the resonance electrocouplings along with their hadronic decay 
parameters. Use of coupled-channel approaches have allowed us to determine such contributions in the fitting 
to the experimental data. Therefore, global analyses of all exclusive meson photo- and electroproduction data 
within the framework of coupled-channel approaches will reveal information on the resonance structure in terms 
of quark-core and meson-baryon cloud contributions at different distance scales. 

The Nir and ir + ir~p electroproduction channels are strongly coupled through final-state interactions. The data 
from experiments with hadronic probes have shown that the ttN — > irirN reactions are the second biggest exclusive 
contributors to inclusive irN interactions. Therefore, data on the mechanisms contributing to single- and charged- 
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double-pion electroproduction off protons are needed for the development of global multi-channel analyses for 
the extraction of j v NN* electrocouplings within the framework of coupled-channel approaches. A consistent de- 
scription of hadronic interactions between the irN and ttttN asymptotic states is critical for the reliable extraction 
of j v NN* electrocouplings within the framework of coupled-channel approaches. 

II.E. Dynamical Coupled Channel Model 

In this section, we report on the development and results of the EBAC-DCC approach spanning the period from 
January 2006 through March 2012. This analysis project has three primary components: 

1. perform a dynamical coupled-channels analysis on the world data on meson production reactions from the 
nucleon to determine the meson-baryon partial-wave amplitudes, 

2. extract the N* parameters from the determined partial-wave amplitudes, 

3. investigate the interpretations of the extracted N* properties in terms of the available hadron models and 
Lattice QCD. 

The Excited Baryon Analysis Center (EBAC) is conducting dynamical coupled-channel (DCC) analyses of 
Jefferson Lab data and other relevant data in order to extract N* parameters and to investigate the reaction mecha- 
nisms for mapping out the important components of the A* structure as a function of distance or Q 2 . This work is 
predicated upon the dynamical model for the A(1232) resonance [59], which was developed by the Argonne Na- 
tional Laboratory-Osaka University (ANL-Osaka) collaboration |60|. In the EBAC extension to the ANL-Osaka 
formulation |59|, the reaction amplitudes T a ^(p,p'; E) for each partial wave are calculated from the following 
coupled-channels integral equations, 

q 2 dqV an (p,q)G 7 {q, E)T J: p(q,p' , E) , (1) 
V a , = v a ,p + ^ E _ Mif , (2) 

N" 

where a, (3, 7 = 7JV, irN, r/N, KY, u)N, and ttttN which has irA, pN, aN resonant components, v a< p are meson- 
exchange interactions deduced from the phenomenological Lagrangian, rjv*,/3 describes the excitation of the nu- 
cleon to a bare N* state with a mass M* , and G 7 (q, E) is a meson-baryon propagator. The DCC model, defined 
by Eqs. ([T]l and Q, respects unitarity for both two- and three-body reactions. 

This dynamical coupled-channel model was used initially in fitting ttN reactions from elastic scattering to 
extract parameters associated with the strong-interaction parts of V a ,p, in Eq. (|2jl and corresponding electromagnetic 
components of V a ,p came from fits to the •yp — > tt°p, tt + ti and p(e, e'ir°' + )N data in the invariant mass range of 
W < 2 GeV. To simplify the analysis during the developmental stage (2006-2010), the KY and ujN channels 
were not included in these fits. 

The resulting six-channel model was then tested by comparing the the predicted ttN, 7 AT — > irirN production 
cross sections with data. In parallel to analyzing the data, a procedure to analytically continue Eqs. ([T} and (|2]) 
to the complex-energy plane was developed to allow for extracting the positions and residues of several nucleon 
resonances. In the following, we present a sample of some of our results in these efforts. 

ILEA. Results for single-pion production reactions 

In fitting the ttN elastic-scattering channel, we found that one or two bare N* states were needed for each partial 
wave. The coupling strengths of the N* — > MB vertex interactions T^* b with M B — irN, 77 A, it A, pN, aN 
were then determined in the x 2- fits to the data and these results can be found in Ref. 16TI . 



T a , p {p,p';E) =V a ,p(p,p') + y2 / 
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FIG. 12. The DCC results |26| of total cross sections (a), differential cross sections (da/dQ,), and photon asymmetry (E) of 
7P — > ir°p (upper parts), jp — > 7r + n (lower parts). 



Our next step was to determine the bare 7iV — > N* interaction Tn'^n by fitting the data from jp — > tt°p and 
7P — » Tr + n reactions. 

Because we did not adjust any parameter which had already been fixed in earlier fits to the irN elastic scattering, 
we found |[26l that our fits to the data were sound only up to invariant masses not exceeding W = 1.6 GeV. 
In Fig. 12 are shown our results for total cross sections (a), differential cross sections (d<r/d£l), and the photon 
asymmetry (E). The Q 2 dependence of the IV* ,7.1V vertex functions were then determined [62 1 by fitting the 
p(e, e'ir°)p and p(e, e'-K + )n data up to W = 1.6 GeV and Q 2 = 1.5 (GeV/c) 2 . 



II.E.2. Results for two-pion production reactions 

As delineated above, the dynamical coupled-channel model was constructed from fitting single-pion data. We 
then tested the efficacy of this model by examining to what extent the model could describe the ttN — > ttttN and 
the 7iV — > ttttN data. Interestingly, at the near-threshold region of W < 1.4 GeV, we found |63, 64 1 that the 
predicted total cross sections are in excellent agreement with the data. At higher W, the predicted irN —> irirN 



cross sections describe the major features of the available data reasonably well, as is shown in Fig. 13 Here, we 
further see the important role that the effects from coupled channels play. The predicted r yp —> tt^tt~p, tt tt°p 
cross sections, however, exceed the data by about a factor of two, while the fits describe, more or less, the overall 
shapes of the two-particle invariant-mass distributions. 



II.E.3. Resonance Extractions 



We define resonances as the eigenstates of the Hamiltonian with the outgoing waves being the respective decay 
channels as is described in Refs. ET1I661 . One can then show that the nucleon resonance positions are the poles 
Mr of meson-baryon scattering amplitudes as calculated from Eqs. ([T} and Q on the Riemann surface in the 
complex-E 1 plane. The coupling of meson-baryon states with the resonances can be determined by the residues 
Rn*,mb at the pole positions. Our procedures for determining Mr and Rn*,mb are further explained in our 
recent work (see: Refs. I2HI65H671 - ). 
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FIG. 13. The predicted |63| total cross sections of the ttN —¥ nnN are compared with data. The dashed curves come from 
switching off the coupled-channel effects in the DCC model of Ref. [59). 
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FIG. 14. The trajectories of the evolution of three nucleon resonances in Pn from the same bare TV* state. The results are from 
Ref. 163. 



With our method of analytic continuation into the complex plane [21 , 66 1, we are able to analyze the dynamical 
origins of the nucleon resonances within the framework of the EBAC dynamical coupled-channel model [59 1. 
This was done by examining how the resonance positions move as each of the coupled-channel couplings are 
systematically switched off. For example, as illustrated in Fig. 14 for the Pu states, this exercise revealed that two 
poles in the Roper region and the next-higher pole are associated with the same bare state on the Riemann surface. 



ILEA. Prospects and Path Forward 



During the developmental stage of the DCC analysis by the EBAC collaboration in 2006-2010, the DCC model 
parameters were determined by separately analyzing the following data sets: 7^ —> irN USTj . ^yN —> irN ll26l . 
N(e, e'n)N ll62l . irN — > ttttN ll63l . and — > irnN l64l . The very extensive data on KA and KT, production, 
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FIG. 15. Preliminary results (red bars) of the determined N* spectrum are compared with 4-star (blue bands) and 3-star (brown 
bands) states listed by the Particle Data Group. 



however, were not included in this analysis. To afford the highest-precision extraction of nucleon resonances, 
it behooves us to perform a combined and simultaneous coupled-channels analysis with all meson production 
reactions included. 

In the summer of 2010 we initiated a comprehensive eight-channel combined analysis of the world's data that 
include strange mesons in the final state, i.e. irN, —> irN, ijN, KA, KY,. The EBAC collaboration came 
to an end on March 31, 2012 and, in its place, the ANL-Osaka collaboration has taken over this DCC analysis 
task. Preliminary results of the full eight-channel combined analysis of the excited nucleon spectrum are shown 
in Fig. [15] We expect to have completed the analysis by early 2013. The ANL-Osaka analysis will then proceed 
to extract the -fN — > N* form factors from the anticipated JLab data on meson electroproduction, extending 
the momentum reach to much higher Q 2 . Further, we will explore the interpretations of the extracted resonance 
parameters in terms of the available hadron models, such as the Dyson-Schwinger-Equation model, constituent 
quark model, and Lattice QCD. Making these connections with these hadron models is needed to complete the 
DCC project with conclusive results, as is discussed in Refs.[59, 60 65l l6"8l . 



II.F. Future developments 

The CLAS collaboration has provided a wealth of data, much of which is still being analyzed. These rich 
data sets have impacted and expanded the N* program, through which reaction models can now be tuned to 
extract the j v NN* electrocouplings for CLAS 12 experiments with Q 2 > 5.0 GeV 2 , thereby enabling deeper N* 
studies [4]. Preliminary CLAS data on charged-double-pion electroproduction for photon virtualities in the range 
of 2.0 < Q 2 < 5.0 GeV 2 have recently become available [69 1 . They span the entire N* excitation region for 
W < 2.0 GeV and the statistics allow for 115 bins in W and Q 2 . The data consist of nine one-fold differential 
cross sections as is shown in Figs[6]and[7] The extension of the JM approach to higher Q 2 values up to 5.0 GeV 2 
covering the entire N* excitation region is in progress and will be completed within two years after the publication 
of this document. 

After the completion of this data analysis, electrocouplings of the Pn(1440) and Di3(1520) resonances will 
become available from both the Ntt and Tr + ir~p electroproduction channels for 0.2 < Q 2 < 5.0 GeV 2 . We 
will then have reliable information on the electrocouplings for these two states over a full range of distances that 
correspond to transitioning, wherein the quark degrees of freedom in the resonance structure dominate. The studies 
of the N* meson-baryon dressing as described in Refs. Il26ll70l strongly suggest a nearly negligible contribution 
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from the meson-baryon cloud to the Ai/% electrocouplings of the Z?i3(1520) resonance for Q 2 > 1.5 GeV 2 . 
Therefore, theoretical interpretations of already available and future CLAS results on Axii electrocouplings of the 
Di3(1520) resonance are of particular interest for approaches that are capable of describing the quark content of 
resonances based on QCD. 

Analysis of the CLAS 7r + 7r~p electroproduction data [69 1 within the framework of the JM approach will deliver 
the first information on electrocouplings of most of the high-lying excited proton states (M > 1.6 GeV) for 
2.0 < Q 2 < 5.0 GeV 2 . This information will allow us to considerably extend our knowledge on how strong 
interactions generate excited proton states having different quantum numbers. 

There will also be analyses of the available and future CLAS results on electrocouplings of all prominent N* 
states for Q 2 > 2.0 GeV 2 within the framework of the Light Cone Sum Rule approach as outlined in Chapter [v] 
that will constrain the quark-distribution amplitudes of the various N* states. Access to the quark-distribution 
amplitudes in the N* structure is of particular importance, since these amplitudes can be evaluated from QCD 
employing lattice calculations. 
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FIG. 16. Predictions of the longitudinal cross section of p° (left) and cj> (right) production versus W at Q 2 — 4 GeV 2 . For 
references to data see I71H73I and references therein. 



Information on the Q 2 evolution of non-resonant mechanisms as obtained from analyses of the CLAS data 
on single- and charged-double-pion electroproduction at Q 2 < 5.0 GeV 2 will serve as the starting point for the 
development of reaction models that will make it possible to determine the j v NN* electrocouplings from fitting 
the anticipated CLAS 12 data for Q 2 from 5.0 to 12.0 GeV 2 . 

A consistent description of a large body of observables in the Nir exclusive channels achieved within the frame- 
work of two conceptually different approaches as outlined in Section II. B. 1 and with the success of the JM model in 
describing of tt + tt~p electroproduction off protons all serve to demonstrate that the meson-baryon degrees of free- 
dom play a significant role at photon virtualities of Q 2 < 5.0 GeV 2 . Further development to the reaction models is 
needed in analyzing these exclusive channels for the anticipated CLAS 12 data, where the quark degrees of freedom 
are expected to dominate. The reaction models for the description of Tr + n , ir°p, and ir + ir~p electroproduction off 
protons for Q 2 > 5.0 GeV 2 should explicitly account for contributions from these quark degrees of freedom. At 
present, however, there is no overarching theory of hadron interactions that will offer any "off-the-shelf" approach 
at these particular distance scales, where the quark degrees of freedom dominate, but are still well inside the regime 
of the non-perturbative strong interaction. Given the state of hadronic theory, we are pursuing a phenomenological 
way for evaluating the non-resonant mechanisms for the higher-Q 2 regime. We will explore the possibilities of 
implementing different models that employ quark degrees of freedom by explicitly comparing the predictions from 
these models directly to the data. First we will start from models that employ handbag diagrams for parameterizing 
non-resonant single-pion electroproduction and we will then extend this work for a proper description of ir + n~p 
electroproduction off protons 

For the kinematics accessible at the JLab energy upgrade, one reaches the region where a description of the 
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processes of interest in terms of quark degrees of freedom applies. In this case, the calculation of cross sections 
and other observables can be performed within the handbag approach, which is based on QCD factorization of 
the scattering amplitudes in hard subprocesses, pion electroproduction off quarks, and the generalized parton 
distributions (GPDs) for p —> p or p — >• A transitions. 

In recent years, the data on the electroproduction of vector and pseudoscalar mesons have been analyzed ex- 
tensively. In particular, as described in Refs. ll7T1473l . a systematic analysis of these processes in the kinematical 
region of large Q 2 (> 3 GeV 2 ) and W larger than about 4 GeV but having small Bjorken-x (i.e. small skewness) 
has led to a set of GPDs (H, E, H, Ht, ■ ■ • ) that respect all theoretical constraints - polynomiality, positivity, 
parton distributions , and nucleon form factors. These GPDs are also in reasonable agreement with moments cal- 
culated within lattice QCD f74l and with the data on deeply virtual Compton scattering in the aforementioned 
kinematical region ll75l . On the other hand, applications in the kinematical region accessible presently from the 
6-GeV JLab data as characterized by rather large values of Bjorken-x and small W, in general, do not lead to 
agreement with experiment. Predictions for p° electroproduction, for instance, fails by an order of magnitude, 



whereas for (f> electroproduction it works quite well, as can be seen in Fig. 16 For the JLab energy upgrade, one 
can expect fair agreement between experiment and predictions for meson electroproduction evaluated from these 
set of GPDs 02). 

To describe the electroproduction of nucleon resonances in meson-baryon intermediate states in the n + n~p 
exclusive electroproduction, one needs the p — > N* transition GPDs. In principle, these GPDs are new unknown 
functions. Therefore, straightforward predictions for reactions like -f*p —> ttN* are not possible at present. In the 
large N c limit, however, one can at least relate the p —> A + GPDs to the flavor diagonal p —> p ones since the 
nucleon and the A are eigenstates of the same object, the chiral soliton 17711781 . The proton-proton GPDs always 
occur in the isovector combination F^- 3 ' = F u — F d where F is a proton-proton GPD. With the help of flavor 
symmetry one can further relate the p — > A + GPDs to all other octet-decuplet transitions. Using these theoretical 
considerations, the observables for j*p — > ttN* can be estimated. One should be aware, however, that the quality 
of the large N c and SU(3)f relations are unknown; corrections of the order of 20 to 30% are to be expected. 
One also should bear in mind that pions electroproduced by transversely-polarized virtual photons must further 
be taken into account as has been shown in Refs. Il72ll73l . Within the handbag approach, the contributions from 
such photons are related to the transversity (helicity-flip) GPDs. Despite this complication, an estimate of hard 
exclusive resonance production seems feasible. 

A well-developed program on resonance studies at high photon virtualities |7 1 will allow us to determine electro- 
couplings of several high-lying N* states with dominant Nirir decays (see the Table |ljl from the data on charged- 
double-pion electroproduction channel. However, reliable extraction of these electrocouplings for these states 
should be supported by independent analyses of other exclusive electroproduction channels having different non- 
resonant mechanisms. The r\p and KA electroproduction channels may well improve our knowledge on electro- 
couplings of the isospin 1/2 Pi3(1720) state due to isospin filtering in these exclusive channels. The studies of KY, 
and -qirN electroproduction may further offer access to the electrocouplings of the £>33(1700) and F3s(1905) reso- 
nances. More detailed studies on the feasibility of incorporating these additional exclusive channels for evaluating 
the electrocouplings of high-lying resonances are, in any case, a clear and present need. 
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III. ILLUMINATING THE MATTER OF LIGHT-QUARK HADRONS 

III A. Heart of the problem 

Quantum chromodynamics (QCD) is the strong-interaction part of the Standard Model of Particle Physics. 
Solving this theory presents a fundamental problem that is unique in the history of science. Never before have 
we been confronted by a theory whose elementary excitations are not those degrees-of-freedom readily accessible 
through experiment; i.e., whose elementary excitations are confined. Moreover, there are numerous reasons to 
believe that QCD generates forces which are so strong that less-than 2% of a nucleon's mass can be attributed 
to the so-called current-quark masses that appear in QCD's Lagrangian; viz., forces capable of generating mass 



from nothing (see Sec. III.C I. This is the phenomenon known as dynamical chiral symmetry breaking (DCSB). 
Elucidating the observable predictions that follow from QCD is basic to drawing the map that explains how the 
Universe is constructed. 

The need to determine the essential nature of light-quark confinement and dynamical chiral symmetry breaking 
(DCSB), and to understand nucleon structure and spectroscopy in terms of QCD's elementary degrees of freedom, 
are two of the basic motivations for an upgraded JLab facility. In addressing these questions one is confronted with 
the challenge of elucidating the role of quarks and gluons in hadrons and nuclei. Neither confinement nor DCSB is 
apparent in QCD's Lagrangian and yet they play the dominant role in determining the observable characteristics of 
real-world QCD. The physics of hadronic matter is ruled by emergent phenomena, such as these, which can only 
be elucidated and understood through the use of nonperturbative methods in quantum field theory. This is both the 
greatest novelty and the greatest challenge within the Standard Model. Essentially new ways and means must be 
found in order to explain precisely via mathematics the observable content of QCD. 

Bridging the gap between QCD and the observed properties of hadrons is a key problem in modern science. 
The international effort focused on the physics of excited nucleons is at the heart of this program. It addresses 
the questions: Which hadron states and resonances are produced by QCD, and how are they constituted? The N* 
program therefore stands alongside the search for hybrid and exotic mesons and baryons as an integral part of the 
search for an understanding of the strongly interacting piece of the Standard Model. 



III.B. Confinement 



Regarding confinement, little is known and much is misapprehended. It is therefore important to state clearly that 
the static potential measured in numerical simulations of quenched lattice-QCD is not related in any known way 
to the question of light-quark confinement. It is a basic feature of QCD that light-quark creation and annihilation 
effects are fundamentally nonperturbative; and hence it is impossible in principle to compute a potential between 
two light quarks IT791 1801 . Thus, in discussing the physics of light-quarks, linearly rising potentials, flux-tube 
models, etc., have no connection with nor justification via QCD. 

A perspective on confinement drawn in quantum field theory was laid out in Ref. [82| and exemplified in Sec. 2 
of Ref. |30l . It draws upon a long list of sources; e.g., Refs. [83-86 1, and, expressed simply, relates confinement to 
the analytic properties of QCD's Schwinger functions, which are often called Euclidean-space Green functions or 
propagators and vertices. For example, one reads from the reconstruction theorem [ 87 , 88 1 that the only Schwinger 
functions which can be associated with expectation values in the Hilbert space of observables; namely, the set of 
measurable expectation values, are those that satisfy the axiom of reflection positivity. This is an extremely tight 
constraint whose full implications have not yet been elucidated. 

There is a deep mathematical background to this perspective. However, for a two-point function; i.e., a prop- 
agator, it means that a detectable particle is associated with the propagator only if there exists a non-negative 
spectral density in terms of which the propagator can be expressed. No function with an inflexion point can be 
written in this way. This is readily illustrated and Fig. [17] serves that purpose. The simple pole of an observable 
particle produces a propagator that is a monotonically-decreasing convex function, whereas the evolution depicted 



in the middle-panel of Fig. 17 is manifest in the propagator as the appearance of an inflexion point at P > 0. 



To complete the illustration, consider A(k ), which is the single scalar function that describes the dressing of a 
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FIG. 17. Left panel - An observable particle is associated with a pole at timelike-P 2 , which becomes a branch point if, e.g., the 
particle is dressed by photons. Middle panel - When the dressing interaction is confining, the real-axis mass-pole splits, moving 
into pairs of complex conjugate singularities. No mass-shell can be associated with a particle whose propagator exhibits such 
singularity structure. The imaginary part of the smallest magnitude singularity is a mass-scale, \i a , whose inverse, d a — is 
a measure of the dressed-parton's fragmentation length. Right panel - A(k 2 ), the function that describes dressing of a Landau- 
gauge gluon propagator, plotted for three distinct cases. A bare gluon is described by A(fc 2 ) = 1/fc 2 (the dashed line), which 
is convex on k 2 £ (0, oo). Such a propagator has a representation in terms of a non-negative spectral density. In some theories, 
interactions generate a mass in the transverse part of the gauge-boson propagator, so that A(fc 2 ) = l/(fc 2 + m 2 ), which can 
also be represented in terms of a non-negative spectral density. In QCD, however, self-interactions generate a momentum- 
dependent mass for the gluon, which is large at infrared momenta but vanishes in the ultraviolet [81 1. This is illustrated by the 
curve labeled "IR-massive but UV-massless." With the generation of a mass-function, A(k 2 ) exhibits an inflexion point and 
hence cannot be expressed in terms of a non-negative spectral density |30|. 



Landau-gauge gluon propagator. Three possibilities are exposed in the right-panel of Fig. 17 The inflexion point 
possessed by M(p 2 ), visible in Fig. 18 entails, too, that the dressed-quark is confined. 



With the view that confinement is related to the analytic properties of QCD's Schwinger functions, the question 
of light-quark confinement may be translated into the challenge of charting the infrared behavior of QCD's univer- 
sal /3-function. (The behavior of the /3-function on the perturbative domain is well known.) This is a well-posed 
problem whose solution is a primary goal of hadron physics; e.g., Refs. ll89lj9D . It is the /3-function that is respon- 
sible for the behavior evident in Figs. [17] and [18] and thereby the scale-dependence of the structure and interactions 
of dressed-gluons and -quarks. One of the more interesting of contemporary questions is whether it is possible to 
reconstruct the /3-function, or at least constrain it tightly, given empirical information on the gluon and quark mass 
functions. 

Experiment-theory feedback within the TV* -programme shows promise for providing the latter J2J|4j[9]|. This is 
illustrated through Fig. 19 which depicts the running-gluon-mass, analogous to M(p) in Fig. 18 and the running- 
coupling determined by analyzing a range of properties of light-quark ground-state, radially-excited and exotic 
scalar-, vector- and flavored-pseudoscalar-mesons in the rainbow-ladder truncation, which is leading order in a 
symmetry -preserving DSE truncation scheme 11921 . Consonant with modern DSE- and lattice-QCD results lIBD . 
these functions derive from a gluon propagator that is a bounded, regular function of spacelike momenta, which 
achieves its maximum value on this domain at k 2 = [91. 93 . 94 1 , and a dressed-quark-gluon vertex that does not 
possess any structure which can qualitatively alter this behavior Il95ll96l . In fact, the dressed-gluon mass drawn 
here produces a gluon propagator much like the curve labeled "IR-massive but UV-massless" in the right-panel of 
Fig.[T7] 

Notably, the value of M g = m s (0) ~ 0.7 GeV is typical 11931 [941 ; and the infrared value of the coupling, 
aiu,(Mg)/ir = 2.2, is interesting because a context is readily provided. With nonperturbatively-massless gauge 
bosons, the coupling below which DCSB breaking is impossible via the gap equations in QED and QCD is 
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FIG. 18. Dressed-quark mass function, M(p): solid curves - DSE results, 1291 , "data" - lattice-QCD simulations 1281 . 
(NB. m = 70MeV is the uppermost curve. Current-quark mass decreases from top to bottom.) The constituent mass arises 
from a cloud of low-momentum gluons attaching themselves to the current-quark: DCSB is a truly nonperturbative effect that 
generates a quark mass from nothing; namely, it occurs even in the chiral limit, as evidenced by the m = curve. 




k 2 (GeV 2 ) k 2 (GeV 2 ) 

FIG. 19. Left panel - Rainbow-ladder gluon running-mass; and right panel - rainbow-ladder effective running-coupling, both 
determined in a DSE analysis of properties of light-quark mesons. The dashed curves illustrate forms for these quantities that 
provide the more realistic picture 1 89 . 103 1. (Figures drawn from Ref. 1 89).) 



a c /n « 1/3 [97-99 1 . In a symmetry -preserving regularization of a vector x vector contact-interaction used in 
rainbow-ladder truncation, a c /ir ss 0.4; and a description of hadron phenomena requires a/ir 1 [100|. With 
nonperturbatively massive gluons and quarks, whose masses and couplings run, the infrared strength required to 
describe hadron phenomena in rainbow-ladder truncation is unsurprisingly a little larger. Moreover, whilst a di- 
rect comparison between am, and a coupling, ctQLat, inferred from quenched-lattice results is not sensible, it is 
nonetheless curious that a.QL,at(0) < 0^(0) lt9TTl . It is thus noteworthy that with a more sophisticated, nonper- 
turbative DSE truncation 110111102 1. some of the infrared strength in the gap equation's kernel is shifted from the 
gluon propagator into the dressed-quark-gluon vertex. This cannot materially affect the net infrared strength re- 
quired to explain observables but does reduce the amount attributed to the effective coupling. (See, e.g., Ref. 11021 . 
wherein a(Mg) = 0.23 ir explains important features of the meson spectrum.) 
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III.C. Dynamical chiral symmetry breaking 



Whilst the nature of confinement is still debated, Fig. 18 shows that DCSB is a fact. This figure displays the 



current-quark of perturbative QCD evolving into a constituent-quark as its momentum becomes smaller. Indeed, 
QCD's dressed-quark behaves as a constituent-like-quark or a current-quark, or something in between, depending 
on the momentum with which its structure is probed. 

Dynamical chiral symmetry breaking is the most important mass generating mechanism for visible matter in the 
Universe. This may be illustrated through a consideration of the nucleon. The nucleon's a-term is a Poincare- and 
renormalization-group-invariant measure of the contribution to the nucleon's mass from the fermion mass term in 
QCD's Lagrangian lfT04l : 

a N A ==° |(m„ + m d )(N(P + K)\J(K)\N(P)) a 0.06 m N , (3) 

where J(K) is the dressed scalar vertex derived from the source \u(x)u(x) + d(x)d(x)] and mjv is the nucleon's 
mass. Some have imagined that the non-valence s-quarks produce a non-negligible contribution but it is straight- 
forward to estimate ifSOl 1 1 05H 

a s N = 0.02 -0.04 m N . (4) 

Based on the strength of DCSB for heavier quarks H 1041 , one can argue that they do not contribute a measurable 
er-term. It is thus plain that more than 90% of the nucleon's mass finds its origin in something other than the 
quarks' current-masses. 

The source is the physics which produces DCSB. As we have already mentioned, Fig. [18] shows that even in 
the chiral limit, when = = cr s N , the massless quark-parton of perturbative QCD appears as a massive 
dressed-quark to a low-momentum probe, carrying a mass-scale of approximately (l/3)mjv- A similar effect is 
experienced by the gluon-partons: they are perturbatively massless but are dressed via self-interactions, so that 



they carry an infrared mass-scale of roughly (2/3)tojv, see Fig. 19 In such circumstances, even the simplest 
symmetry-preserving Poincare-co variant computation of the nucleon's mass will produce m° N ps 3 Mq, where 
Mq sa 0.35 GeV is a mass-scale associated with the infrared behavior of the chiral-limit dressed-quark mass- 
function. The details of real-world QCD fix the strength of the running coupling at all momentum scales. That 
strength can, however, be varied in models; and this is how we know t hat if t he interaction strength is reduced, the 
nucleon mass tracks directly the reduction in Mq (see Fig. 
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and Sec. 



III.D i. Thus, the nucleon's mass is a visible 



measure of the strength of DCSB in QCD. These observations are a contemporary statement of the notions first 
expressed in Ref. 11061 . 

It is worth noting in addition that DCSB is an amplifier of explicit chiral symmetry breaking. This is why the 
result in Eq. ^ is ten-times larger than the ratio rh/m,N, where rh is the renormalization-group-invariant current- 
mass of the nucleon's valence-quarks. The result in Eq. Q is not anomalous: the nucleon contains no valence 
strangeness. Following this reasoning, one can view DCSB as being responsible for roughly 98% of the proton's 
mass, so that the Higgs mechanism is (almost) irrelevant to light-quark physics. 



The behavior illustrated in Figs. 17 - T9]has a marked influence on hadron elastic form factors. This is established, 
e.g., via comparisons between Refs. [108-1 12 1 and Refs. [100] 11131 II 141 . Owing to the greater sensitivity of 
excited states to the long-range part of the interaction in QCD 11891 11031 II 151 . we expect this influence to be 
even larger in the Q 2 -dependence of nucleon-to-resonance electrocouplings, the extraction of which, via meson 
electroproduction off protons, is an important part of the current CLAS program and studies planned with the 
CLAS12 detector j4l l5l IT HTOl . In combination with well-constrained QCD-based theory, such data can potentially, 
therefore, be used to chart the evolution of the mass function on 0.3 < p < 1.2, which is a domain that bridges the 
gap between nonperturbative and perturbative QCD. This can plausibly assist in unfolding the relationship between 
confinement and DCSB. 

In closing this subsection we re-emphasize that the appearance of running masses for gluons and quarks is a 
quantum field theoretical effect, unrealizable in quantum mechanics. It entails, moreover, that: quarks are not Dirac 
particles; and the coupling between quarks and gluons involves structures that cannot be computed in perturbation 
theory. Recent progress with the two-body problem in quantum field theory [101] has enabled these facts to be 
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FIG. 20. Evolution with current-quark mass of the ratio mrf /[3M], which varies by less-than 1% on the domain depicted. 
The calculation is described in Ref. 11071 . NB. The current-quark mass is expressed through the computed value of m 2 : 
— 0.49 GeV 2 marks the s-quark current-mass. 




i — r*—r 



J L 



N N* N- N* 



• • 

<* 



A A A- A*- Q. o* a- a*- 



FIG. 2 1 . Comparison between DSE-computed hadron masses (filled circles) and: bare baryon masses from the Excited Baryon 
Analysis Center (EBAC), |65] (filled diamonds) and JulichJ 1 17 1 (filled triangles); and experiment \ 22\,filled-squares. For the 
coupled-channels models a symbol at the lower extremity indicates that no associated state is found in the analysis, whilst a 
symbol at the upper extremity indicates that the analysis reports a dynamically-generated resonance with no corresponding bare- 
baryon state. In connection with f2-baryons the open-circles represent a shift downward in the computed results by 100 MeV. 
This is an estimate of the effect produced by pseudoscalar-meson loop corrections in A-like systems at a s-quark current-mass. 



established [ 1 16 1. One may now plausibly argue that theory is in a position to produce the first reliable symmetry- 
preserving, Poincare-invariant prediction of the light-quark hadron spectrum II 1021 . 



III.D. Mesons and Baryons: Unified Treatment 



Owing to the importance of DCSB, it is only within a symmetry-preserving, Poincare-invariant framework that 
full capitalization on the results of the N* -program is possible. One must be able to correlate the properties of 
meson and baryon ground- and excited-states within a single, symmetry-preserving framework, where symmetry- 
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FIG. 22. EBAC examined the Pi i -channel and found that the two poles associated with the Roper resonance and the next 
higher resonance were all associated with the same seed dressed-quark state. Coupling to the continuum of meson-baryon final 
states induces multiple observed resonances from the same bare state. In EBAC's analysis, all PDG-identified resonances were 
found to consist of a core state plus meson-baryon components. (Adapted from Ref. 1 65 1.) 



preserving includes the consequence that all relevant Ward-Takahashi identities are satisfied. This is not to say that 
constituent-quark-like models are worthless. As will be seen in this article, they are of continuing value because 
there is nothing better that is yet providing a bigger picture. Nevertheless, such models have no connection with 
quantum field theory and therefore not with QCD; and they are not "symmetry-preserving" and hence cannot 
veraciously connect meson and baryon properties. 

An alternative is being pursued within quantum field theory via the Faddeev equation. This analogue of the 
Bethe-Salpeter equation sums all possible interactions that can occur between three dressed-quarks. A tractable 
equation [ 1 18 1 is founded on the observation that an interaction which describes color-singlet mesons also gener- 
ates nonpointlike quark-quark (diquark) correlations in the color-anti triplet channel 01 191 . The dominant correla- 
tions for ground state octet and decuplet baryons are scalar (0 + ) and axial-vector (1 + ) diquarks because, e.g., the 
associated mass-scales are smaller than the baryons' masses and their parity matches that of these baryons. On 
the other hand, pseudoscalar (0~) and vector (1~) diquarks dominate in the parity-partners of those ground states 
OTl 1 1071 . This approach treats mesons and baryons on the same footing and, in particular, enables the impact of 
DCSB to be expressed in the prediction of baryon properties. 

Incorporating lessons learnt from meson studies [ 120|, a unified spectrum of u, d, s-quark hadrons was obtained 
using symmetry-preserving regularization of a vector x vector contact interaction 1311110711 . These studies simul- 
taneously correlate the masses of meson and baryon ground- and excited-states within a single framework; and 
in comparison with relevant quantities, they produce rms < 15%, where rms is the root-mean-square-relative- 



error/degree-of freedom. As indicated by Fig. 21 they uniformly overestimate the PDG values of meson and 
baryon masses [22 1. Given that the truncation employed omits deliberately the effects of a meson-cloud in the 
Faddeev kernel, this is a good outcome, since inclusion of such contributions acts to reduce the computed masses. 



Following this line of reasoning, a striking result is agreement between the DSE-computed baryon masses 11107 



and the bare masses employed in modern coupled-channels models of pion-nucleon reactions [65 117), see Fig. 21 
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TABLE II. Bare masses (GeV) determined in an Argonne-Osaka coupled-channels analysis of single- and double-pion electro- 
production reactions compared with DSE results for the mass of each baryon's dressed-quark core. The notation is as follows: 
Pn corresponds to the iV(1440); Sn to the iV(1535) and the second state in this partial wave; P33 to the A and the next state 
in this partial wave; and D33 to the parity partner of the A. The rms-|rel. error] = 9.4 ± 5.7%. 
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The Roper resonance is very interesting. The DSE studies [31, 107 1 produce a radial excitation of the nucleon at 
1.82 ± 0.07 GeV. This state is predominantly a radial excitation of the quark-diquark system, with the diquark 
correlations in their ground state. Its predicted mass lies precisely at the value determined in the analysis of 
Ref. 11651 . This is significant because for almost 50 years the Roper resonance has defied understanding. 

Discovered in 1963/64 [ 121], the Roper appears to be an exact copy of the proton except that its mass is 50% 
greater. The mass was the problem: hitherto it could not be explained by any symmetry-preserving QCD-based 
tool. That has now changed. Combined, see Fig. 22 Refs. l3"Tll65l[l071 demonstrate that the Roper resonance is 
indeed the proton's first radial excitation, and that its mass is far lighter than normal for such an excitation because 
the Roper obscures its dressed-quark-core with a dense cloud of pions and other mesons. Such feedback between 
QCD-based theory and reaction models is critical now and for the foreseeable future, especially since analyses of 
CLAS data on nucleon-resonance electrocouplings suggest strongly that this structure is typical; i.e., most low- 
lying N* -states can best be understood as an internal quark-core dressed additionally by a meson cloud [18). This 
is highlighted further by a comparison between the DSE results and the bare masses obtained in the most complete 
Argonne-Osaka coupled-channels analysis to date, see Table [IT 112211 

Additional analysis within the framework of Refs. 113111107 1 suggests a fascinating new possibility for the Roper, 



which is evident in Table. Ill The nucleon ground state is dominated by the scalar diquark, with a significantly 
smaller but nevertheless important axial-vector diquark component. This feature persists in solutions obtained 
with more sophisticated Faddeev equation kernels (see, e.g., Table 2 in Ref. Ill 101 ). From the perspective of 
the nucleon's parity partner and its radial excitation, the scalar diquark component of the ground-state nucleon 
actually appears to be unnaturally large. Expanding the study to include baryons containing one or more s-quarks, 
the picture is confirmed: the ground state N, A, S, 5 are all characterized by ~ 80% scalar diquark content OTI . 
whereas their parity partners have a 50 — 50 mix of J = 0, 1 diquarks. 

One can nevertheless understand the structure of the octet ground-states. As with so much else in hadron physics, 
the composition of these flavor octet states is intimately connected with DCSB. In a two-color version of QCD, 
scalar diquarks are Goldstone modes II 1231 11241 . (This is a long-known result of Pauli-Gursey symmetry.) A 
"memory" of this persists in the three-color theory, for example: in low masses of scalar diquark correlations; and 
in large values of their canonically normalized Bethe-Salpeter amplitudes and hence strong quark+quark— diquark 
coupling within the octet ground-states. (A qualitatively identical effect explains the large value of the ttN cou- 
pling constant and its analogues involving other pseudoscalar-mesons and octet-baryons.) There is no commensu- 
rate enhancement mechanism associated with the axial-vector diquark correlations. Therefore the scalar diquark 
correlations dominate within octet ground-states. 

Within the Faddeev equation treatment, the effect on the first radial excitations is dramatic: orthogonality of the 
ground- and excited-states forces the radial excitations to be constituted almost entirely from axial-vector diquark 
correlations. It is critical to check whether this outcome survives with Faddeev equation kernels built from a 
momentum-dependent interaction. 

This brings us to another, very significant observation; namely, the match between the DSE-computed level 
ordering and that of experiment, something which has historically been difficult for models to obtain (see, e.g., the 
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TABLE III. Diquark content of the baryons' dressed-quark cores, computed with a symmetry-preserving regularization of a 
vector x vector contact interaction 11251 . 
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discussion in Ref. 1 126 1) and is not achieved in contemporary numerical simulations of lattice-regularized QCD 
(see, e.g., Ref. 01271 ). In particular, the DSE calculations produce a parity-partner for each ground-state that is 
always more massive than its first radial excitation so that, in the nucleon channel, e.g., the first J p = ~ state 

lies above the second J p = | + state. 

A veracious expression of DCSB in the meson spectrum is critical to this success. One might ask why and 
how? It is DCSB that both ensures the dressed-quark-cores of pseudoscalar and vector mesons are far lighter 
than those of their parity partners and produces strong quark+antiquark— meson couplings, which are expressed 
in large values for the canonically normalized Bethe-Salpeter amplitudes (Table 3 in Ref. fTTlD . The remnants 
of Pauli-Giirsey symmetry described previously entail that these features are carried into the diquark sector: as 
evident in Fig. 3 and Table 5 of Ref. PTTl and their comparison with Fig. 2 and Table 3 therein. The inflated 
masses but, more importantly, the suppressed values of the Bethe-Salpeter amplitudes for negative-parity diquarks, 
in comparison with those of positive-parity diquarks, guarantee the computed level ordering: attraction in a given 
channel diminishes with the square of the Bethe-Salpeter amplitude (see App. C in Ref. [3TJ). Hence, an approach 
within which DCSB cannot be realized or a simulation whose parameters are such that the importance of DCSB is 
suppressed will both necessarily have difficulty reproducing the experimental ordering of levels. 

The computation of spectra is an important and necessary prerequisite to the calculation of nucleon transition 
form factors, the importance of which is difficult to overestimate given the potential of such form factors to assist 
in charting the long-range behavior of QCD's running coupling. To place this in context, Refs. 11001 11071 II 131 
114 1 explored the sensitivity of a range of hadron properties to the running of the dressed-quark mass-function. 



These studies established conclusively that static properties are not a sensitive probe of the behavior in Figs. 18 



19 viz., regularized via a symmetry-preserving procedure, a vector x vector contact-interaction predicts masses, 
magnetic and quadrupole moments, and radii that are practically indistinguishable from results obtained with the 
most sophisticated QCD-based interactions available currently [89, 128|. 



III.E. Nucleon to Resonance Transition Form Factors 



The story is completely different, however, with the momentum-dependence of form factors; e.g., in the case of 
the pion , the difference between the form factor obtained with M(p) = constant and that derived from M{p 2 ) in 
Fig. 18 is dramatically apparent for Q 2 > M 2 (p — 0) II 131 . The study of diquark form factors in Ref. 11001 has 



enabled another reference computation to be undertaken; namely, nucleon elastic and nucleon-to-Roper transition 
form factors Il68l . It shows that axial-vector-diquark dominance of the Roper, Table |HT] has a material impact on 
the nucleon-to-Roper transition form factor. 

We choose to illustrate the analysis of Ref. Il68l via Fig. [23] The figure displays results obtained using a 
light-front constituent-quark model Il24l . which employed a constituent-quark mass of 0.22 GeV and identical 
momentum-space harmonic oscillator wave functions for both the nucleon and Roper (width = 0.38 GeV) but with 
a zero introduced for the Roper, whose location was fixed by an orthogonality condition. The quark mass is smaller 
than that which is typical of DCSB in QCD but a more significant difference is the choice of spin-flavour wave 
functions for the nucleon and Roper. In Ref. [24] they are simple SU (6) x 0(3) S- wave states in the three-quark 
center-of-mass system, in contrast to the markedly different spin-flavour structure produced by Faddeev equation 
analyses of these states. 
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FIG. 23. Helicity amplitudes for the 7*p — > Pn(1440) transition, with x — Q 2 /m 2 N : A 1 / 2 (upper panel); and S1/2 (lower 
panel). Solid curves - DSE computation of Ref. |68 1, obtained using a contact interaction but amended here via an estimate of 
the impact of the dressed-quark mass in Fig.[T8] which softens the x > 1 behavior without affecting x < 1; dashed curves - 
the light-front constituent quark model results from Ref. [24 1 ; long-dash-dot curves - the light-front constituent quark model 
results from Ref. 1 129 1; short-dashed curves - a smooth fit to the bare form factors inferred in Ref. 12 1 1 1621 \l 301 : and data - 
Refs. fT71[T8ll23l. 

Owing to this, in Fig. [23] we also display the light-front quark model results from Ref. Ml 291 . It is stated 
therein that large effects accrue from "configuration mixing;" i.e., the inclusion of 5[/(6)-breaking terms and 
high-momentum components in the wave functions of the nucleon and Roper. In particular, that configuration 
mixing yields a marked suppression of the calculated helicity amplitudes in comparison with both relativistic and 
non-relativistic results based on a simple harmonic oscillator Ansatz for the baryon wave functions, as used in 
Ref. (31. 

There is also another difference; namely, Ref. |129| employs Dirac and Pauli form factors to describe the in- 
teraction between a photon and a constituent-quark 11311 . As apparent in Fig. 2 of Ref. B129I , they also have a 
noticeable impact, providing roughly half the suppression on 0.5 < Q 2 /GeV 2 < 1.5. The same figure also 
highlights the impact on the form factors of high-momentum tails in the nucleon and Roper wave functions. 

In reflecting upon constituent-quark form factors, we note that the interaction between a photon and a dressed- 
quark in QCD is not simply that of a Dirac fermion 01 161I132T4137I . Moreover, the interaction of a dressed-quark 
with the photon in Ref. l68l is also modulated by form factors, see Apps. A3, C6 therein. On the other hand, the 
purely phenomenological form factors in Refs. 11129111311 are inconsistent with a number of constraints that apply 
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to the dressed-quark-photon vertex in quantum field theory; e.g., the dressed-quark's Dirac form factor should 
approach unity with increasing Q 2 and neither its Dirac nor Pauli form factors may possess a zero. Notwithstanding 
these observations, the results from Ref. Ml 291 are more similar to the DSE curves than those in Ref. Il24ll . 

In an interesting new development, the study of Ref. [24] has been updated [138]. The new version models the 
impact of a running dressed-quark mass within the light-front formulation of quantum mechanics and yields results 
that are also closer to those produced by the DSE analysis. 

Helicity amplitudes can also be computed using the Argonne-Osaka Collaboration's dynamical coupled- 
channels framework |59|. In this approach, one imagines that a Hamiltonian is defined in terms of bare baryon 
states and bare meson-baryon couplings; the physical amplitudes are computed by solving coupled-channels equa- 
tions derived therefrom; and the parameters characterizing the bare states are determined by requiring a good fit 
to data. In connection with the -f*p —> P 11 (1440) transition, results are available for both helicity amplitudes 



ET1 l62l [1301 . The associated bare form factors are reproduced in Fig. 23 for Q 2 < 1.5GeV 2 we depict a smooth 



interpolation; and for larger Q an extrapolation based on perturbative QCD power laws (Ai ~ l/Q ~ £1) 



The bare form factors are evidently similar to the results obtained in Ref. 11681 and in Ref. [129]: both in mag- 
nitude and Q 2 -evolution. Regarding the transverse amplitude, Ref. ||21~I argues that the bare component plays an 
important role in changing the sign of the real part of the complete amplitude in the vicinity of Q 2 = 0. In this case 
the similarity between the bare form factor and the DSE results is perhaps most remarkable - e.g., the appearance 
of the zero in A 1 , and the Q 2 = magnitude of the amplitude (in units of 10~ 3 GcV^ 1 ^ 2 ) 

Ref.(2! Ref. (HO Ref. EH [62j [J30l Ref. (61 

A i(0) -35.1 -32.3 -18.6 -16.3 ' (5) 

These similarities strengthen support for an interpretation of the bare-masses, -couplings, etc., inferred via coupled- 
channels analyses, as those quantities comparable with hadron structure calculations that exclude the meson-baryon 
coupled-channel effects which are determined by multichannel unitarity conditions. 

An additional remark is valuable here. The Argonne-Osaka Collaboration computes electroproduction form 
factors at the resonance pole in the complex plane and hence they are complex-valued functions. Whilst this is 
consistent with the standard theory of scattering 111391 , it differs markedly from phenomenological approaches that 
use a Breit-Wigner parametrization of resonant amplitudes in fitting data. As concerns the j*p — > Pn(1440) 
transition, the real parts of the Argonne-Osaka Collaboration's complete amplitudes are qualitatively similar to the 
results in Refs. [IV, 18 23 1 but the Argonne-Osaka Collaboration's amplitudes also have sizeable imaginary parts. 
This complicates a direct comparison between theory and extant data. 



III.F. Prospects 

A compelling goal of the international theory effort that works in concert with the A* -program is to understand 
how the interactions between dressed-quarks and -gluons create nucleon ground- and excited-states, and how 
these interactions emerge from QCD. This compilation shows no single approach is yet able to provide a unified 
description of all N* phenomena; and that intelligent reaction theory will long be necessary as a bridge between 
experiment and QCD-based theory. Nonetheless, material progress has been made since the release of the White 
Paper on "Theory Support for the Excited Baryon Program at the JLab 12-GeV Upgrade" [5|, in developing 
strategies, methods and approaches to the physics of nucleon resonances. Some of that achieved via the Dyson- 
Schwinger equations is indicated above. Additional contributions relevant to the N* program are: verification of 
the accuracy of the diquark truncation of the quark-quark scattering matrix within the Faddeev equation 111 121 : and 
a computation of the A — > irN transition form factor 1 140 1. 

A continued international effort is necessary if the goal of turning experiment into a probe of the dressed- 
quark mass function and related quantities is to be achieved. In our view, precision data on nucleon-resonance 
transition form factors provides a realistic means by which to constrain empirically the momentum evolution of 
the dressed-quark mass function and therefrom the infrared behavior of QCD's /3-function; in particular, to locate 
unambiguously the transition boundary between the constituent- and current-quark domains that is signalled by the 



sharp drop apparent in Fig. 18 That drop can be related to an inflexion point in QCD's /3-function. Contemporary 
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theory indicates that this transition boundary lies at p 2 ~ 0.6 GeV 2 . Since a probe's input momentum Q is 
principally shared equally amongst the dressed-quarks in a transition process, then each can be considered as 
absorbing a momentum fraction Q/3. Thus in order to cover the domain p 2 € [0.5, 1.0] GeV 2 one requires 
Q 2 € [5, 10] GeV 2 ; i.e., the upgraded JLab facility. 

In concrete terms, a DSE study of the N —> 7V(1535) transition is underway, using the contact-interaction, for 



comparison with data 11141 1171 and other computations 1 141 1; and an analysis of the TV — > A transition has begun, 
with the aim of revealing the origin of the unexpectedly rapid Q 2 -evolution of the magnetic form factor in this 
process. 

At the same time, the Faddeev equation framework of Ref. 11 101 . is being applied to the N — > 7V(1440) transi- 
tion. The strong momentum dependence of the dressed-quark mass function is an integral part of this framework. 
Therefore, in this study it will be possible, e.g., to vary artificially the position of the marked drop in the dressed- 
quark mass function and thereby identify experimental signatures for its presence and location. In addition, it will 
provide a crucial check on the results in Table III It is notable that DCSB produces an anomalous electromagnetic 
moment for the dressed-quark. This is known to produce a significant modification of the proton's Pauli form fac- 
tor at Q 2 < 2 GeV 2 II 1421 . It is also likely to be important for a reliable description of F 2 * in the nucleon-to-Roper 
transition. 

The Faddeev equation framework of Ref. [ 1 10 1 involves parametrizations of the dressed-quark propagators that 
are not directly determined via the gap equation. An important complement would be to employ the ab initio 
rainbow -ladder truncation approach of Ref. 01111 11121 in the computation of properties of excited-state baryons, 
especially the Roper resonance. Even a result for the Roper's mass and its Faddeev amplitude would be useful, 
given the results in Table III In order to achieve this, however, technical difficulties must be faced and overcome. 
Here there is incipient progress, made possible through the use of generalized spectral representations of propaga- 
tors and vertices. 

In parallel with the program outlined herein, an effort is beginning with the aim of providing the reaction theory 
necessary to make reliable contact between experiment and predictions based on the dressed-quark core. While 
rudimentary estimates can and will be made of the contribution from pseudoscalar meson loops to the dressed- 
quark core of the nucleon and its excited states, a detailed comparison with experiment will only follow when the 
DSE-based results are used to constrain the input for dynamical coupled channels calculations. 
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IV. iV* PHYSICS FROM LATTICE QCD 
IV.A. Introduction 

Quantum ChromoDynamics (QCD), when combined with the electroweak interactions, underlies all of nuclear 
physics, from the spectrum and structure of hadrons to the most complex nuclear reactions. The underlying sym- 
metries that are the basis of QCD were established long ago. Under very modest assumptions, these symmetries 
predict a rich and exotic spectrum of QCD bound states, few of which have been observed experimentally. While 
QCD predicts that quarks and gluons are the basic building blocks of nuclear matter, the rich structure that is exhib- 
ited by matter suggests there are underlying collective degrees of freedom. Experiments at nuclear and high-energy 
physics laboratories around the world measure the properties of matter with the aim to determine its underlying 
structure. Several such new experiments worldwide are under construction, such as the 12-GeV upgrade at Jeffer- 
son Lab's electron accelerator, its existing the experimental halls, as well as the new Hall D. 

To provide a theoretical determination and interpretation of the spectrum, ah initio computations within lattice 
QCD have been used. Historically, the calculation of the masses of the lowest-lying states, for both baryons and 
mesons, has been a benchmark calculation of this discretized, finite-volume computational approach, where the 
aim is well-understood control over the various systematic errors that enter into a calculation; for a recent review, 
see II 1431 . However, there is now increasing effort aimed at calculating the excited states of the theory, with several 
groups presenting investigations of the low-lying excited baryon spectrum, using a variety of discretizations, num- 
bers of quark flavors, interpolating operators, and fitting methodologies (Refs. IU44H1471 ). Some aspects of these 
calculations remain unresolved and are the subject of intense effort, notably the ordering of the Roper resonance 
in the low-lying nucleon spectrum. 

The Hadron Spectrum Collaboration, involving the Lattice Group at Jefferson Lab, Carnegie Mellon University, 
University of Maryland, University of Washington, and Trinity College (Dublin), is now several years into its pro- 
gram to compute the high-lying excited- state spectrum of QCD, as well as their (excited-state) electromagnetic 
transition form factors up to Q 2 ~ 10 GeV 2 . This program has been utilizing "anisotropic" lattices, with finer tem- 
poral than spatial resolution, enabling the hadron correlation functions to be observed at short temporal distances 
and hence many energy levels to be extracted 1148111491 . Recent advances suggest that there is a rich spectrum of 
mesons and baryons, beyond what is seen experimentally. In fact, the HSC's calculation of excited spectra, as well 
as recent successes with GPUs, were featured in Selected FY 10 Accomplishments in Nuclear Theory in the FY12 
Congressional Budget Request. 



IV.B. Spectrum 

The development of new operator constructions that follow from continuum symmetry constructions has al- 
lowed, for the first time, the reliable identification of the spin and masses of the single-particle spectrum at a 
statistical precision at or below about 1%. In particular, the excited spectrum of isovector as well as isoscalar 
mesons (Refs. B15H41531 ) shows a pattern of states, some of which are familiar from the qq constituent quark 
model, with up to total spin J — A and arranged into corresponding multiplets. In addition, there are indications 
of a rich spectrum of exotic J PC states, as well as a pattern of states interpretable as non-exotic hybrids 11541 . 
The pattern of these multiplets of states, as well as their relative separation in energy, suggest a phenomenology of 
constituent quarks coupled with effective gluonic degrees of freedom. In particular, the pattern of these exotic and 
non-exotic hybrid states appears to be consistent with a bag-model description and inconsistent with a flux-tube 
model fl54l . 

Recently, this lattice program has been extended into the baryon spectrum, revealing for the first time, the 
excited-state single-particle spectrum of nucleons and Deltas along with their total spin up to J = | in both positive 
and negative parity 1 1 27 1 . The results for the lightest-mass ensemble are shown in Fig. [24] There was found a high 
multiplicity of levels spanning across J p which is consistent with SU(6) ® 0(3) multiplet counting, and hence 
with that of the non-relativistic qqq constituent quark model. In particular, the counting of levels in the low-lying 
negative -parity sectors are consistent with the non-relativistic quark model and with the observed experimental 
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FIG. 24. Results from Ref. 1 127] showing the spin-identified spectrum of Nucleons and Deltas from the lattices at = 
396 MeV, in units of the calculated £1 mass. The states identified through their spectral overlaps are shown, and the full 
number of states expected from SU (6) ® 0(3) counting are found. 



states 1 22 1 . The spectrum observed in the first-excited positive-parity sector is also consistent in counting with 
the quark model, but the comparison with experiment is less clear, with the quark model predicting more states 
than are observed experimentally, spurring phenomenological investigations to explain the discrepancies (e.g., see 
Refs. 11221 [T26l[l55til59"l ). 

In addition, it was found that there is significant mixing among each of the allowed multiplets, including the 
20-plet that is present in the non-relativistic qqq quark model but does not appear in quark-diquark models II 1 5711 
(see in particular Ref. B160I ). These results lend credence to the assertion that there is no "freezing" of degrees 
of freedom with respect to those of the non-relativistic quark model. These qualitative features of the calculated 
spectrum extend across all three of the quark-mass ensembles studied. Furthermore, no evidence was found for the 
emergence of parity-doubling in the spectrum 1 161 1. 

The results for the baryon spectrum investigations from Ref. [ 127 1 suggest that to faithfully describe the excited 
spectrum requires the use of non-local operator constructions. Fig. 25 shows a comparison of the results for the 
i + 



Nucleon J = I spectrum taken from Ref. [ 127 1 and with some other calculations in full QCD from Refs. II 146 



2 

11501 . Up to some lattice scale ambiguity, it is clear there are a distinctly different number of states found at 
comparable pion masses. Namely, there sis four nearly degenerate excited states found at approximately 2.2 GeV, 
and three nearly degenerate states near 2.8 GeV. The new results suggest the observed excited J = | + states are 
admixtures of radial excitations as well as Z?-wave and anti-symmetric P-wave structures, and the inclusion of 
operators featuring such structures is essential to resolve the degeneracy of states. 

It was argued that the extracted N and A spectrum can be interpreted in terms of single-hadron states, and based 
on investigations in the meson sector II 1521 and initial investigations of the baryon sector at a larger volume 11271 . 
little evidence was found for multi-hadron states. To study multi-particle states, and hence the resonant nature 
of excited states, operator constructions with a larger number of fermion fields are needed. Such constructions 
are in progress |162|, and it is believed that the addition of these operators will lead to a denser spectrum of 
states. With suitable understanding of the discrete energy spectrum of the system, the Liischer formalism [163| 
and its inelastic extensions (for example, see Ref. |164|) can be used to extract the energy dependent phase shift 
for a resonant system, such as has been performed for the I = 1 p system [165|. The energy of the resonant 
state is determined from the energy dependence of the phase shift. It is this resonant energy that is suitable for 
chiral extrapolations. Suitably large lattice volumes and smaller pion masses are needed to adequately control the 
systematic uncertainties in these calculations. 
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FIG. 25. Comparison of results for the nucleon J = | channel. The results shown in grey are from Ref. 1 150 1, while those 
in orange are from Ref. [ 146 1. Note that data are plotted using the scale-setting scheme in the respective papers. Results from 
Ref. 11271 are shown in red (the ground state), green and blue. At the lightest pion mass, there is a clustering of four states as 
indicated near 2 GeV, while there are three nearly degenerate states 2.7 GeV. Operators featuring the derivative constructions 
discussed in Ref. 1 127 1 feature prominently in these excited states, suggesting previous results are insensitive to these excited 
states because the operator bases used were incomplete. 



IV.C. Electromagnetic transition form factors 



The measurement of the excited-to-ground state radiative transition form factors in the baryon sector provides a 
probe into the internal structure of hadrons. Analytically, these transition form factors can be expressed in terms 
of matrix elements between states (N(pf)\V f2 (q)\N* (pi)) where is a vector (or possibly axial-vector) current 
with some four-momentum q = pj — pi between the final (pf) and initial (pi) states. This matrix element can be 
related to the usual form factors F^(q 2 ) and F 2 *(g 2 ). However, the exact meaning as to the initial state |7V*) is 
the source of some ambiguity since in general it is a resonance. In particular, how is the electromagnetic decay 
disentangled from that of some Nir hadronic contribution? 

Finite-volume lattice-QCD calculations are formulated in Euclidean space, and as such, one does not directly 
observe the imaginary part of the pole of a resonant state. However, the information is encoded in the volume and 
energy dependence of excited levels in the spectrum. Liischer's formalism 116311 and its many generalizations show 
how to relate the infinite-volume energy-dependent phase shifts in resonant scattering to the energy dependence 
of levels determined in a continuous but finite-volume box in Euclidean space. In addition, infinite-volume matrix 
elements can be related to those in finite-volume [ 166] up to a factor which can be determined from the derivative 
of the phase shift. 
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FIG. 26. Exploratory evaluation of the F 1 11 (Q 2 ) and -F^" 11 (Q 2 ) form factors for the transition from the ground state proton 
to the excited Pn(1440) state carried out within the framework of unquenched LQCD [32| on the Nf = 2+1 anisotropic 
lattices in comparison with the experimental data from CLAS (black bullets) 1 17|. LQCD results shown by green diamonds, 
red squares, and golden triangles are obtained for pion masses of 390 MeV, 450 MeV, and 875 MeV, whose volumes are 3, 2.5, 
2.5 fm, respectively. 



For the determination of transition form factors, what all this means in practice is that one must determine the 
excited-state transition matrix element from each excited level in the resonant region of a state, down to the ground 
state. The excited levels and the ground state might each have some non-zero momentum, arising in some Q 2 
dependence. In finite volume, the transition form factors are both Q 2 and energy dependent, the latter coming 
from the discrete energies of the states within the resonant region. The infinite-volume form factors are related 
to these finite-volume form factors via the derivative of the phase shift as well as another kinematic function. 
Sitting close to the resonant energy, in the large volume limit the form factors become independent of the energy 
as expected. 

The determination of transition form factors for highly excited states was first done in the charmonium sector 
with quenched QCD H1671I1681 . Crucial to these calculations was the use of a large basis of non-local operators 
to form the optimal projection onto each excited level. In a quenched theory, the excited charmonium states are 
stable and have no hadronic decays, thus there is no correction factor. 

The determination of the electromagnetic transitions in light-quark baryons will eventually require the determi- 
nation of the transition matrix elements from multiple excited levels in the resonance regime, the latter determined 
through the spectrum calculations in the previous section. However, as a first step, the Q 2 dependence of tran- 
sition form factors between the ground and first-excited state can be investigated within a limited basis. These 
first calculations of the Ff^(Q 2 ) excited transition levels, in Refs. Il32lll691 already have shown many interesting 
features. 

The first calculations of the P\x — > "fN transition form factors were performed a few years ago using the 
quenched approximation [169]. Since then, these calculations have been extended to full QCD with two light 
quarks and one strange quark (Nf = 2 + 1) using the same anisotropic lattice ensembles as for the spectrum 
calculations. Preliminary results [32J_of the Q 2 dependence of the first-excited nucleon (the Roper) to the ground- 



state proton, Ff ^, are shown in Fig. 26 These results focus on the low-Q 2 region. At the unphysical pion masses 



used, some points are in the time-like region. What is significant in these calculations with full-QCD lattice 
ensembles is that the sign of F2 at low Q 2 has flipped compared to the quenched result, which had relatively 
mild Q 2 dependence at similar pion masses. These results suggest that at low Q 2 the pion-cloud dynamics are 
significant in full QCD. 

The results so far are very encouraging, and the prospects are quite good for extending these calculations. The 
use of the larger operator basis employed in the spectrum calculations, supplemented with multi -particle operators, 
and including the correction factors from the resonant structure contained in phase shifts, should allow for the 
determination of multiple excited-level transition form factors up to about Q 2 3 GeV 2 . 
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FIG. 27. Pion form factor utilizing an extended basis of smearing functions to increase the range of Q 2 with multiple 
pion masses at 580, 875, 1350 MeV. The experimental points are shown as (black) circles while the lowest gray band is the 
extrapolation to the physical pion mass using lattice points from Ref. 1 1701 . 



IV.D. Form factors at Q 2 w 6 GeV 2 



The traditional steps in a lattice form-factor calculation involve choosing suitable creation and annihilation 
operators with the quantum numbers of interest, and typically where the quark fields are spatially smeared so as to 
optimize overlap with the state of interest, often the ground state. These smearing parameters are typically chosen 
to optimize the overlap of a hadron at rest or at low momentum. As the momentum is increased, the overlap of 
the boosted operator with the desired state in flight becomes small and statistically noisy. One method to achieve 
high Q 2 is to decrease the quark smearing, which has the effect of increasing overlap onto many excited states. 
By choosing a suitably large basis of smearing, one can then project onto the desired excited state at high(er) 
momentum. This technique can extend the range of Q 2 in form-factor calculations until lattice discretization 
effects become dominant. An earlier version of this technique (with smaller basis) was used for a quenched 



calculation of the Roper transition form factor reaching about 6 GeV 2 11169111711 . Figure 27 shows an example 
from Nf = 2 + 1 at 580, 875, 1350 MeV pion masses using extended basis to extract pion form factors with Q 2 
reaching nearly 7 GeV 2 117011 for the highest-mass ensemble. The extrapolated form factor at the physical pion 
mass shows reasonable agreement with JLab precision measurements. Future attempts will focus on decreasing 
the pion masses and exploring Q 2 -dependence of pion form factors for yet higher Q 2 . 

As before, these form-factor calculations need to be extended to use a larger operator basis of single and multi- 
particle operators to overlap with the levels within the resonant region of the excited state, say the Roper. These 
operator constructions are suitable for projecting onto excited states with high momentum, as demonstrated in 
Ref. |162|. Future work will apply these techniques to form-factor calculations. 



IV.E. Form factors at high Q 2 > 10 GeV 2 

At very high Q 2 , lattice discretization effects can become quite large. A costly method to control these effects 
is to go to much smaller lattice spacing, basically a ~ An alternative method that was been devised long 

ago is to use renormalization-group techniques 11721 . and in particular, step-scaling techniques introduced by the 
ALPHA collaboration. The step-scaling method was initially applied to compute the QCD running coupling and 
quark masses. The technique was later extended to handle heavy-quark masses with a relativistic action 1173111741 . 
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The physical insight is that the heavy-quark mass dependence of ratios of observables is expected to be milder than 
the observable itself. For form factors, the role of the large heavy-quark mass scale is now played by the large 
momentum scale Q. Basically, the idea is to construct ratios of observables (form factors) such that the overall Q 2 
dependence is mild, and that suitable products of these ratios, evaluated at different lattice sizes and spacings, can 
be extrapolated to equivalent results at large volume and fine lattice spacing. The desired form factor is extracted 
from the ratios. 

The technique, only briefly sketched here, is being used now in a USQCD lattice-QCD proposal by D. Renner 
(Ref. 117511 ) to compute the pion form-factor at large Q 2 , and the technique is briefly discussed in Ref. II 1 701 - In 
principle, the same technique can be used to compute excited-state transition form factors, and although feasibility 
has yet to be established, it seems worth further investigation. 

IV.F. Outlook 

There has been considerable recent progress in the determination of the highly excited spectrum of QCD using 
lattice techniques. While at unphysically large pion masses and small lattice volumes, already some qualitative 
pictures of the spectrum of mesons and baryons is obtained. With the inclusion of multi-hadron operators, the 
outlook is quite promising for the determination of the excited spectrum of QCD. Anisotropic lattice configurations 
with several volumes are available now for pion masses down 230 MeV. Thus, it seems quite feasible to discern 
the resonant structure for at least a few low-lying states of mesons and baryons, of course within some systematic 
uncertainties, in the two-year timeframe. One of the more open questions is how to properly handle multi-channel 
decays which becomes more prevalent for higher-lying states. Some theoretical work has already been done using 
coupled-channel methods, but more work is needed and welcomed. 

With the spectrum in hand, it is fairly straightforward to determine electromagnetic transition form factors for 
the lowest few levels of N*, and up to some moderate Q 2 of a few GeV 2 , in the two-year time-frame. Baryon 
form factors will probably continue to drop purely disconnected terms from the current insertion. Meson transition 
form factors, namely an exotic to non-exotic meson will be the first target in the short time-frame (less than two 
years), with the aim to determine photo-couplings. It might well be possible that with the new baryon operator 
techniques developed, the transition form factors can be extracted to Q 2 ss 6 GeV 2 . Going to an isotropic lattice 
with a small lattice spacing, it seems feasible to reach higher Q 2 , say 10 GeV 2 , and this could be available in less 
than five years. To reach Q 2 ^> 10 GeV 2 will probably require step-scaling techniques. The high-Q 2 limit is of 
considerable interest since it allows for direct comparisons with perturbative methods. 
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V. LIGHT-CONE SUM RULES: A BRIDGE BETWEEN ELECTROCOUPLINGS AND DISTRIBUTION 

AMPLITIDES OF NUCLEON RESONANCES 

We expect that at photon virtualities from 5 to 10 GeV 2 of CLAS12 the electroproduction cross sections of nu- 
clear resonances will become amenable to the QCD description in terms of quark partons, whereas the description 
in terms of meson-baryon degrees of freedom becomes much less suitable than at smaller momentum transfers. 
The major challenge for theory is that quantitative description of form factors in this transition region must include 
nonperturbative contributions. In Ref. 114111 we have suggested to use a combination of light-cone sum rules (LC- 
SRs) and lattice calculations. To our opinion this approach presents a reasonable compromise between theoretical 
rigour and the necessity to make phenomenologically relevant predictions. 



V.A. Light-cone wave functions and distribution amplitudes 

The quantum-mechanical picture of a nucleon as a superposition of states with different number of partons 
assumes the infinite momentum frame or light-cone quantization. Although a priori there is no reason to expect 
that the states with, say, 100 partons (quarks and gluons) are suppressed as compared those with the three valence 
quarks, the phenomenological success of the quark model allows one to hope that only a first few Fock components 
are really necessary. In hard exclusive reactions which involve a large momentum transfer to the nucleon, the 
dominance of valence states is widely expected and can be proven, at least within QCD perturbation theory II 1761 
11771 . 

The most general parametrization of the three-quark sector involves six scalar light-cone wave functions 11781 
11791 which correspond to different possibilities to couple the quark helicities Aj and orbital angular momentum 
L z to produce the helicity-1/2 nucleon state: Xi + X 2 + A3 + L z = 1/2. In particular if the quark helicities A; 
sum up to 1/2, then zero angular momentum is allowed, L = 0. The corresponding contribution can be written 
as EgHlTJl : 

fc 2 ))l4(z3, fc 3 )> - \4(x 2 , k 2 ))\4(x 3 , fc 3 )> . (6) 

Here ^^(xi, ki) is the light-cone wave function that depends on the momentum fractions Xi and transverse mo- 
menta h L of the quarks, \u\(xi, ki)) is a quark state with the indicated momenta and color index a, and e abc is the 
fully antisymmetric tensor; arrows indicate helicities. The integration measure is defined as 

[dx] — / dx\dx 2 dx 3 (5( Xj — l) , 
Jo 

J [d 2 k] = (16tt 3 )- 2 J dhdk 2 dk 3 5( ki) • (7) 

In hard processes the contribution of ^(xi, ki) is dominant whereas the other existing three-quark wave functions 
give rise to a power-suppressed correction, i.e. a correction of higher twist. 

The light-front description of a nucleon is very attractive for model building, but faces conceptual difficulties 
that do not allow the calculation of light-cone wave functions from first principles, at least at present. In particular 
there are subtle issues with renormalization and gauge dependence. An alternative approach has been to describe 
nucleon structure in terms of distribution amplitudes (DA) corresponding to matrix elements of nonlocal gauge- 
invariant light-ray operators. The classification of DAs goes in twist rather than number of constituents as for the 
wave functions. For example the leading-twist-three nucleon (proton) DA is defined by the matrix element 11801 : 

{Q\e^ k (u\{ ai n)cMM n )) ^ d lM\N{p)) = ~^f NP ■ n ft u\(p) f[dx] cp N ( Xi ) , (8) 
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where grt(-l-) = (1/2)(1 ± 75)5 are quark fields of given helicity, p^, p 2 = m 2 N , is the proton momentum, un(p) 
the usual Dirac spinor in relativistic normalization, an auxiliary light-like vector n 2 = and C the charge-con- 
jugation matrix. The Wilson lines that ensure gauge invariance are inserted between the quarks; they are not shown 
for brevity. The normalization constant fx is defined in such a way that 

[dx] ip N (xi) = 1. (9) 

In principle, the complete set of nucleon DAs carries full information on the nucleon structure, same as the com- 
plete basis of light-cone wave functions. In practice, however, both expansions have to be truncated and usefulness 
of a truncated version, taking into account either a first few Fock states or a few lowest twists, depends on the 
physics application. 

Using the wave function in Eq. (|6]l to calculate the matrix element in Eq. ^ it is easy to show that the DA 
tPN(xi) is related to the integral of the wave function ^f^(xi, ki) over transverse momenta, which corresponds to 
the limit of zero transverse separation between the quarks in the position space 11761 : 

fN(p)(PN(xi,fi) ~ J [d 2 k}^ N (xi,ki) . (10) 

\k\<n 

Thus, the normalization constant /jv can be interpreted as the nucleon wave function at the origin (in position 
space). 

Higher-twist three-quark DAs are related, in a loose sense, with similar integrals of the wave functions including 
extra powers of the transverse momentum, and with contributions of the other existing wave functions which 
correspond to nonzero quark orbital angular momentum. 

As always in a field theory, extraction of the asymptotic behavior produces divergences that have to be regulated. 
As the result, the DAs become scheme- and scale-dependent. In the calculation of physical observables this 
dependence is cancelled by the corresponding dependence of the coefficient functions. The DA Lp^(x i: fi) can be 
expanded in the set of orthogonal polynomials V n k(%i) defined as eigenfunctions of the corresponding one-loop 
evolution equation: 

00 N 

<pjv(£i,/*) = 120X1^2^3 

n=0 k=0 

where 

\dx\x 1 X2X' i Vnk{x i )T n ik< = NnAn'Skk' (12) 



and 



- ( M = ^)(^) W °". (B, 

2., 



Here Af n k are convention-dependent normalization factors, (3q — 11— |n/ and j n h the corresponding anomalous 



dimensions. The double sum in Eq. (Hi goes over all existing orthogonal polynomials V n k{xi), k — 0, . . . , n, 
of degree n. Explicit expressions for the polynomials V n k(xi) for n — 0, 1, 2 and the corresponding anomalous 
dimensions can be found in Ref. Ill 811 . 

In what follows we will refer to the coefficients c n k(fj,o) as shape parameters. The set of these coefficients to- 
gether with the normalization constant /jv(Mo) a * a reference scale /io specifies the momentum fraction distribution 
of valence quarks on the nucleon. They are nonperturbative quantities that can be related to matrix elements of 
local gauge-invariant three-quark operators (see below). 

In the last twenty years there had been mounting evidence that the simple-minded picture of a proton with the 
three valence quarks in an S-wave is insufficient, so that for example the proton spin is definitely not constructed 
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from the quark spins alone. If the orbital angular momenta of quarks and gluons are nonzero, the nucleon is 
intrinsically deformed. The general classification of three-quark light-cone wave functions with nonvanishing 
angular momentum has been worked out in Refs. 1178111791 . In particular the wave functions with L z = ±1 play 
a decisive role in hard processes involving a helicity flip, e.g. the Pauli electromagnetic form factor F 2 (Q 2 ) of 
the proton |182|. These wave functions are related, in the limit of small transverse separation, to the twist-four 
nucleon DAs introduced in Ref. 11 1 SOU : 

(0|e' ijfe (4( ai n)C^j(a 2 n))^4(a3n)|iV(p)) = ~p ■ nfu\, (p) J[dx] e -if»I>*«< 

■p-njt m N u\ (p) J [dx] e- ip - n ^ *««« 
f N ^ WW (x t )-\?^(x t ) , 
■p-ni m N u\{p) [ [dx] e - tp - n ^ x * a > 



' fc (uJ(a 1 n)C , ^7_ L ^(a 2 n)) ^^dl(a 3 n)\N(p)) 



(0\e^ k (uj ( ai n)Cfiiu]{a 2 n)) k (a 3 n)\N(p)) 



1 

~2 
x 

A 2 

12 



X ^ (Xi) , 



(14) 



where $ A ' (xi) and <3> 



4' WW (xi) are the so-called Wandzura-Wilczek contributions, which can be expressed 
in terms of the leading-twist DA (ppf(xi) 11811 . The two new constants A^ and A^ are defined in such a way that 
the integrals of the "genuine" twist-4 DAs $4, ^4, S4 are normalized to unity, similar to Eq. |9]). They are related 
to certain normalization integrals of the light-cone wave functions for the three-quark states with L z = ±1, see 
Ref. 032 for details. 

Light-cone wave functions and DAs of all baryons, including the nucleon resonances, can be constructed in a 
similar manner, taking into account spin and flavor symmetries. They can be constructed for all baryons of arbitrary 
spin without any conceptual complications, although it will become messy. The problem is only that "construct" 
means basically that one can enumerate different independent components and find their symmetries. To calculate 
them nonperturbatively is becoming increasingly difficult, however. This extension is especially simple for the 
parity doublets of the usual J p = | + octet since the nonlocal operators entering the definitions of nucleon DAs do 
not have a definite parity. Thus the same operators couple also to N* (1535) and one can define the corresponding 
leading-twist DA by the same expression as for the nucleon: 



e ijk (ui(ain)Criuj{a 2 n)) t4(a 3 n)\N*(p)) - \f N *p-nfiu\,{p)f[dx] e -^"5>^ 



(15) 



where, of course, p 2 — m 2 N ,. The constant fx* has a physical meaning of the wave function of A^*(1535) at the 



origin. The DA 0at» (x^ is normalized to unity (|9|) and has an expansion identical to ( 1 1 

00 N 



<Pn 



(xi,fi) = 120^1x2x3 c^ fc {n)Vnk{xi) 



(16) 



n=0 fc=0 



albeit with different shape parameters c^ k . 

Similar as for the nucleon, there exist three independent subleading twist-4 distribution amplitudes for the 
iV*(1535) resonance: , and . Explicit expressions are given in Ref. H141L 



V.B. Moments of distribution amplitudes from lattice QCD 

The normalization constants /, Ai, A2 and the shape parameters c n k are related to matrix elements of local 
three-quark operators between vacuum and the baryon state of interest, and can be calculated using lattice QCD. 
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FIG. 28. Probability amplitude /at, /jv* to find the three valence quarks in the nucleon and TV* (1535) at the same space-time 
point (wave function at the origin). 



Investigations of excited hadrons using this method are generally much more difficult compared to the ground 
states. On the other hand, the states of opposite parity can be separated rather reliably as propagating forwards 
and backwards in euclidian time. For this reason, for the time being we concentrate on the study of the ground 
state baryon octet J p = | + , and the lowest mass octet with negative parity, J p = | , N* (1535) being the prime 
example. 

Following the exploratory studies reported in Refs. 1U41I IT83 184 1 QCDSF collaboration is investing significant 
effort to make such calculations fully quantitative. The calculation is rather involved and requires the following 
steps: (1) Find lattice (discretized) operators that transform according to irreducible representations of spinorial 
group H(4); (2) Calculate non-perturbative renormalization constants for these operators; (3) Compute matrix 
elements of these operators on the lattice from suitable correlation functions, and (4) Extrapolate m T — > mP hys , 
lattice volume V — > oo and lattice spacing a — » 0. 

Irreducibly transforming H (4) multiplets for three-quark operators have been constructed in Ref. II 1 851 . Non- 
perturbative renormalization and one-loop scheme conversion factors RI-MOM— s- MS have been calculated in 
Ref. £186|. A consistent perturbative renormalization scheme for the three-quarks operators in dimensional regu- 
larization has been found 1 187 1 and the calculation of two-loop conversion factors using this scheme is in progress. 

The matrix elements of interest are calculated from correlation functions of the form (O a /3 1 (x)Af(y) T ), where 
Af is a smeared nucleon interpolator and O is a local three-quark operator with up to two derivatives, and applying 
the parity "projection" operator (l/2)(l±m74/-E) 11881 . In this way we get access to the normalization constants, 
the first and the second moments of the distribution amplitudes. Calculation of yet higher moments is considerably 
more difficult because one cannot avoid mixing with operators of lower dimension. 

The correlation functions were evaluated using Nf = 2 dynamic Wilson (clover) fermions on several lattices 
and a range of pion masses > 180 MeV. Our preliminary results for the wave functions the nucleon and 



iV*(1535) at the origin are summarized in Fig. 28 [189|. The extrapolation of the results for the nucleon to the 



physical pion mass and infinite volume as well as the analysis of the related systematic errors are in progress. An 



example of such an analysis is shown in Fig. 29 



This analysis will be done using one-loop chiral perturbation theory. The necessary expressions have been 



worked out in Ref. 11911 . Whereas the pion mass dependence of nucleon couplings is generally in agreement with 
expectations, we observe a large difference (up to a factor of three) in TV* (1535) couplings calculated with heavy 
and light pions: All couplings drop significantly in the transition region where the decay N* — > Ntt opens up. This 
effect can be due to the change in the structure of the wave function, but also to contamination of our N* (1535) 
results by the contribution of the ttN scattering state, or some other lattice artefact. This is one of the issues that 
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0.10 0.15 0.20 

FIG. 29. The chiral extrapolation of /jv to the physical (light) quark masses. The red points are lattice data and the blue points 
are corrected for finite volume effects. The green bands are the 1- and 1-a errors, respectively. The left-most black "data point" 
at the physical mass shows the recently updated estimate from QCD sum rule calculations 1 190 1. 
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FIG. 30. Leading-twist distribution amplitudes of the nucleon (left) and TV* (1535) (right) in barycentric coordinates xi + X2 + 
x 3 = 1. 



have to be clarified in future. 

We also find that the wave function of the TV* (1535) resonance is much more asymmetric compared to the 
nucleon: nearly 50% of the total momentum is carried by the u-quark with the same helicity. This shape is 
illustrated in Fig. 30 where the leading-twist distribution amplitudes of the nucleon (left) and N* (1535) (right) are 
shown in barycentric coordinates xi + X2 + x$ = 1; Xi are the momentum fractions carried by the three valence 
quarks. 

Our plans for the coming 2-3 years are as follows. The final analysis of the QCDSF lattice data using two 
flavors of dynamic fermions is nearly completed and in future we will go over to Nf = 2 + 1 studies, i.e. include 
dynamic strange quarks. The generation of the corresponding gauge configurations is in progress and first results 
are expected in one year from now. We will continue the studies of the lowest mass states in the J p = 1 /2 + and 
J p = l/2~ baryon octets. In particular the distribution amplitudes of the A and £ baryons will be studied for 
the first time. At a later stage we hope to be able to do similar calculations for the J p = 3/2 ± decuplets. We 
are working on the calculation of two-loop conversion factors RI-MOM— s- MS using the renormalization scheme 
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suggested in II 1871 and plan to employ them in the future studies. Main attention will be payed to the analysis of 
various sources of systematic uncertainties. With the recent advances in the algorithms and computer hardware 
the quark mass and finite volume extrapolations of lattice data have become less of a problem, which allows us to 
concentrate on more subtle issues. Our latest simulations for small pion masses make possible, for the first time, 
to study the transition region where decays of resonances, e.g. N* —> Nir, become kinematically allowed. We 
have to understand the influence of finite resonance width on the calculation of operator matrix elements and to 
this end plan to consider p-meson distribution amplitudes as a simpler example. We will also make detailed studies 
of meson (pion) distribution amplutides in order to understand better the lattice discretization errors and work out 
a concrete procedure to minimize their effect. The full programm is expected to last five years and is part of the 
proposal for the renewal of the Transregional Collaborative Research Centre (SFB/Transregio 55 "Hadron Physics 
with Lattice QCD") which will be submitted to the German Research Council (DFG) in April 2012. 



V.C. Light-cone distribution amplitudes and form factors 



The QCD approach to hard reactions is based on the concept of factorization: one tries to identify the short 
distance subprocess which is calculable in perturbation theory and take into account the contributions of large 
distances in terms of nonprerturbative parton distributions. 

The problem is that in the case of the baryon form factors the hard perturbative QCD (pQCD) contribution is 
only the third term of the factorization expansion. Schematically, one can envisage the expansion of, say, the Dirac 
electromagnetic nucleon form factor Fi(Q 2 ) of the form 



F^Q 2 ) ~ A(Q 2 ) + 



a s (Q 2 )\ BiQ 2 } , fa s (Q 2 ) 
Q 2 



C_ 



(17) 



where C is a constant determined by the nucleon DAs, while A(Q 2 ) and B(Q 2 ) are form-factor-type functions 
generated by contributions of low virtualities, see Fig. 31 The soft functions A(Q 2 ) and B(Q 2 ) are purely non- 
perturbative and cannot be further simplified e.g. factorized in terms of DAs. In the light-cone formalism, they 
are determined by overlap integrals of the soft parts of hadronic wave functions corresponding to large transverse 
separations. Various estimates suggest that A{Q 2 ) < 1/Q 6 , B(Q 2 ) < 1/Q 4 and at very large Q 2 they are further 
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FIG. 31. Structure of QCD factorization for baryon form factors. 



suppressed by the Sudakov form factor. To be precise, in higher orders in a s (Q) there exist double-logarithmic 
contributions ~ 1/Q 4 II 1921 that are not factorized in the standard manner; however, also they are suppressed by 
the Sudakov mechanism |192, 193]. Thus, the third term in ( 17 1 is formally the leading one at large Q 2 to power 
accuracy. 

The main problem of the pQCD approach I176 J1771 is a numerical suppression of each hard gluon exchange 
by the a s /Tr factor which is a standard perturbation theory penalty for each extra loop. If, say, a s /ir ~ 0.1, the 
pQCD contribution to baryon form factors is suppressed by a factor of 100 compared to the purely soft term. As 
the result, the onset of the perturbative regime is postponed to very large momentum transfers since the factorizable 
pQCD contribution 0(1/Q 4 ) has to win over nonperturbative effects that are suppressed by extra powers of 1/Q 2 , 
but do not involve small coefficients. There is an (almost) overall consensus that "soft" contributions play the 
dominant role at present energies. Indeed, it is known for a long time that the use of QCD-motivated models for the 
wave functions allows one to obtain, without much effort, soft contributions comparable in size to experimentally 
observed values. Also models of generalized parton distributions usually are chosen such that the experimental 
data on form factors are described by the soft contributions alone. A subtle point for these semi-phenomenological 
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approaches is to avoid double counting of hard rescattering contributions "hidden" in the model-dependent hadron 
wave functions or GPD parametrizations. 

One expects that the rapid development of lattice QCD will allow one to calculate several benchmark baryon 
form factors to sufficient precision from first principles. Such calculations are necessary and interesting in its own 
right, but do not add to our understanding of how QCD actually "works" to transfer the large momentum along 
the nucleon constituents, the quarks and gluons. The main motivation to study "hard"processes has always been 
to understand hadron properties in terms of quark and gluon degrees of freedom; for example, the rationale for 
the continuing measurements of the total inclusive cross section in deep inelastic scattering is to extract quark and 
gluon parton distributions. Similar, experimental measurements of the electroproduction of nucleon resonances at 
large momentum transfers should eventually allow one to get insight in their structure on parton level, in particular 
momentum fraction distributions of the valence quarks and their orbital angular momentum encoded in DAs, and 
this task is obscured by the presence of large "soft" contributions which have to be subtracted. 

Starting in Ref. II 1941 and in subsequent publications we have been developing an approach to hard exclusive 
processes with baryons based on light-cone sum rules (LCSR) B1951 11961 . This technique is attractive because 
in LCSRs "soft" contributions to the form factors are calculated in terms of the same DAs that enter the pQCD 
calculation and there is no double counting. Thus, the LCSRs provide one with the most direct relation of the 
hadron form factors and distribution amplitudes for realistic momentum transfers of the order of 2 — 10 GeV 2 that 
is available at present, with no other nonperturbative parameters. It is also sufficiently general and can be applied 
to many hard reactions. 

The basic object of the LCSR approach is the correlation function 



dxe^ x (N*(P)\T{r,(0)j(x)}\0) 



(18) 



in which j represents the electromagnetic (or weak) probe and rj is a suitable operator with nucleon quantum 
numbers. The nucleon resonance in the final state is explicitly represented by its state vector \N*(P)), see a 
schematic representation in Fig. 32 When both the momentum transfer q 2 = —Q 2 and the momentum (P 1 ) 2 = 




FIG. 32. Schematic structure of the light-cone sum rule for electroproduction of nucleon resonances. 

(P + q) 2 flowing in the r\ vertex are large and negative, the asymptotic of the correlation function is governed by 
the light-cone kinematics x 2 — > and can be studied using the operator product expansion (OPE) T{rj(0)j(x)} ~ 
Ci(x)Oi(0) on the light-cone x 2 = 0. The x 2 -singularity of a particular perturbatively calculable short-distance 
factor Ci{x) is determined by the twist of the relevant composite operator Ou whose matrix element (N* |Oj(0)|0) 
is given by an appropriate moment of the N* DA. Next, one can represent the answer in form of the dispersion 
integral in (P 1 ) 2 and define the nucleon contribution by the cutoff in the quark-antiquark invariant mass, the so- 
called interval of duality so (or continuum threshold). The main role of the interval of duality is that it does not 
allow large momenta \k 2 \ > so to flow through the 77-vertex; to the lowest order 0{oP s ) one obtains a purely soft 
contribution to the form factor as a sum of terms ordered by twist of the relevant operators and hence including both 
the leading- and the higher-twist nucleon DAs. Note that, in difference to the hard mechanism, the contribution of 
higher-twist DAs is only suppressed by powers of the interval of duality sq ~ 2 GeV 2 (or by powers of the Borel 
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Q 2 

FIG. 33. The LCSR calculation for the helicity amplitudes A\ /2 (Q 2 ) and Si /2 (Q 2 ) for the electroproduction of the N* (1535) 
resonance using the lattice results for the lowest moments of the AT* (1535) DAs. The curves are obtained using the central 
values of the lattice parameters, and the shaded areas show the corresponding uncertainty. Figure taken from Ref. 11411 . 
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FIG. 34. LCSR results for the magnetic proton form factor (normalized to the dipole formula) for a realistic model of nu- 
cleon distribution amplitudes 1 197 1 . Left panel: Leading order (LO); right panel: next-to-leading order (NLO) for twist-three 
contributions. Figure adapted from Ref. [ 1 98 1 . 

parameter if one applies some standard QCD sum rule machinery), but not by powers of Q 2 . This feature is in 
agreement with the common wisdom that soft contributions are not constrained to small transverse separations. 

We stress that LCSRs are not based on any nonperturbative model of the nucleon structure, but rather present a 
relation between the physical observables (form factors) and baryon wave functions at small transverse separation 
(distribution amplitudes). 

Historically, LCSRs were developed in Refs. 11 1 951 11961 in an attempt to overcome difficulties of the Shifman- 
Vainstein-Zakharov QCD sum rule approach 1 199 1 for exclusive processes dominated by the light-cone kinematics. 
In the last 20 years LCSRs have been applied extensively to the exclusive B-decays and remain to be the only 
nonperturbative technique that allows one to calculate the corresponding form factors directly at large recoil. In 
fact the value of the CKM matrix element V u b quoted by the Particle Data Group as the one extracted from 
exclusive semileptonic decay B — > irlvi is largely based on the recently updated LCSR calculations of the form 
factor f? (0) B2001 12011 (although the lattice QCD calculations have become competitive). Another important 
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FIG. 35. LCSR results for the electric to magnetic proton form factor ratio for a realistic model of nucleon distribution am- 
plitudes 11971 . Left panel: Leading order (LO); right panel: next-to-leading order (NLO) for twist-three contributions. Figure 
adapted from Ref. fT98l . 



application of LCSRs was for calculation of the electromagnetic pion form factor. More references and further 
details can be found in the review articles [202, 2031. 

LCSRs for meson form factors have achieved a certain degree of maturity. One lesson is that they are fully 
consistent with pQCD and factorization theorems. In particular the LCSRs also contain terms generating the 
asymptotic pQCD contributions. In the pion case, it was explicitly demonstrated that the contribution of hard 
rescattering is correctly reproduced in the LCSR approach as a part of the 0(a s ) correction. It should be noted 
that the diagrams of LCSR that contain the "hard" pQCD contributions also possess "soft" parts, i.e., one should 
perform a separation of "hard" and "soft" terms inside each diagram. As a result, the distinction between "hard" 
and "soft" contributions appears to be scale- and scheme-dependent. Most of the LCSRs for meson decays have 
been derived to the next-to-leading-order (NLO) accuracy in the strong coupling. The first NLO LCSR calculations 
were done in 1997-1998 and since then the NLO accuracy has become standard in this field. The size of NLO 
corrections depends on the form factor in question but typically is of the order of 20%, for the momentum transfers 
of interest. 

Derivation of LCSRs for exclusive reactions involving baryons is, conceptually, a straightforward generalization 
of the LCSRs for mesons. On the other hand, there are a few new technical issues that had to be resolved, and 
also the calculations become much more challenging. The development so far was mainly to explore the existing 
possibilities and identify potential applications. Following the first application to the electromagnetic and axial 
form factors of the nucleon in Refs. II 1941 1 1971 , LCSRs have been considered for the 7*7V — > A transition |204|, 
heavy baryon decays (see |205| and references therein) and various transitions between baryons in the octet and 
the decuplet (e.g. 020611 ). In the work 11411 we have suggested to use the same approach to the study of electro- 
production of resonances at large momentum transfers and in particular TV* (1535). Since the structure of sum 
rules for the nucleon elastic form factors and electroproduction of N* (1535) is very similar, the difference in form 
factors should expose directly the difference in the wave functions, which is of prime interest. The results for the 
helicity amplitudes Ax^iQ 2 ) and <Si/2(Q 2 ) using the lattice results for the lowest moments of the N* (1535) DAs 
appear to be in a good agreement with the existing data, see Fig. [33] 

All existing LCSRs for baryons are written to the leading order in the strong coupling which corresponds, 
roughly speaking, to the parton model level description of deep-inelastic scattering. Combined with realistic 
models of DAs the existing sum rules yield a reasonable description of the existing data to the expected 30-50% 
accuracy. In order to match the accuracy of the future experimental data and also of the next generation of lattice 
results, the LCSRs will have to be advanced to include NLO radiative corrections, as it has become standard for 
meson decays. 

The first step towards LCSRs to the NLO accuracy was done in Ref. [198] where the 0(a s ) corrections are 
calculated for the (leading) twist-three contributions to the sum rules for electromagnetic (elastic) nucleon form 



factors derived in M194II1971I . The results are shown in Fig. 34 and Fig. 35 
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The NLO corrections are large and their effect increases with Q 2 which may be counterintuitive. This behavior 
is, however, expected on general grounds because the leading regions for large momentum transfers corresponding 
to the ERBL (Efremov-Radyushkin-Brodsky-Lepage) collinear factorization appear at the NNLO level only, i.e. 
0(a 2 s ). The corrections for the Ge/Gm ratio are larger than for the magnetic form factor Gm itself, which is 
again expected since the electric form factor suffers from cancellations between chirality-conserving and chirality- 
violating contributions. 

Large NLO corrections can be compensated by the change in the nucleon DA, similar as it happens with parton 
distributions — e.g. the small-x behavior of the LO and NLO gluon distribution is very different — but such an 
analysis would so far be premature since NLO corrections have not been calculated so far for the contributions of 
twist-four DAs that take into account the effects of orbital angular momentum. 

In addition, it is necesary to develop a technique for the resummation of "kinematic" corrections to the sum rules 
that are due to nonvanishing masses of the resonances. The corresponding corrections to the total cross section 
of the deep-inelastic scattering are known as Wandzura-Wilczek corrections and can be resummed to all orders in 
terms of the Nachtmann variable; we are looking for a generalization of this method to non-forward kinematics 
which is also important in a broader context 12071 . 

With these improvements, we expect that the LCSR approach can be used to constrain light-cone DAs of the 
nucleon and its resonances from the comparison with the electroproduction data. These constraints can then be 
compared with the lattice QCD calculations. In order to facilitate this comparison, a work is in progress to derive 
general expressions for the necessary light-cone sum rules to the NLO accuracy. The project is to have the LCSRs 
available as a computer code allowing one to calculate elastic electromagnetic and axial form factors and also a 
range of transition form factors involving nucleon resonances from a given set of distribution amplitudes. Although 
gross features of the wave functions of resonances can definitely be extracted from such an analysis, the level of 
details "seen" in sum rule calculations will have to be tested on case by case basis. For this reason we are also 
working on similar calculations for the "gold-plated" decays like 7* — > 777, 7* — > 7/7, see [208], where the 
theoretical uncertainties are expected to be small. 
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VI. QUARK-HADRON DUALITY AND TRANSITION FORM FACTORS 

VI.A. Historical perspective 

Understanding the structure and interactions of hadrons at intermediate energies is one of the most challenging 
outstanding problems in nuclear physics. While many hadronic observables can be described in terms of effective 
meson and baryon degrees of freedom at low energies, at energies 3> the nucleon mass M perturbative QCD has 
been very successful in describing processes in terms of fundamental quark and gluon (parton) constituents. 
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FIG. 36. Proton F 2 structure function data from Jefferson Lab Hall C I209H211I . SLAC EHHED, and NMC EH) at 
Q 2 = 0.5, 1.5,3 and 5.5 GeV 2 , compared with an empirical fit 12151 to the transverse and longitudinal resonance cross 
sections (solid), and a global fit to DIS data (dashed). (Figure from Ref. 12161 .) 



A connection between the low and high energy realms is realized through the remarkable phenomenon of quark- 
hadron duality, where one often finds dual descriptions of observables in terms of either explicit partonic degrees 
of freedom, or as averages over hadronic variables. In principle, with access to complete sets of either hadronic or 
partonic states, the realization of duality would be essentially trivial, effectively through a simple transformation 
from one complete set of basis states to another. In practice, however, at finite energies one is typically restricted 
to a limited set of basis states, so that the experimental observation of duality raises the question of not why duality 
exists, but rather how it arises where it exists, and how we can make use of it. 

Historically, duality in the strong interaction physics represented the relationship between the description of 
hadronic scattering amplitudes in terms of s-channel resonances at low energies, and i-channel Regge poles at 
high energies H219I1 . The merger of these dual descriptions at intermediate energies remained a prized goal of 
physicists in the decade or so before the advent of QCD. Progress towards synthesizing the two descriptions was 
made with the development of finite energy sum rules (FESRs) [220, 221], 

dv v n 3m A{y, t) = [ dv v n A asy (v, t) [FESR] (19) 
Jo 

relating the imaginary part of the amplitude A at finite energy to the asymptotic high energy amplitude „4 a sy> 
where s, t and u are the usual Mandelstam variables and v = (s — u)/A. The asymptotic amplitude „4 asy is 
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FIG. 37. Ratio of proton F% structure functions integrated over specific resonance regions, relative to the global fit of parton 
distributions from Alekhin et al. I217II218I . (Figure from Ref. 120^1 .) 



then extrapolated into the v < v max region and compared with the measured amplitude A through Eq. ( 19 1. The 
assumption made here is that beyond some maximum energy v > v max the scattering amplitude can be represented 
by its asymptotic form, calculated within Regge theory. 



The FESRs are generalizations of superconvergence relations in Regge theory relating dispersion integrals over 
the amplitudes at low energies to high-energy parameters. They constitute a powerful tool allowing one to use 
experimental information on the low energy cross sections for the analysis of high energy scattering data. Con- 
versely, they can be used to connect low energy parameters (such as resonance widths and couplings) to parameters 
describing the behavior of cross sections at high energy. It was in the context of FESRs, in fact, that the early ex- 
pressions of Bloom-Gilman duality were made in the early 1970s [222 223 1, suitably extended to lepton scattering 
kinematics. 
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VLB. Duality in nucleon structure functions 



One of the most dramatic realizations of duality in nature is in inclusive electron-nucleon scattering, usually 
referred to as "Bloom-Gilman" duality, where structure functions averaged over the resonance region are found 
to be remarkably similar to the leading twist structure functions describing the deep-inelastic scattering (DIS) 
continuum [ 209 216 222-226 1. As Fig. 36 illustrates, the resonance data are seen to oscillate around the scaling 
curve and slide along it with increasing Q\ 

An intriguing feature of the lepton scattering data is that the duality appears to be realized not just over the entire 
resonance region as a whole, W < 2 GeV, where W 2 = M 2 + Q 2 (l — x)/x, but also in individual resonance 
regions. This is illustrated in Fig. [37] for the ratios of structure functions integrated over specific intervals of W 
at fixed Q 2 , with the 1st, 2nd, 3rd and 4th resonance regions defined by 1.3 < W 2 < 1.9 GeV 2 , 1.9 < W 2 < 
2.5 GeV 2 , 2.5 < W 2 < 3.1 GeV 2 , and 3.1 < W 2 < 3.9 GeV 2 , respectively. The "DIS" region in Fig 
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is 

defined to be 3.9 < W 2 < 4.5 GeV 2 . In all cases the duality is realized at the < 10 — 15% level, suggesting that 
Bloom-Gilman duality exists locally as well as globally. 

Understanding the microscopic origin of quark-hadron duality has proved to be a major challenge in QCD. 
Until recently the only rigorous connection with QCD has been within the operator product expansion (OPE), 
in which moments (or a; -integrals) of structure functions are expanded as a series in inverse powers of Q 2 . The 
leading, 0(1) term is given by matrix elements of (leading twist) quark-gluon bilocal operators associated with 
free quark scattering, while the 0(1/Q 2 ) and higher terms correspond to nonperturbative (higher twist) quark- 
gluon interactions. In the language of the OPE, duality is then synonymous with the suppression of higher twist 
contributions to the moments B2271 . 

This close relationship between the leading twist cross sections and the resonance-averaged cross sections sug- 
gests that the total higher twist contributions are small at scales Q 2 ~ 1 GeV 2 . This implies that, on average, 
nonperturbative interactions between quarks and gluons are not dominant at these scales, and that a highly non- 
trivial pattern of interferences emerges between the resonances (and the nonresonant background) to effect the 
cancellation of the higher twist contributions. The physics of parton distributions and nucleon resonances is there- 
fore intimately connected. In fact, in the limit of a large number of colors, the spectrum of hadrons in QCD is one 
of infinitely narrow resonances [228 1, which graphically illustrates the fact that resonances are an integral part of 
scaling structure functions. 

The phenomenological results raise the question of how can a scaling structure function be built up entirely 
from resonances, each of whose contribution falls rapidly with Q 2 12291 ? A number of studies using various 
nonperturbative models have demonstrated how sums over resonances can indeed yield a Q 2 independent function 
(see Ref. [226 1 for a review). The key observation is that while the contribution from each individual resonance 
diminishes with Q 2 , with increasing energy new states become accessible whose contributions compensate in such 
a way as to maintain an approximately constant strength overall. At a more microscopic level, the critical aspect 
of realizing the suppression of the higher twists is that at least one complete set of even and odd parity resonances 
must be summed over for duality to hold B2301 . Explicit demonstration of how this cancellation takes place was 
made in the SU(6) quark model and its extensions [230 -2321. 

One of the ultimate goals of duality studies is to determine the extent to which resonance region data can be used 
to learn about leading twist structure functions. At present, most global analyses of parton distribution functions 
impose strong cuts on Q 2 and W 2 for lepton scattering data in order to exclude the region where higher twists 
and other subleading effects are important. By relaxing the cuts to just inside the traditional resonance region, 
W > 1.7 GeV, the CTEQ-Jefferson Lab (CJ) collaboration could increase the statistics of the DIS data by a 
factor ~ 2 [233, 234 1! Not only were the fits found to be stable with the weaker cuts, the larger database led to 
significantly reduced errors, up to 40-60% at large x, where data are scarce. Future plans include extending these 
cuts to even lower values of W, which demands better understanding of the resonance region and the procedures 
for systematically averaging over the resonance structure functions. 

The determination of parton distributions at large x is vital not just for understanding the dynamics of valence 



quarks in the nucleon B235I |236], which are currently obscured by nuclear corrections in deuterium DIS data 
needed to extract the structure function of the free neutron. It is also critical in applications to experiments at high 
energy colliders, where uncertainties at large x in the d quark distribution in particular feeds down to lower x at 
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higher Q 2 [237| and can have important consequences for searches for new particles, such as W' and Z' bosons 
|238|. Thus in an indirect way, better knowledge of the nucleon resonance region can have a profound impact on 
physics at the LHC! 



VI.C. Duality in inclusive meson production 

Extending the concept of duality to less inclusive reactions, we can ask whether the semi-inclusive production 
of mesons displays a similar relation between partonic and resonance-based descriptions. Such studies have only 
recently been performed, for ratios of semi-inclusive tt + to tt~ cross sections measured at Jefferson Lab as a 
function of z = E^/v, where E n is the pion energy and v is the energy transfer to the target I239L 

The data displayed a smooth behavior in z, consistent with earlier observations at higher energies at CERN 
1 240], prompting suggestions that factorization of semi-inclusive cross sections into scattering and fragmentation 
sub-processes may hold to relatively low energies. Such factorization was found in fact in simple quark models by 
explicitly summing over N* resonances in the s-channel of 7* iV —> ttN scattering 12321 . 

At the quark level, the (normalized) semi-inclusive cross section for the production of pions from a nucleon 
target can be factorized (at leading order in a s ) into a product of a parton distribution function describing the hard 
scattering from a parton in the target, and the probability of the struck parton fragmenting into a specific hadron, 

^ c< £ e 2 q q N (x) D*(z), (20) 

where e q is the quark charge, and D q is the fragmentation function for quark q to produce a pion with energy 
fraction z. As pointed out by Close and Isgur M230I . duality between structure functions represented by (incoher- 
ent) parton distributions and by a (coherent) sum of squares of form factors can be achieved by summing over 
neighboring odd and even parity states. In the SU(6) model this is realized by summing over states in the 56 + 
(L = 0, even parity) and 70~ (L = 1, odd parity) multiplets, with each representation weighted equally. 




FIG. 38. Momentum spectrum of produced hadrons in the inclusive hadron production reaction 7*iV — > MX. From 
Ref. 12411 . 

The pion production cross sections at the hadronic level are constructed by summing coherently over excited 
nucleon resonances (N*) in the s-channel intermediate state and in the final state (iVJ) of jN — > N* — > irN^, 
where both N£ and iV^ belong to the 56 + and 70~ multiplets. Within this framework, the probabilities of the 
jN — > 7riV| transitions can be obtained by summing over the intermediate states N± spanning the 56 + and 70~ 
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multiplets, with the differential cross section 



da 

dxdz 



oc 



E 

N* 



N: 



(21) 



Here R 



■yN—yN* 



is the — > N* transition form factor, which depends on the masses of the virtual photon and 

where M| is the invariant 



excited nucleon (M*), and T>n^->-n^tt 1S a function representing the decay N£ — > ttN^ 
mass of the final state N%. 

Summing over the N£ states in the 56 + and 70 _ multiplets, one finds ratios of unpolarized ir~ to ir + semi- 
inclusive cross sections consistent with the parton model results for ratios of parton distributions satisfying SU(6) 
symmetry [230-232]. Duality was also found to be realized in more realistic scenarios with broken SU(6) sym- 



metry, with sums over resonances able to reproduce parton model semi -inclusive cross section ratios [232). The 
absence of strong resonant enhancement on top of the smooth background is indeed one of the notable features of 
the Jefferson Lab Hall C data 1123911 , in accord with expectations from duality. 



VI.D. Exclusive-inclusive connection 



The general folklore in hadronic physics is that duality works more effectively for inclusive observables than for 
exclusive, due to the presence in the latter of fewer hadronic states over which to average. For exclusive processes, 
such as the production of a meson M in coincidence with and a baryon B, eN — > eM B, duality may be more 
speculative. Nevertheless, there are correspondence arguments formulated long ago which relate the exclusive 
cross sections at low energy to inclusive production rates at high energy. The exclusive-inclusive connection dates 
back to the early dates of DIS and the discussion of scaling laws in high energy processes. Bjorken & Kogut [241 ] 
proposed the correspondence relations by demanding the continuity of the dynamics as one goes from one (known) 
region of kinematics to another (which is unknown or poorly known). 




FIG. 39. Deeply virtual (nonforward) Compton scattering at the hadronic level, with excitation of nucleon resonances N* in 
the intermediate state. (Adapted from Ref. 12421 .) 



For processes such as 7* N — > AIB, the correspondence principle relates properties of exclusive (resonant) final 
states with inclusive particle spectra for the corresponding reaction 7* N — > MX. This is illustrated in Fig. 38 for 
a typical inclusive momentum spectrum Ed 3 a/dp 3 , where E and p are the energy and momentum of the observed 
final state particle AI. As p increases, the inclusive continuum gives way to to the region dominated by resonances. 
The correspondence argument postulates that the resonance contribution to the cross section should be comparable 
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to the continuum contribution extrapolated from high energy into the resonance region 

dp 



d P E ^ 



incl res d Pr 



(22) 

1 



where the integration region over the inclusive cross section includes contributions up to a missing mass Mx, 
with Pmin = Pmax ~ M| /4p max . The correspondence relation (??) is another manifestation of the FESR in 
Eq. ( p~9] >, in which the cross section in the resonance region for p m j n < p < p max is dual to the high-energy cross 
section extrapolated down to the same region. 

The inclusive cross section d 3 cr/dp 3 is generally a function of the longitudinal momentum fraction x, the trans- 
verse momentum pr, and the invariant mass squared s, 



Ed 3 (T 

a dp 



, -f(x,p' T ,sQ z ). (23) 



At large s or large Q 2 this effectively reduces to a function of only x and p\, 

f{x,p 2 T ,sQ 2 ) -> f(x,p%) , s^oo. (24) 

The continuity relation (??) implies that there should be no systematic variation of either side of the equation with 
external parameters. 

Applications of the exclusive-inclusive correspondence have also been made to real Compton scattering cross 
sections from the proton at large center of mass frame angles 12431 . as well as to hard exclusive pion photoproduc- 
tion [244-246], and more recently to deeply virtual Compton scattering, ep —> ejp 12421 . The latter in particular 
used a simple model with scalar constituents confined by a harmonic oscillator potential to show how sums over 



intermediate state resonances, Fig. 39 lead to destructive interference between all but the elastic contribution, and 
the emergence of scaling behavior for the associated generalized parton distributions (GPDs). 

Future work will build on these exploratory studies, generalizing the calculations to include spin- 1/2 quarks and 
non-degenerate multiplets, as well as incorporating nonresonant background within same framework 1 247] 12481 . 
Extension to flavor non-diagonal transitions will also establish a direct link between transition form factors and the 
GPDs, with duality providing a crucial link between the hadronic and partonic descriptions. 
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VII. THE N* ELECTROCOUPLING INTERPRETATION WITHIN THE FRAMEWORK OF CONSTITUENT 

QUARK MODELS 

VILA. Introduction 

The study of the electromagnetic excitation of the nucleon resonances is expected to provide a good test for our 
knowledge concerning the internal structure of baryons. From a fundamental point of view, the description of the 
resonance spectrum and excitation should be performed within a Quantum ChromoDynamics (QCD) approach, 
which, however, does not allow up to now to extract all the hadron properties in a systematic way. Therefore, one 
has to rely on models, such as the Constituent Quark Models (CQM). In CQMs quarks are considered as effective 
internal degrees of freedom and can acquire a mass and a finite size. Phenomenological results of these models 
provide useful constraints based on experimental data for the development of QCD-based approaches such as 
LQCD and DSEQCD. For instance the choice of basis configurations in the recent LQCD studies of N* spectrum 
[ 3 1 was motived by quark model results. 

In the following we report some results of recent approaches using the constituent quark idea in the framework 
of various light front (LF) formulations of the quark wave function (Sees. II-IV) and a discussion on the use 
of CQM for the interpretation of resonance electrocouplings at high Q 2 with particular attention to some future 
perspectives (Sec. V). 



VII.B. Covariant quark-diquark model for the N and TV* electromagnetic transition form factors 

The study of hadron structure using the fundamental theory, Quantum ChromoDynamics, can in practise be 
done only in the large Q 2 regime or, by means of lattice simulations, in the unphysical quark masses regime 0. 
For this reason one has to rely on effective descriptions either with the degrees of freedom of QCD (quarks and 
gluons) within the Dyson-Schwinger framework 0, or in terms of the degrees of freedom observed at low Q 2 , the 
meson cloud and the light baryon core, using a dynamical coupled-channels reaction (dynamical models or DM) 
framework 151 12491 . The DSEQCD helps to understand the transition between the perturbative regime of QCD 
and the low Q 2 regime, where the quarks acquire masses and structure dynamically due to the gluon dressing, 
although the meson degrees of freedom are not included till the moment 0. Dynamical models, on the other 
hand, help to explain the transition between the low Q 2 picture, in terms of a finite size baryon and the surrounding 
meson cloud, and the intermediate region when Q 2 > 2 GeV 2 , where the baryon core effects become increasingly 
important [ 249 1 . To complete the picture a parameterization of the structure of the baryon core is required, and a 
possibility is to use the meson-baryon dressing model to extract from the data the contributions of the core, that 
can be interpreted as a 3-valence quark system Il62ll70l . 

Alternative descriptions comprise effective chiral perturbation theory, that can be used to interpolate lattice QCD 
results but is restricted to the low Q 2 regime, perturbative QCD that works only at very large Q 2 with a threshold 
that is still under discussion, QCD sum rules and constituent quark models that can include also chiral symmetry 
and/or unquenched effects 0. 

CQMs include the gluon and quark-antiquark polarization in the quark substructure (that also generates the 
constituent quark mass) with effective inter-quark interactions [5|. There are different versions according to the 
inter-quark interaction potential and the kinematic considered (nonrelativistic, or relativistic). Among the relativis- 
tic descriptions there are, in particular, different implementations of relativity based on the Poincare invariance 
0. 

We discuss now with some detail the covariant quark-diquark model, also known as the spectator quark model, 
and present some of its results. Contrarily to other CQMs, this model is not based on a wave equation determined 
by some complex and nonlinear potential. For that reason, the model is not used to predict the baryonic spectrum. 
Instead, the wave functions are built from the baryon internal symmetries only, with the shape of the wave functions 
determined directly by the experimental data, or lattice data for some ground state systems [250 1. 

In the covariant spectator theory (CST) 12511 the 3-body baryon systems are described in terms of a vertex 
function T where 2 quarks are on-mass-shell [252-254]. In this approach confinement ensures that the vertex T 
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FIG. 40. Nucleon form factors. Model II of Ref. [254], Left panel: n P Ge p /Gm p ratio, including the Jefferson Lab data. Right 
panel: neutron electric form form factor. 



vanishes when the 3 quarks are simultaneously on-mass-shell, and the singularities associated with the propagator 
of the off-mass-shell quark are canceled by the vertex T [252, 253]. The baryon state can then be described by a 
wave function ^f(P, k) = {m q — ft — ie)^ 1 T(P, k), where P is the baryon momentum, m q the quark mass and k 
the quark four- momentum [253 2541. 

The CST formulation is motivated by the fact that in impulse approximation only one quark interacts with the 
photon, while the two other quarks are spectators. Therefore, by integrating over the relative momentum of these 
two quarks, one can reduce the 3-quark system to a quark-diquark system, where the effective diquark has an 
averaged mass mj 12531 12541 . In these conditions the baryon is described by a wave function for the quark- 
diquark, with individual states associated with the internal symmetries (color, flavor, spin, momentum, etc.). The 
electromagnetic interaction current is given in impulse approximation by the coupling of the photon with the off- 
mass-shell quark, while the diquark acts as a spectator on-mass-shell particle [250 254 255 1 . 

The photon-quark interaction is parameterized by using the vector meson dominance (VMD) mechanism, based 
on a combination of two poles associated with vector mesons: a light vector meson (mass m v = m p ~ m u ) 
and an effective heavy meson with mass M/, = 2M, where M is the nucleon mass, which modulates the short 
range structure [250 254, 255 1. The free parameters of the current were calibrated for the SU(3) sector by nucleon 
electromagnetic form factor data [ 254 1 and with lattice QCD simulations associated with the baryon decuplet 12551 . 
A parameterization based on VMD has the advantage in the generalization to the lattice QCD regime 0255H2581 
and also for the time-like region (Q 2 < 0) 12591 . 

The covariant spectator quark model was applied to the description of the nucleon elastic form factors using a 
simple model where the quark-diquark motion is taken in the S-state approximation [254|. The nucleon data were 
used to fix the quark current as well as the radial wave function [254]. A specific model with no explicit pion cloud 
effects, except the effects included in the VMD parameterization is presented in Fig. |40] This parameterization, 
based only on the valence quark degrees of freedom, was extended successfully for the nucleon on the lattice 
regime 1741 . 

The model was also applied to the first nucleon resonance the A(1232), in particular to the 7^ — > A(1232) 
transition. Within a minimal model where the A is described as an S-state of 3-quarks with the total spin and 
isospin 3/2, one obtains, for dominant transition form factor G* M (Q) < 2.071 < 2.07, where I < 1, is the overlap 
integral between the nucleon and A radial wave functions (both are S-states) in the limit Q 2 = B260I . This simple 
relation, which is a consequence of the normalization of the nucleon and A quark wave functions, illustrates the 
incapability in describing the 7^ — !• A(1232), with quark degrees of freedom only, since the experimental result 
is G* M (Q) ~ 3. The discrepancy, is common to all constituent quark models, and is also a manifestation of the 
importance of the pion excitation which contributes with about 30-40% of the strength of the reaction [5 70, 249]. 
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FIG. 41. Nucleon electromagnetic transition for spin 3/2 resonances. Left panel: G* M / (3Gr>) (Gd is the nucleon dipole form 
factor) for the 7iV A(1232) reaction 12611 . Right panel: G* M for the jN -> A(1600) reaction 12621 . In both cases the 
dashed line gives the valence quark contribution and the solid line the full result. 



The model can however explain the quark core contribution in the transition, as extracted from the data using the 
EBAC model ll70l . when the meson-baryon cloud is subtracted [260 1 . The comparison of the model with the EBAC 
estimate is presented in Fig. 41 (left panel, dashed line), and also with the G* M data, when meson-baryon cloud is 



included (solid line). The model was also extended successfully to the reaction in the lattice regime 12561 12571 . 
The description of the quadmpole form factors G* E (electric) and G* c (Coulomb) is also possible once small D 
state components are included [257 261 ]. In that case, the lattice QCD data can be well described by an extension 
of the model with an admixture of D-states less than 



1% [257 1, but the experimental data are fairly explained 
only when the meson-baryon cloud and valence quark degrees of freedom are combined 12571 . Finally, the model 
was also applied to the first radial excitation of the A(1232), the A(1600) resonance 12621 . In this case no extra 
parameters are necessary, and the meson-baryon cloud effects are largely dominant at low Q 2 . The results for G* M 
are presented in the right panel of Fig. 41 In both systems the valence quark effects are dominant for Q 2 > 2 



GeV 2 . 

The model was also extended to the spin 1/2 state A(1440) (Roper), interpreted as the first radial excitation 
of the nucleon [263 1. The A?" (1440) shared with the nucleon the spin and isospin structure, differing in the radial 
wave function. Under that assumption we calculated the transition form factors for the jN —> A(1440) reaction 
based exclusively on the valence quark degrees of freedom 12631 . As an example, we present the Dirac-type 
form factor in Fig. 42 (left panel). The model is also consistent with the lattice data [ 263 1 . The covariant 



spectator quark model was also applied to the chiral partner of the nucleon A(1535) (negative parity) under two 
approximations: a pointlike diquark and a quark core restricted to spin 1/2 states [264 1 . Under these approximations 
the jN — > JV(1535) transition form factors were calculated for the region Q 2 ^> 0.23 GeV 2 [264 1. The result for 
Fi is presented in Fig. 42 (right panel). In both reactions the results are consistent with the data for Q 2 > 1.5 



GeV 2 [263 , 264 1, except for F 2 * for the reaction with A(1535). Our results support the idea that the valence quark 
dominance for the intermediate and high Q 2 region, but also the necessity of the meson excitations for the lower 
Q 2 region (Q 2 < 2 GeV 2 ). The form factor F 2 * for the jN — > A(1535) reaction is particularly interesting from 
the perspective of a quark model, since the data suggest that F 2 * « for Q 2 > 2 GeV 2 , contrarily to the result of 
the spectator quark model. These facts suggest that the valence quark and meson cloud contributions have opposite 
signs and cancel in the sum 12641 [265 1 . The direct consequence of the result for F 2 is the proportionality 
between the amplitudes A\/% and 6*1/2 for Q 2 > 2 GeV 2 12661 . 

Other applications of the covariant spectator quark model are the elastic electromagnetic form factors of the 
baryon octet (spin 1/2) [258, 267], and the baryon decuplet (spin 3/2) [255, 268-270], as well as the electromag- 
netic transition between octet and decuplet baryons, similarly to the jN —> A(1232) reaction 127 1 1 . The study of 
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FIG. 42. Dirac-type form factors F* for — > N* transitions. Left panel: jN — >• Pn(1440) reaction | 263 |. Right panel: 
7 7V -> 5xi (1535) reaction 12641 . 

the octet electromagnetic structure in the nuclear medium is also in progress \212\. 

Future work will establish how higher angular momentum states in the wave function, namely P and D states, 
may contribute to the nucleon form factors. This work will be facilitated by the results in Ref. B273I . where it was 
already possible to constrain those terms of the wave function by existing deep inelastic scattering data. 

Extensions for higher resonances are underway for Pn(1710), Z?i3(1520) and Sii(1650). The last two cases 
depend on the inclusion of an isospin 1/2, spin 3/2 core in a state of the total angular momentum 1/2. These states 
are expected to be the same as that in the part of the nucleon structure 112731 . 

In future developments, the quality and quantity of the future lattice QCD studies will be crucial to constrain the 
parameterization of the wave functions, and clarify the effect of the valence quarks and meson cloud, following 
the successful applications to the lattice QCD regime for the nucleon 12581 12731 . 7./V — > A(1232) transition 025711 
and Roper 112631 . 

In parallel, the comparison with the estimate of the quark core contributions performed by the EBAC group 
preferentially for Q 2 > 2 GeV 2 fl62"ll70l . will be also very useful in the next two years. To complement the quark 
models, the use of dynamical models and/or effective chiral models M274I to estimate the meson cloud effects 
are also very important. This is particularly relevant for the 7iV — > iV(1535) reaction. From the experimental 
side, new accurate measurements in the low Q 2 region as well as the high Q 2 region, as will be measured in 
the future after the Jefferson Lab 12 GeV upgrade, will be crucial, for the purposes of either to test the present 
parameterizations at high Q 2 , or to calibrate the models for new calculations at even larger Q 2 . The clarifications 
between the different analysis of the data such as EBAC, CLAS, SAID, MAID, Jiilich and Bonn-Gatchina, will 
also have an important role 1131 [T7l l39l l40l I275ti277l . 



VII.C. Nucleon electromagnetic form factors and electroexcitation of low lying nucleon resonances up to 
Q 2 = 12 GeV 2 in a light-front relativistic quark model 

VII.C. 1. Introduction 

In recent decade, with the advent of the new generation of electron beam facilities, there is dramatic progress in 
the investigation of the electroexcitation of nucleon resonances with significant extension of the range of Q 2 . The 
most accurate and complete information has been obtained for the electroexcitation amplitudes of the four lowest 
excited states, which have been measured in a range of Q 2 up to 8 and 4.5 GeV 2 for the A(1232)i-33, N (1535)5n 
and iV(1440)Pn, N(1520)D 13 , respectively (see reviews HUHIT1 and the recent update |[T8l[T38l . At relatively 
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small Q 2 , nearly massless Goldstone bosons (pions) can produce significant pion-loop contributions. However it is 
expected that the corresponding hadronic component, including meson-cloud contributions, will be losing strength 
with increasing Q 2 . The Jefferson Lab 12 GeV upgrade will open up a new era in the exploration of excited 
nucleons when the ground state and excited nucleon's quark core will be fully exposed to the electromagnetic 
probe. 

Our goal is to predict 3q core contribution to the electroexcitation amplitudes of the resonances A(1232)P33, 
iV(1440)Pn, 7V(1520)-Di3, and AT(1535)5n. The approach we use is based on light-front (LF) dynamics which 
realizes Poincare invariance and the description of the vertices N(N*) — »• 3q,Nw in terms of wave functions. 
The corresponding LF relativistic model for bound states is formulated in Refs. [278 -281 ]. The parameters of the 
model for the 3q contribution have been specified via description of the nucleon electromagnetic form factors in the 
approach that combines 3q and pion-cloud contributions. The pion-loop contributions to nucleon electromagnetic 
form factors have been described according to the LF approach of Ref. M282I . 



VII.C.2. Quark core contribution to transition amplitudes 

The 3q contribution to the 7* TV — > N(N*) transitions has been evaluated within the approach of Refs. [280 
281 1 where the LF relativistic quark model is formulated in infinite momentum frame (IMF). The IMF is chosen 
in such a way, that the initial hadron moves along the z-axis with the momentum P —> 00, the virtual photon 

/ 2 _ 2 r-)2 2,2 r-j2 \ 

momentum is fc M = I m ° ut ^p" , Q ± , — ™p n ) , the final hadron momentum is P' = P + k, and 

Q 2 = —k 2 = Q^; rai n and m out are masses of the initial and final hadrons, respectively. The matrix elements of 
the electromagnetic current are related to the 3g-wave functions in the following way: 

~ < N(N*),S' Z \J°£\N,S Z > - eSi J *'+Q 4 *dr, (25) 

where S z and S' z are the projections of the hadron spins on the z-direction, Qi (i = a, b, c) are the charges of the 
quarks in units of e, e 2 /47r = a, W and are wave functions in the vertices N(N*) — > 3q, and dT is the phase 
space volume: 

dr = (27r) -6^qcX^^ c (26) 

The quark momenta in the initial and final hadrons are parameterized via: 

p t = x t V + qu . , p' t = Xl V + , (27) 
Pq i± =P'q^=0, Sq i± = Sq^ = 0, q! i± = q iX - iftQ ±1 (28) 
= 1, y a =x a -l, y b = x b , y c = x c . (29) 

Here we have supposed that quark a is an active quark. 

The wave function is related to the wave function in the c.m.s. of the system of three quarks through Melosh 
matrices |283|: 

* = U + (p a )U + (p b )U + (p c )y fs M<la,<lb,<lc), (30) 

where we have separated the flavor-spin-space part of the wave function f ss in the c.m.s. of the quarks and its 
spatial part $(q a , qb, q c )- The Melosh matrices are defined by 

tti \ m q + M x l + iei m aiq im 
V(m„ + M oXl ) 2 +q^ 

where m q is the quark mass. The flavor-spin-space parts of the wave functions are constructed according to 
commonly used rules [25 , 284]. To construct these parts we need also the z-components of quark momenta in the 
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c.m.s. of quarks. They are defined by: 

q - = 2 [ XiMo - ) ' = 2 ~ J ' ( } 

where M and Mq are invariant masses of the systems of initial and final quarks: 



„ q 2 , + m 2 )2 q'ji + m 2 



To study sensitivity to the form of the quark wave function, we employ two widely used forms of the spatial 
parts of wave functions: 

$! ~ «cp(-M 2 /6af), $ 2 ~ ezp [-(q 2 + q 2 + q 2 )/2a 2 ] , (34) 
used, respectively, in Refs. j|278ti28ll and fl56l . 

VII. C. 3. Nucleoli 

The nucleon electromagnetic form factors were described by combining the 3q-core and pion-cloud contribu- 
tions to the nucleon wave function. With the pion loops evaluated according to Ref. |282|, the nucleon wave 
function has the form: 

\N >= 0.95|3<7 > +0.313|7Vtt >, (35) 

where the portions of different contributions were found from the condition the charge of the proton be unity: 
F\ P (Q 2 = 0) = 1. The value of the quark mass at Q 2 = has been taken equal to m q (0) = 0.22 GeV from the 
description of baryon and meson masses in the relativized quark model 11561 12851. Therefore, the only unknown 
parameters in the description of the 3q contribution to nucleon formfactors were the quantities a,i and a 2 in Eqs. 



(34 1. These parameters were found equal to 

at = 0.37 GeV, a 2 = 0.405 GeV (36) 



from the description of the magnetic moments at Q 2 = (see Fig. 43 I. The parameters (36 1 give very close 
magnitudes for the mean values of invariant masses and momenta of quarks at Q 2 = 0: < Mq >« 1.35 GeV 2 
and < q 2 0.1 GeV 2 , i = a, b.c. 

The constant value of the quark mass gives rapidly decreasing form factors Ge p (Q 2 ), Gmp(Q 2 ), and Gmu(Q 2 ) 
(see Fig. 43 1. The wave functions (34i increase as m q decreases. Therefore, to describe the experimental data we 



have assumed the Q -dependent quark mass that decreases with increasing Q . 

TO ( 1)(Q2)= 0-22 GeV (2) 2 0.22 GeV 

« WJ l + Q 2 /60GcV 2 ' 9 W ' l + Q 2 /10GeV 2 { ' 

for the wave functions $i and $2, respectively. Momentum dependent quark mass allowed us to obtain good 
description of the nucleon electromagnetic form factors up to Q 2 = 16 GeV 2 . From Fig. 44 it can be seen that at 
Q 2 > 2 GeV 2 , these form factors are dominated by the 3q-core contribution. 



VII.C.4. Nucleon resonances A(1232)P 33 , JV(1440)Pn, iV(1520)L>i 3 , and iV(1535)Si 



No calculations are available that allow for the separation of the 3q and Nn (or nucleon-meson) contributions 
to nucleon resonances. Therefore, the weights c* (c* < 1) of the 3q contributions to the resonances: \N* >= 
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FIG. 43. Nucleoli electromagnetic form factors. The solid curves correspond to the results obtained taking into account two 
contributions to the nucleon (Eq. \35\ : the pion-cloud |282| and the 3q core with the running quark masses < |37[ > for the wave 
functions $1 (black curves) and $2 (red curves) in Eqs. |34b . The black and red dashed curves are the results obtained for 
the nucleon taken as a pure 3q state with the parameters \36\ and constant quark mass. Dotted curve for Geti(Q 2 ) is the pion 
cloud contribution 12821 . Data are from Refs. I286U294I . 



c* \3q > +..., are unknown. We determine these weights by fitting to experimental amplitudes at Q 2 = 2 — 3 GeV 2 , 
assuming that at these Q 2 the transition amplitudes are dominated by the 3g-core contribution, as is the case for 
the nucleon. Then we predict the transition amplitudes at higher Q 2 (see Figs. 45p8 1. 



As it is shown in Refs. 1307113081 . there are difficulties in the utilization of the LF approaches 11251 12781 12791 
I309| |310I for the hadrons with spins J > 1. These difficulties can be avoided if Eq. ( 25 1 is used to calculate only 
those matrix elements that correspond to S' z = J 13071 . This restricts the number of transition form factors that 
can be calculated for the resonances A(1232)P33 and iV(1520)£>i3, and only two transition form factors can be 
investigated for these resonances: Gi(Q 2 ) and (^(Q 2 ) (the definitions can be found in review [10]). For these 
resonances we can not present the results for the transition helicity amplitudes. The results for the resonances with 
J = |: AT(1440)Pn and AT(1535)5n, are presented in terms of the transition helicity amplitudes. 
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FIG. 44. Nucleon electromagnetic form factors. The legend for the black and red solid curves is as for Fig. [43] Dotted curves 
are the pion cloud contributions | 282 |. 



VII. C. 5. Discussion 



The important feature of the obtained predictions for the resonances is the fact that at Q 2 > 2 — 3 GeV 2 both 
investigated amplitudes for each resonance are described well by the 3q contribution by fitting the only parameter, 
that is the weight of this contribution to the resonance. These predictions need to be checked at higher Q 2 . 

The results for the resonances allow us also to make conclusions on the size and form of expected pion-cloud 
and/or meson-baryon contributions to the amplitudes. According to our predictions for the 3q contributions, one 
can expect that pion-cloud contributions to the form factor G2{Q 2 ) for the A(1232)P33, to Si/ 2 amplitude for the 
iV(1440)Pn, and to the form factor G\(Q 2 ) for the 7V(1520)-Di3 are small. Large contributions are expected to 
the longitudinal amplitude for the ^(1535)511 and to the form factor G 2 {Q 2 ) for the iV(1520)Di3. The expected 
pion-cloud contributions to the form factor G\ (Q 2 ) for the A(1232)P33 and to A 1 / 2 amplitude for the ^(lSSS)^!! 

*Mp(n) 



have Q 2 behavior similar to that in the nucleon formfactors G Mp(n)(Q 2 )- In Fig. 46 by dotted curves we show 



estimated pion-cloud contribution to Ai/ 2 amplitude for the Roper resonance. It can be seen that non-trivial Q 2 - 
dependence of this contribution can be expected. 

The remarkable feature that follow from the description of the nucleon electromagnetic formfactors in our ap- 
proach is the decreasing quark mass with increasing Q 2 . This is in qualitative agreement with the QCD lattice 
calculations and with Dyson-Schwinger equations [28 29l 131 11 where the running quark mass is generated dy- 
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FIG. 45. The j*p -> A(1232)P 33 transition form factors; Gi (Q 2 ) ~ G M - G B . Weight factors are c%\ = 0.67 ± 0.04 and 

(2) r i 

c N , — 0.72 ± 0.04 for the wave functions $i (black curves) and $2 (red curves) in Eqs. 1 34 1. Solid circles correspond to the 
amplitudes extracted from the CLAS data by JLab group [ 17 1, bands represent model uncertainties of these results. The results 
from other experiments are: open triangles [295 -297 1; open crosses 1 298 - 300 1; open rhombuses 1 301 1; open boxes 13021 ; and 
open circles I303II304I . 



namically. The mechanism that generates the running quark mass can generate also the quark anomalous magnetic 
moments and form factors. This should be incorporated in model calculations. Introducing quark form factors will 
cause a faster Q 2 fall-off of electromagnetic form factors in quark models. This will force m q (Q 2 ) to drop faster 
with Q 2 to describe the data. 



VII.D. Light-Front Holographic QCD 

The relation between the hadronic short-distance constituent quark and gluon particle limit and the long-range 
confining domain is yet one of the most challenging aspects of particle physics due to the strong coupling nature 
of Quantum ChromoDynamics, the fundamental theory of the strong interactions. The central question is how one 
can compute hadronic properties from first principles; i.e., directly from the QCD Lagrangian. The most successful 
theoretical approach thus far has been to quantize QCD on discrete lattices in Euclidean space-time. 113 1 211 Lattice 
numerical results follow from computation of frame-dependent moments of distributions in Euclidean space and 
dynamical observables in Minkowski space-time, such as the time-like hadronic form factors, are not amenable 
to Euclidean lattice computations. The Dyson-Schwinger methods have led to many important insights, such as 
the infrared fixed point behavior of the strong coupling constant, [313] but in practice, the analyses are limited to 
ladder approximation in Landau gauge. Baryon spectroscopy and the excitation dynamics of nucleon resonances 
encoded in the nucleon transition form factors can provide fundamental insight into the strong-coupling dynamics 
of QCD. New theoretical tools are thus of primary interest for the interpretation of the results expected at the new 
mass scale and kinematic regions accessible to the JLab 12 GeV Upgrade Project. 

The AdS/CFT correspondence between gravity or string theory on a higher-dimensional anti-de Sitter (AdS) 
space and conformal field theories in physical space-time 113 141 has led to a semi-classical approximation for 
strongly-coupled QCD, which provides physical insights into its non-perturbative dynamics. The correspondence 
is holographic in the sense that it determines a duality between theories in different number of space-time di- 
mensions. This geometric approach leads in fact to a simple analytical and phenomenologically compelling non- 
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FIG. 46. The y*p -> iV(1440)Pu transition amplitudes. Blue lines correspond to the MAID results 13051 13061 . Dotted 
curves are estimated pion-cloud contributions, c^l — 0.73 ± 0.05, c^l — 0.77 ± 0.05. The open triangles correspond to the 
amplitudes extracted from CLAS 2n electroproduction data |38 |. Other legend is as for Fig. |45| 




perturbative approximation to the full light-front QCD Hamiltonian - "Light-Front Holography". [ 3 1 5 1 Light-Front 
Holography is in fact one of the most remarkable features of the AdS/CFT correspondence. 13141 The Hamiltonian 
equation of motion in the light-front (LF) is frame independent and has a structure similar to eigenmode equations 
in AdS space. This makes a direct connection of QCD with AdS/CFT methods possible. [ 3 1 5 1 Remarkably, the 
AdS equations correspond to the kinetic energy terms of the partons inside a hadron, whereas the interaction terms 
build confinement and correspond to the truncation of AdS space in an effective dual gravity approximation. 031511 
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FIG. 48. The 7*p — > 7V(1535)Sn transition amplitudes. The amplitudes extracted from the CLAS and JLab/Hall C data 
on ep — > er/p are: the stars |52|, the open boxes 1 14|, the open circles 1531 . the crosses |54|, and the rhombuses 1171 1371 . 
c%l = 0.88 ± 0.03, c$ = 0.94 ± 0.03. Other legend is as for Fig. ^ 



One can also study the gauge/gravity duality starting from the bound-state structure of hadrons in QCD quantized 
in the light-front. The LF Lorentz-invariant Hamiltonian equation for the relativistic bound-state system is 

P M P^(P)) = (P+P--P 2 X ) \^(P))=M 2 \HP)),P ± =P°±P 3 , (38) 

where the LF time evolution operator P~ is determined canonically from the QCD Lagrangian. 0161 To a first 
semi-classical approximation, where quantum loops and quark masses are not included, this leads to a LF Hamil- 
tonian equation which describes the bound-state dynamics of light hadrons in terms of an invariant impact variable 
( 113 15 1 which measures the separation of the partons within the hadron at equal light-front time r = x° + x 3 . 113 1711 
This allows us to identify the holographic variable z in AdS space with an impact variable £. 113151 The resulting 
Lorentz-invariant Schrodinger equation for general spin incorporates color confinement and is systematically im- 
provable. 

Light-front holographic methods were originally introduced H 3 1 8 1 1 3 1 9l by matching the electromagnetic current 
matrix elements in AdS space [320] with the corresponding expression using LF theory in physical space time. It 
was also shown that one obtains identical holographic mapping using the matrix elements of the energy-momentum 
tensor [321 ] by perturbing the AdS metric around its static solution. 032211 

A gravity dual to QCD is not known, but the mechanisms of confinement can be incorporated in the gauge/gravity 
correspondence by modifying the AdS geometry in the large infrared (IR) domain z ~ 1/Aqcd> which also sets 
the scale of the strong interactions. 1132311 In this simplified approach we consider the propagation of hadronic 
modes in a fixed effective gravitational background asymptotic to AdS space, which encodes salient properties of 
the QCD dual theory, such as the ultraviolet (UV) conformal limit at the AdS boundary, as well as modifications 
of the background geometry in the large z IR region to describe confinement. The modified theory generates the 
point-like hard behavior expected from QCD, H3241 13251 instead of the soft behavior characteristic of extended 
objects. 13231 
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VII.D.l. Nucleon Form Factors 



In the higher dimensional gravity theory, hadronic amplitudes for the transition A —> B correspond to the 
coupling of an external electromagnetic (EM) field A M (x, z) propagating in AdS space with a fermionic mode 
^p(x, z) given by the left-hand side of the equation below 

d 4 x dz y/g y B ,p> 0, z) &m T A Af (x, z)^ A ,p(x, z) ~ 

(P'-P- q) e^j B (P>), a'\J^ A (P),a), 

where the coordinates of AdSs are the Minkowski coordinates x^ and z labeled x M — [x^,z), with M,N = 
1, • • • 5, gis the determinant of the metric tensor and efj is the vielbein with tangent indices A, B = 1, • • ■ ,5. The 
expression on the right-hand side represents the QCD EM transition amplitude in physical space-time. It is the EM 
matrix element of the quark current J M = e q q-j^q, and represents a local coupling to pointlike constituents. Can 
the transition amplitudes be related for arbitrary values of the momentum transfer ql How can we recover hard 
pointlike scattering at large q from the soft collision of extended objects? |320| Although the expressions for the 
transition amplitudes look very different, one can show that a precise mapping of the J + elements can be carried 
out at fixed LF time, providing an exact correspondence between the holographic variable z and the LF impact 
variable £ in ordinary space-time. 113 181 

A particularly interesting model is the "soft wall" model of Ref. [ 326 1, since it leads to linear Regge trajectories 
consistent with the light-quark hadron spectroscopy and avoids the ambiguities in the choice of boundary con- 
ditions at the infrared wall. In this case the effective potential takes the form of a harmonic oscillator confining 
potential k 4 z 2 . For a hadronic state with twist t = N + L (N is the number of components and L the internal 
orbital angular momentum) the elastic form factor is expressed as a t — 1 product of poles along the vector meson 
Regge radial trajectory (Q 2 — —q 2 > 0) B319I1 

F(Q 2 ) = s w 1 x 7 o (39) 

where M 2 — > An 2 (n + 1/2) . For a pion, for example, the lowest Fock state - the valence state - is a twist-2 state, 



and thus the form factor is the well known monopole form. The remarkable analytical form of Eq. ( 39 1, expressed 
in terms of the p vector meson mass and its radial excitations, incorporates the correct scaling behavior from the 
constituent's hard scattering with the photon [324 325 1 and the mass gap from confinement. 



VII.D.2. Computing Nucleoli Form Factors in Light-Front Holographic QCD 

As an illustrative example we consider in this section the spin non-flip elastic proton form factor and the form 
factor for the j*p —> A r (1440)Pn transition measured recently at JLab. In order to compute the separate features 
of the proton an neutron form factors one needs to incorporate the spin-flavor structure of the nucleons, prop- 
erties which are absent in the usual models of the gauge/gravity correspondence. This can be readily included 
in AdS/QCD by weighting the different Fock-state components by the charges and spin-projections of the quark 
constituents; e.g., as given by the SU(6) spin-flavor symmetry. 

Using the SU(6) spin-flavor symmetry the expression for the spin-non flip proton form factors for the transition 
n, L n'L is 13311 

F?n,L^AQ 2 )= R4 J %^ L (z)V{Q,z)^\z), (40) 
where we have factored out the plane wave dependence of the AdS fields 

- ^(^tW^ 7 '^ 1 ^ 2 ' ^'*' /2 - ,41) 
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FIG. 49. Dirac proton form factors in light-front holographic QCD. Left: scaling of proton elastic form factor Q 4 F-f (Q 2 ). 
Right: proton transition form factor Ff N ^ N * (Q 2 ) for the j*p — > iV(1440)Pn transition. Data compilation from Diehl 1 327 1 
(left) and CLAS tt and 2tt electroproduction data |T7||38l|328 329 1 (right). 




The bulk-toboundary propagator V(Q, z) has the integral representation 0321 



V(Q,z) = k 2 z 2 



dx 



{1-xf 



-k 2 z 2 x/ (1— a;) 



(42) 



with V(Q = 0, z) = V(Q, z = 0) = 1. The orthonormality of the Laguerre polynomials in (41 1 implies that the 
nucleon form factor at Q 2 = is one if n = n' and zero otherwise. Using ( 42 1 in ( 40 1 we find 
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for the elastic proton Dirac form factor and 



F? N ^(Q 2 ) = 4 7 „,w —> (44 > 



S 14 H| 



for the EM spin non-flip proton to Roper transition form factor. The results (|43| and (04]), compared with available 



data in Fig. 49 correspond to the valence approximation. The transition form factor d44b is expressed in terms of 
the mass of the p vector meson and its first two radial excited states, with no additional parameters. The results in 
Fig. 49 are in good agreement with experimental data. The transition form factor to the TV ( 1440) Pi 1 state shown 



in Fig. 49] corresponds to the first radial excitation of the three-quark ground state of the nucleon. In fact, the 
Roper resonance A r (1440)Pn and the iV(1710)Pn are well accounted in the light-front holographic framework 
as the first and second radial states of the nucleon family, likewise the A(1600)p33 corresponds to the first radial 
excitation of the A family as shown in Fig.[50]for the positive-parity light-baryons. 13331 In the case of massless 
quarks, the nucleon eigenstates have Fock components with different orbital angular momentum, L — and 
L = 1, but with equal probability. In effect, in AdS/QCD the nucleons angular momentum is carried by quark 
orbital angular momentum since soft gluons do not appear as quanta in the proton. 

Light-front holographic QCD methods have also been used to obtain general parton distributions (GPDs) in 
Ref. J334|, and a study of the EM nucleon to A transition form factors and nucleon to the Sn (1535) negative parity 
nucleon state has been carried out in the framework of the Sakai and Sugimoto model in Refs. |335| and [336 1 
respectively. It is certainly worth to extend the simple computations described here and perform a systematic study 
of the different transition form factors measured at JLab. This study will help to discriminate among models and 
compare with the new results expected from the JLab 12 GeV Upgrade Project, in particular at photon virtualities 
Q 2 > 5 GeV 2 , which correspond to the experimental coverage of the CLAS12 detector. 



VILE. Constituent Quark Models and the interpretation of the nucleon form factors 



Various Constituent Quark Models (CQM) have been proposed in the past decades after the pioneering work 
of Isgur and Karl (IK) 1155L Among them let us quote the relativized Capstick-Isgur model (CI) [ 156], the alge- 
braic approach (BIL) H337L the hypercentral CQM (hCQM) 113381 . the chiral Goldstone Boson Exchange model 
(xCQM) 115811 and the Bonn instanton model (BN) [339 340 1. They are all able to fairly reproduce the baryon 
spectrum, which is the first test to be performed before applying any model to the description of other baryon 
properties. The models, although different, have a simple general structure, since, according to the prescription 
provided by the early Lattice QCD calculations H341I . the three-quark interaction V% q is split into a spin-flavor in- 
dependent part Vi nv , which is 5'J7(6)-invariant and contains the confinement interaction, and a S'J7(6)-dependent 
part V s f, which contains spin and eventually flavor dependent interactions 



V 3q = V mv + V sf (45) 

In Tab. |IV] a summary of the main features of various Constituent Quark Models is reported. 

After having checked that these models provide a reasonable description of the baryon spectrum, they have been 
applied to the calculation of many baryon properties, including electrocouplings. One should however not forget 
that in many cases the calculations referred to as CQM calculations are actually performed using a simple h.o. wave 
function for the internal quark motion either in the non relativistic (HO) or relativistic (relHO) framework. The 
former (HO) applies to the calculations of refs. 134211 and 12841 . while the latter (relHO) is valid for ref. E5l . The 
relativized model (CI) of ref. |156] is used for a systematic calculation of the transition amplitudes in ref. [343] 
and, within a light front approach in refs. 13441 and ref. 11291 for the transitions to the A and Roper resonances 
respectively. In the algebraic approach [337], a particular form of the charge distribution along the string is assumed 
and used for both the elastic and transition form factors; the elastic form factors are fairly well reproduced, but 
there are problems with the transition amplitudes, specially at low Q 2 . There is no helicity amplitude calculation 
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TABLE IV. Illustration of the features of various CQMs 



CQM 


Kin. Energy 




Vsf 


ref. 


Isgur-Karl 


non rel. 


h.o. + shift 


OGE 


EST 


Capstick-Isgur 


rel 


string +coul-like 


OGE 


1 156 1 


U(7)B.I.L. 


relM 2 


vibr + L 


Gursey-Rad 


|337| 


Hypercentral G.S. 


non rel/ rel 


0(6): lin + hyp.coul 


OGE 


1 338 1 


Glozman-Riska 


rel 


h.o. / linear 


GBE 


|158| 


Bonn 


rel 


linear + 3 body 


instanton 


|339| 



with the GBE model, whereas the BN model has been also used for the helicity amplitudes |345|, with particular 
attention to the strange baryons M346L Finally, the hCQM has produced predictions for the transverse excitation 
of the negative parity resonances l27ll and also for the main resonances, both for the longitudinal and transverse 
excitation [ 347 1 . 

In some recent approaches the CQ idea is used to derive relations between the various electromagnetic form 
factors, relations which, after having fitted one selected quantity, say the elastic proton form factor (Sec. 2) or the 
helicity amplitude at intermediate Q 2 (Sec. 3), are used to predict the other quantities of interest. A remarkable 
prediction of both the proton elastic form factor and the proton transition to the Roper resonance is provided by 
the light-front holographic approach (Sec. 4). 



200 

150 

sC" 100 

I 50 
o 

C 
o 

S -50 

Q. 

CO 

5 -100 

-150 
-200 



h.o. A 



hCQM A 1/2 
hCQM A 3/2 




2 3 
Q 2 (GeV 2 ) 



FIG. 51. The Di3(1520) helicity amplitudes (upper part) and A1/2 (lower part) predicted by the hCQM (full curves), in 
comparison with the data of refs. [17 18 38 1 and the PDG values 1221 at the photon point. The h.o.results for two different 
values of the proton r.m.s. radius (0.5 fm and 0.86/m)A are also shown. 



The works briefly illustrated above have shown that the three-quark idea is able to fairly reproduce a large variety 
of observables, in particular the helicity amplitudes at medium Q 2 , however, a detailed comparison with data shows 
that, besides the fundamental valence quarks, other issues are or presumably will be of relevant importance for the 
interpretation of the transition amplitudes. These issues are: relativity, meson cloud and quark-antiquark pair 
effects, and quark form factors. 

A consistent relativistic treatment is certainly important for the description of the elastic nucleon form factors. 
In fact, in the non-relativistic hCQM |338|, the proton radius compatible with the spectrum is too low, about 
0.5 fm, and the resulting form factors [348 1 are higher than data. However, the introduction of the Lorentz boosts 
improves the description of the elastic form factors 13481 and determines a ratio ^ p G^/G p M lower than 1 H349I . 
Using a relativistic formulation of the hCQM in the Point Form approach, in which again the unknown parameters 
are fitted to the spectrum, the predicted elastic nucleon form factors are nicely close to data [350]. Furthermore, if 
one introduces quark form factors, an accurate description of data is achieved [ 350 1 . Since such form factors are 
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fitted, this means that they contain, in an uncontrolled manner, all the missing contributions. 

Applying hCQM to the excitation of higher resonances demonstrated that the inclusion of relativity is less 
crucial, since the Lorentz boosts affect only slightly the helicity amplitudes [351 1. A quite different situation 
occurs for the excitation to the A, which is a spin-isospin excitation of the nucleon and as such it shares with the 
nucleon the spatial structure. In this case relativity is certainly important, however it does not seem to be sufficient 
even within LF approaches. In fact, the good results of the Rome group [344 1 are obtained introducing quark form 
factors, while in Sec. 2 the quark wave function fitted to the elastic nucleon form factor leads to a lack of strength 
at low Q 2 in the A excitation. In Sec. 3 a pion cloud term is present from the beginning in the nucleon form factor, 
nevertheless the transition to the A is too low at low Q 2 . 

Of course, the future data at high Q 2 will force, at least for consistency reasons, to use a relativistic approach 
also for the other resonances. 
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FIG. 52. The ratio F p at high Q calculated using the theoretical form factors of Ref. 1 350 1. 



At medium-low Q 2 the behavior of the helicity amplitudes is often described quite well, also in a non relativistic 
approach [347 1. An example is provided by Fig. 51 where the hCQM results are compared with the more recent 



Jlab data. In Fig. 51 there are also the h.o. results, which do not seem to be able to reproduce the data. The good 



agreement achieved by the hCQM has a dynamical origin. Let us remind that in hCQM the SU (6)-invariant part 



of the quark potential of Eq. (45 i is 



V™ QM =~l + ™ (46) 

(x = p 2 + X 2 is the hyperradius) however the main responsible of the medium-high Q 2 behavior of the helicity 
amplitudes is the hypercoulomb interaction — |. In fact, in the analytical version of hCQM presented in ref. 13521 . 
it is shown that the helicity amplitudes provided by the — - term are quite similar to the ones calculated with the 
fullhCQM. 

The main problem with the description provided by CQM (non relativistic or relativistic) is the lack of strength 
at low Q 2 , which is attributed, with general consensus, to the missing meson cloud or quark-antiquark pair effects 
lf]T1l27ll57ll62l . In fact, it has been shown within a dynamical model lfTTl l57 62] that the meson-cloud contributions 
are relevant at low Q 2 and tend to compensate the lack of strength of unquenched three-quark models 13531 . 

To conclude, a fully relativistic and unquenched hCQM is not yet available and work is now in progress in this 
direction, but certainly it will be a valuable tool for the interpretation of the helicity amplitudes at high Q 2 . 

However, also taking into account the one pion contribution there seems to be some problem. In Sec. 2, the quark 
wave function is chosen in order to reproduce the proton form factor, in this way all possible extra contributions 
(meson cloud, quark form factors,....) are implicitly included, but the description of the N — A transition needs an 
extra pion term. On the other hand, in Sec. 3 it is shown that the pion term explicitly included in the fit to the proton 
is not sufficient for the description of the N — A transition. In fact, the inclusion of a pion cloud term, either fitted 
or calculated (as e.g. in ref. [354]) seems to be too restrictive, since it is equivalent to only one quark-antiquark 
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configuration. If one wants to include consistently all quark-antiquark effects, one has to proceed to unquenching 
the CQM, as it has been done in M355I . Such an unquenching is achieved by summing over all quark loops, that is 
over all intermediate meson-baryon states; the sum is in particular necessary in order to preserve the OZI rule. 

This unquenching has been recently performed also for the baryon sector [ 356 1 . The state for a baryon A is 
written as 



1*4 >= TV 



\A > + f & \BCklJ > 

BClj J 



< BCklJ\T^\A > 
M A ~E B - E c 



(47) 



where B (C) is any intermediate baryon (meson), Eg(Ec) are the corresponding energies, M A is the baryon mass, 
T* is the 3 Pq pair creation operator and k, I and J are the relative momentum, the orbital and total angular momen- 
tum, respectively. Such unquenched model, with the inclusion of the quark-antiquark pair creation mechanism, will 
allow to build up a consistent description of all the baryon properties (spectrum, form factors,...). There are already 
some applications 113561 . in particular it has been checked that, thanks to the summation over all the intermediate 
states prescribed in Eq. (??), the good account of the baryon magnetic moments provided by the standard CQM is 
not verified M356I . Using an interaction containing the quark-antiquark production the resonances acquire a finite 
width, at variance with what happens in all CQMs, allowing a consistent description of both electromagnetic and 
strong vertices. 

The structure of the state in Eq. (??) is more general than the one containing a single pion contribution. The 
influence of the quark-antiquark cloud will be certainly important at low Q 2 , but one can also expect that the 
multiquark components, which are mixed with the standard 3q states as in Eq. (??), may have a quite different 
behavior M3571 4359 1 in the medium-high Q 2 region, leading therefore to some new and interesting behavior also at 
short distances. Actually, there are some clues that this may really happen. First, it has been shown in [ 17 1 that the 
quantity Q 3 A^ 2 seems to become flat in the range around 4 GeV 2 (see Fig. |2j, while the CQM calculations do 
not show any structure. 

A second important issue is the ratio R p = [L v G p e jG v M between the proton form factors. A convenient way of 
understanding its behavior is to consider the ratio Q 2 F p / Ff , which is expected to saturate at high Q 2 I360l l36l1 . 
while it should pass through the value AM 2 /k p in correspondence of a zero for R p M362L The predictions of the 
hCQM 0501 are compared with the Jlab data 12861 12871 [363ti366l in Fig. [52] For a pure three-quark state, even 
in presence of quark form factors as in [ 350 1 the occurrence of a zero seems to be difficult, while an interference 
between three- and multi -quark configurations may be a possible candidate for the generation of a dip in the electric 
form factors B362L 

It is interesting to note that the Interacting Quark Diquark model introduced in Ref. ([367]) and its relativistic 
reformulation ([368 1), both of which do not exhibit missing states in the non strange sector under 2GeV 2 , give 
rise to a ratio i? = /i p G p e / G p m that goes through a zero at around 8GeV 2 after the introduction of quark form 
factors, as calculated in Ref. (|369))- 

Once the quark-antiquark pair creation effects have been included consistently in the CQM, it will be possible to 
disentangle the quark forms factors from the other dynamic mechanisms. The presence of structures with a finite 
dimension has been shown in a recent analysis of deep inelastic electron-proton scattering [370|. 
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